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SECTION  1 


INTRODUCTION  AND  SUMMARY 


l.L.  Introduction 


This  interim  report  is  submitted  in  partial  fulfillment  of  documen¬ 
tation  requirements  for  the  study:  Active  Control  of  Space  Structures-Six 
(ACOSS  Six).  For  the  successful  control  of  large  space  structures  (LSS) ,  it 
is  recognized  by  the  LSS  community  at  large  that  the  fundamental  problem  is 
the  design  of  a  finite-dimensional  compensator  to  control  an  infinite 
dimensional  system  (1).  To  date,  numerous  theoretical  contributions  towards 
"solving"  this  problem  have  been  made;  however,  the  resolution  of  this 
problem  requires  further  theoretical  research  which  must  be  validated  through 
appropriate  design,  analysis  and  experimentation.  As  a  result  of  our  endeavors, 
together  with  exposure  to  a  broad  spectrum  of  knowledge  provided  by  the  LSS 
community,  it  is  our  judgment  that  the  technical  issues  in  LSS  control 
technology  include  the  following: 

(1)  LSS  Modeling  accuracy  should  be  known  to  within  some  specified 
bounds.  Modeling  errors  will  limit  achievable  control  system 
performance.  These  errors  may  be  introduced  through  initially 
assumed  structural  properties  or  the  truncation  process  implicit 
in  the  finite-element  method.  In  space,  LSS  parameters  may 
vary  as  a  function  of  thermal  gradients,  configuration  changes, 
or  depletion  of  consumables.  The  more  stringent  the  mission 
performance  requirements,  the  greater  the  LSS  model  fidelity 
required. 

(2)  Upper  Atmosphere  Models  must  be  improved  and  verified  by  appro¬ 
priate  experiments.  Accurate  knowledge  of  the  external  forces 
(e.g.,  earth  magnetic  and  gravitational  fields,  solar  wind  and 
radiation  pressure,  drag)  acting  upon  a  LSS  is  necessary  to 
satisfy  precision  control  requirements. 

(3)  System  Identification  is  necessary  for  the  purpose  of  LSS  struc¬ 
tural  model  verification.  Parameter  identification  techniques 
must  be  developed  such  that  modal  frequencies,  damping  ratios, 
and  mode  shapes  can  be  accurately  determined.  Consideration 
must  be  given  to  the  type  of  sensors,  onboard  r*.ocessing 
requirements,  data  reduction,  post-processing  requirements,  and 
control  law  configuration. 

(4)  Sensor  and  Actuator  specifications  must  be  determined  in  order 
to  assess  the  applicability  of  existing  hardware,  as  well  as  to 
provide  new  directions  in  research  and  development. 

(5)  Control  Law  Design  Methodology  must  address  the  following: 

(a)  The  model  reduction  process  which  reduces  the  high¬ 
dimensional  finite-element  model  to  a  lower  order  design 
model . 
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(b)  The  design  method  for  determining  reduced-order 
compensators . 

(c)  The  criteria  for  determining  overall  closed-loop 
system  stability. 

(d)  Direct  digital  design  methodologies  and  implementation 
techniques . 

(e)  Sensor  and  actuator  placement  techniques  which  yield 
maximum  observability  and  controllability. 

(f)  Sensor  and  actuator  dynamics. 

(g)  The  dynamic  interaction  of  the  attitude,  vibration,  and 
shape  control  laws. 

The  theoretical  research  documented  in  the  present  report  addresses  control 
law  design  methodologies  with  particular  emphasis  on  model  reduction,  design 
methods  for  determining  reduced-order  compensators,  and  criteria  for 
determining  overall  closed-loop  system  stability. 

1 . 2  Research  Scope 


1.2.1  Scope  of  Theoretical  Research 


The  enabling  technology  necessary  for  successful  active  control  of 
large  precision  space  structures  requires  mastery  of  at  least  the  following 
principal  topics: 

(1)  Reduction  of  high-order,  high-fidelity  structural  models  to  lower- 
order  models  suitable  for  the  design  of  feedback  controllers. 

(2)  Design  of  fixed-order  compensators  which  ensures  stability  for 
the  total  closed-loop  system. 

(3)  Analytical  techniques  for  determining  overall  closed-loop 
system  stability  (including  robustness). 

The  theoretical  results  documented  in  this  report  address  these  topics. 
Given  a  finite-element  modal  representation  of  a  LSS,  a  technique  which 
aggregates  the  vibration  modes  with  respect  to  actuator  and  sensor  influences 
has  been  developed  (Section  4) .  By  proper  selection  of  the  aggregation 
transformation,  this  reduced-model  may  be  made  to  approximate  the  open-loop 
performance  of  the  LSS  more  accurately  than  methods  utilizing  direct  modal 
truncation.  New  results  relative  to  the  extended  Kosut  design  approach  to 
multivariable  output  feedback  controller  design  are  presented  in  Section  2. 

The  stochastic  output  feedback  controller  (Section  7)  is  formulated  as  the 
zero-order  limit  of  the  minimum  variance  fixed-order  output  feedback 
compensator.  Section  5  discusses  pareto-optimal  control  (multiple-objective 
optimization)  with  respect  to  the  design  of  LSS  vibration  controllers. 
Stability  margins  for  observer-based,  fixed-order  compensators  is  investigated 
in  Section  3.  The  general  conditions  for  LSS  closed-loop  stabilization  by 
output  feedback  control  are  derived  (see  Section  6) . 


2 


1.2.2  Scope  of  Applications  Research 

In  order  to  assess  performance,  sensitivity,  and  hardware  requirements 
of  the  various  active  structural  control  methods,  a  simple  but  realistic 
evaluation  model  is  necessary.  Section  9  describes  such  a  model  (Draper 
Model  #2 )  which  is  based  on  realistic  sizes  and  masses.  This  model  contains 
a  simple  optical  system  with  unclassified  performance  measures  and  tolerances 
A  line-of-sight  (LOS)  performance  measure  (Section  8)  is  available  for 
evaluation  purposes.  Section  10  presents  preliminary  results  of  a  sinusoidal 
disturbance  (see  Section  9)  impact  on  the  LOS  performance. 

1 . 3  Summary  of  Principal  Results 

1.3.1  Results  from  Theoretical  Research 


Previously  reported  results  relating  to  the  algebraic  consistency  of, 
and  the  structure  of  solutions  to,  the  Kosut  gain  equation  are  substantially 
improved  (Section  2).  Conditions  for  design-model-independence  of  free 
parameters  occurring  in  such  solutions  are  given.  The  nature  and  extent 
of  performance  degradation  to  be  expected  due  to  decentralizing  the  feedback 
information  structure  is  obtained.  Stable  designs  obtained  by  the  extended 
Kosut  method  are  shown  to  be  robust.  Major  theoretical  problems  that  need  to 
be  addressed  are:  (1)  Characterization  of  systems  for  which  the  extended 
Kosut  approach  leads  to  a  stable  closed-loop  design  (at  the  design  model 
level);  (2)  Synthesize  a  systematic  procedure  for  assigning  values  to  the 
free  parameters  (w. «  ■  available)  so  as  to  improve  stability  and/or  perfor¬ 
mance  of  the  full-order  system  incorporating  the  reduced-order  controller. 

The  stability  margins  for  fixed-order  compensators  (Section  3)  was 
investigated.  The  observer-based  fixed-order  compensator  was  assumed  since 
it  is  representative  of  the  most  general  case.  The  objective  was  to  find 
a  procedure  which  would  determine  the  gain  matrices  such  that  the  full-order, 
closed-loop  system  remained  stable.  If  the  resulting  equations  could  be 
solved  analytically,  then  explicit  stability  conditions  could  be  given.  Due 
to  the  complexity  of  these  equations,  an  analytical  closed-form  expression  is 
not  apparent  at  this  time. 

A  standard  approach  to  LSS  model  reduction  is  to  truncate  the  open- 
loop  modal  representation  directly  by  neglecting  certain  modes.  The  aggregate 
model  development  (Section  4)  provides  a  technique  for  aggregating  the 
vibration  modes  (state  variables)  with  respect  to  actuator  and  sensor 
influences.  The  aggregation  is  accomplished  by  a  similarity  transformation 
which  can  be  implemented  by  Gaussian  elimination  or  Gram-Schmidt  orthogon- 
alization  procedures.  By  proper  choice  of  the  aggregation  transformation, 
the  resulting  model  may  be  reduced  to  a  lower  order  while  retaining  an 
approximate  open-loop  performance  equivalent  to  that  of  the  full-order  model. 
Since  the  aggregation  transformation  is  nonunique,  several  approaches  were 
examined  which  utilized  this  nonuniqueness  for  imposing  additional  properties 
on  the  aggregate  design  model.  The  feasibility  of  these  approaches  was  also 
analyzed . 
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Pareto  optimality  (Section  5)  is  a  practical  notion  of  "multiple- 
objective"  optimality  that  is  useful  t-~  .^uirol  designers  faced  with  multiple 
disparate  performance  obierf5’^^.  A  design  is  Pareto  optimal  if  none  of  the 
multiple  objectives  can  be  further  improved  without  degrading  any  other. 

The  ..^plication  of  Pareto-optimal  control  to  the  design  of  LSS  vibration 
controllers  provides  the  control  designer  with  a  systematic  procedure  for 
making  efficient  design  tradeoffs.  The  underlying  multiple-objective  design 
problems  relevant  to  LSS  active  control  are  shown  via  Canavin's  modal  dash- 
pot  concept  and  Aubrun's  low  authority  control  theory.  Favorable  results 
were  obtained  by  Pareto  optimal  designs  on  the  12-mode  Draper  tetrahedral 
model.  Efficient  tradeoffs  between  high  damping  and  low  output  feedback 
gains  were  demonstrated. 

General  conditions  for  closed-loop  LSS  stabilization  (Section  6)  by 
output  feedback  control  were  developed.  If  these  conditions  are  satisfied  by 
proper  design,  asymptotic  stability  of  the  closed-loop  LSS  system  is  ensured. 
The  resulting  theorems  are  applicable  to  the  general  case  of  LSS  that  contain 
not  only  stable  elastic  modes  but  also  rigid-body  modes  and  unstable  modes. 

A  specific  output  feedback  control  design  technique  is  also  presented  for 
integrated  control  of  elastic  and  rigid-body  modes  which  increases  both  the 
damping  and  frequency  of  selected  modes.  This  controller  is  closed-loop 
asymptotically  stable  as  well  as  robust  with  respect  to  modal  truncation  and 
parameter  variations.  A  technique  is  presented  for  synthesis  of  actuator  in¬ 
fluences  such  that  control  spillover  to  selected  modes  is  prevented.  It  is 
subsequently  shown  that  this  concept  is  readily  extendable  to  placement  of 
actuators  which  minimizes  control  spillover. 

The  minimum  variance  fixed  form  output  feedback  compensator  has  been 
proposed  as  a  possible  solution  to  the  problem  of  designing  reduced-order 
controllers  for  plants  of  large  dimensions.  The  zero-order  limit  of  this 
compensator  is  the  stochastic  output  feedback  controller  (Section  7) .  A 
detailed  discussion  of  the  stochastic  output  feedback  controller  and  its 
relation  to  several  other  compensator  formulations  such  as  the  minimum  order 
dual  observer  based  compensator  and  the  Levine-Athans-Johnson  reduced-order 
compensator  is  presented.  Lower  and  upper  bounds  on  the  optimal  stochastic 
output  feedback  controller  cost  are  derived.  These  bounds  should  be  useful 
in  preliminary  evaluations  of  subsequent  controller  performance.  These 
bounds  may  also  provide  useful  starting  estimates  in  iteratively  computing 
the  optimal  gains.  Presently,  this  computation  represents  the  major 
theoretical  problem  with  the  stochastic  output  feedback  controller. 

1.3.2  Results  from  Applications  Research 

The  Draper  model  #2  (Sections  8,  9)  was  evaluated  with  respect  to  the 
sinusoidal  disturbances  given  in  Section  9.  The  preliminary  results 
(Section  10)  show  that  the  addition  of  damping  to  individual  modes  is  not 
sufficient  to  suppress  a  sinusoidal  disturbance  impact  on  the  steady-state 
LOS  amplitude  error. 
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SECTION  2 


SUBOPT IMAL  OUTPUT  FEEDBACK  -  EXTENDED  KOSUT  METHOD 


2. 1  Introduction 

2.1.1  Synopsis  of  Previous  Research 

The  Kosut  approach  to  multivariable  output  feedback  controller  design 
was  originally  developed  [1]  as  an  alternative  to  optimal  output  feedback 
approaches  [2]  in  which  an  acceptable  solution  to  the  strongly  coupled 
nonlinear  equations  for  the  feedback  gains  was  seldom  obtainable  [3] .  By 
accepting  suboptimal  values  for  the  control  system  performance  measure, 
equations  for  the  feedback  gains  were  obtained  that  could  be  solved  without 
iteration.  Part  of  the  price  paid  for  the  resultant  simplicity  of  the 
design  procedure  was  the  lack  of  a  guarantee  for  stability  of  the  closed- 
loop  system.  As  a  result,  the  method  found  only  limited  application  [4], 

Recently,  a  fresh  look  at  the  Kosut  approach  was  taken  to  determine 
its  feasibility  for  use  as  a  design  tool  for  controllers  of  large  flexible 
spacecraft  [5] .  It  was  immediately  recognized  that  certain  assumptions  on 
which  the  method  is  based  are  incompatible  with  corresponding  properties 
of  reduced-order  modal  models  that  are  needed  for  controller  design  with 
large-scale  systems;  in  particular,  the  state-to-output  observation  matrix 
in  the  reduced-order  model  must  have  maximum  rank.  It  was  observed  that 
this  restriction  on  the  design  technique  is  not  in  fact  necessary.  The 
linear  equation  for  the  feedback  gain  matrix  can  be  solved  in  all  cases; 
in  particular,  when  this  equation  is  rank-deficient,  the  general  solution 
for  the  feedback  gains  contains  free  parameters  proportional  in  number  to 
the  rank-deficiency.  It  was  demonstrated  how  such  free  parameters  can  be 
selected  so  as  to  improve  the  performance  of  the  full-order  closed-loop 
system  relative  to  what  would  be  possible  if  the  rank  restriction  on  the 
observation  matrix  were  retained  [6] .  As  part  of  a  parallel  investigation, 
the  Kosut  approach  incorporating  the  above  extensions  was  compared  to  other 
representative  multivariable  controller  design  approaches  in  the  context 
of  a  specific  example  having  certain  of  the  principal  features  of  a  large 
flexible  spacecraft  model  [5,7].  In  spite  of  the  lack  of  a  guarantee  for 
stability  of  the  closed-loop  design  model,  the  performance  of  the  example 
system  incorporating  a  Kosut  controller  compared  very  favorably  with  that 
using  other  controller  designs.  Finally,  an  algorithm  for  the  numerical 
solution  of  rank-deficient  linear  equations  was  developed  [8]  to  enable 
application  of  the  extended  Kosut  approach  to  high-order  design  models. 

2.1.2  Summary  of  New  Results 

New  results  relative  to  the  extended  Kosut  design  approach  appearing 
in  the  present  report  are  briefly  summarized.  A  general  criterion  for 
stability  of  a  closed  loop  system  -  at  the  design  model  level  -  for  output 
feedback  systems  is  given  in  terms  of  the  basic  parameters  of  the  full-order 
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structural  model,  including  sensor  and  actuator  parameters,  and  the 
feedback  gain  matrix.  A  simple  and  useful,  but  undesirably  restrictive, 
sufficient  condition  for  stability  is  also  given.  Robustness  of  stability 
is  established  for  designs  obtained  by  the  extended  Kosut  method  in  the 
situation  that  only  rate  sensing  is  used.  In  addition,  certain  previously 
reported  results  are  improved.  In  particular,  conditions  are  given  under 
which  free  parameters  occuring  in  the  feedback  gain  matrix  do  not  influence 
the  stability  of  the  design  model.  Under  such  conditions,  a  decision  on 
their  selection  can  be  based  solely  on  considerations  of  stability  and 
performance  of  the  full-order  closed-loop  system.  Also,  a  simplified  proof 
is  given  of  the  general  algebraic  consistency  of  the  linear  equation  for 
the  feedback  gains  [6] . 

It  has  been  noted  previously  [5]  that  the  essential  nature  of  the 
Kosut  design  approach,  including  extensions,  is  not  fundamentally  altered 
when  a  decentralized  information  structure  is  imposed  in  the  feedback  path. 

A  closer  look  at  the  differences  in  detail,  however,  suggests  a  much  more 
satisfactory  way  of  representing  and  understanding  the  general  solution  to 
the  rank-deficient  gain  equation.  A  modified  representation  for  the  general 
solution  is  presented  that  is  more  readily  related  to  the  classical  algebraic 
theory  of  linear  spaces.  The  algorithm  for  numerical  solution  of  such 
equations  is  correspondingly  modified.  A  much  clearer  understanding  of  the 
impact  of  decentralization  is  revealed  by  a  study  of  its  implications  upon 
performance  for  a  specific  system.  Performance  of  the  two-mass  oscillator 
system  incorporating  several  decentralized  Kosut  designs  is  examined. 

Specific  limitations  on  system  performance  are  observed  when  compared  with 
corresponding  centralized  designs. 

The  section  is  concluded  with  a  brief  statement  of  important  questions 
that  remain  to  be  resolved  regarding  the  extended  Kosut  approach. 

2.2  Notation 


In  order  to  be  self-contained,  the  notation  to  be  used  is  briefly 
summarized.  It  is  assumed  that  the  structural  dynamics  of  systems  to  be 
discussed  can  be  adequately  represented  by  a  finite-element  model: 

Mq  +  Kq  *  f  (2-1) 


T 

where  q  =  (q^,...,qn)  is  the  vector  of  generalized  physical  coordinates, 

M:nxn  is  the  (positive  definite)  mass  matrix,  K:nxn  is  the  (positive 

T 

semidef inite)  stiffness  matrix,  and  f  =  (f f  )  is  the  vector  of  forces 

1  n 

applied  to  the  structure.  Equation  (2-1)  may  be  transformed  to  an  equivalent 
representation  in  modal  coordinates: 


U  + 


2 

irn 


*Tf 


T 

<J>  B. u 
A 


(2-2) 
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where  Q  : n*n  is  the  diagonal  matrix  containing  squares  of  the  natural 
2  2 

frequencies  u)^, . .  •  ,0)^,  4>:nxn  is  the  nonsingular  matrix  whose  columns  are 
the  structural  mode  shapes  <J>\ . . .  normalized  such  that 


$TM4>  =»  I  ,  *TK*  =  ft2 

n 


(2-3) 


B,:nxm  is  the  actuator  influence  matrix  for  the  control  u  =  (u,,...,u  )  , 
A  T  1  m 

and  n  =  (ri^,...»nn)  is  the  vector  of  modal  coordinates  defined  by 


q  = 


(2-4) 


Superscript  "T"  denotes  matrix  transpose,  and  the  notation  I  denotes  the 
nxn  identity  matrix.  The  linear  measurement  equation  has  the  form: 


y  =  cpq  +  V 


(2-5) 


where  y  =  (y^,...,y£)  is  the  vector  of  system  outputs,  and  C^rJlxn,  Cy:£xn 
are  the  displacement  and  rate  measurement  coefficient  matrices,  respectively. 

The  structural  model  (2-2)  may  be  placed  in  a  convenient  state  space 
form  as  follows.  Without  loss  of  generality,  it  can  be  assumed  that  the 
components  of  n  have  been  reordered  so  that  the  first  N  (N  <  n)  of  them  are 

T 

the  selected  critical  modes  nc  =  (r> » -  -  ♦  »r»N)  »  and  that  the  remainder  are 
residual  modes  nD  =  (nM, i > • • • >1  )T*  Corresponding  submatrices  ft2:NXN, 
(n-N)x(n-N)  of  £r  and  <I>c:nXN,  $R:nx(n-N)  of  $  are  defined  similarly. 

With  x  =  (n  ,  nr,  nD)  as  a  state  vector,  the  structural  model  (2-2), 

v*  K  K 

(2-4),  (2-5)  becomes: 


Ax  +  Bu 


(2-6) 


y  =  Cx 


(2-7) 
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where 


Using  residual  mode  truncation  as  the  model  reduction  process,  a  reduced 


order  model  with  the  state  xc  .  (nc»  0 c>  is  obtained 


A„x„  +  B-u 


(2-8) 


2. 3  New  Results  for  Centralized  Controllers 

2.3.1  Stability  and  Robustness  for  Output  Feedback  Designs 

2. 3. 1.1  General  Criteria  for  Stability 

As  noted  previously,  one  of  the  principal  objections  to  the  use  of 
the  extended  Kosut  method  for  controller  design  is  the  absence  of  any 
useful  conditions  which  guarantee  closed-loop  stability  even  at  the  design 
model  level.  The  approach  taken  to  developing  such  conditions  has  been  to 
try  to  discover  an  appropriate  subclass  of  the  systems  (2-8),  (2-9)  being 
considered  for  which  the  extended  Kosut  approach  leads  to  a  stable  closed- 
loop  system  at  the  design  model  level. 

The  first  result  characterizes  the  eigenvalues  of  the  closed-loop 
design  model  for  general  output  feedback  systems  in  terms  of  fundamental 
model  parameters. 

Theorem  2-1.  Assume  that  the  reduced-order  model  (2-8) ,  (2-9)  is  connected 
by  an  output  feedback  law:  u  «  Gy^,.  Then  the  eigenvalues  of  the  closed- 

loop  system  are  the  roots  of  the  (2N)th  degree  polynomial  equation: 

det  j  *CT  [-A2M  +  *BAGCy  -  (K  -  V^p)]  *c  [  “  0  D  (2-10) 

The  proof  is  given  in  Section  2.5.  A  remark  about  Equation  (2-10)  is 
appropriate.  Since  the  mode  shapes  $  ,...,<pn  are  linearly  independent, 
the  matrix  <t>c  =  [(J)1  !  •••  |  $N]  :n*N  has  maximum  rank.  It  does  not  follow  that 
equation  (2-10)  can  be  replaced  with  the  (2n)th  degree  polynomial  equation: 

det  [-A2M  +  ABAGCy  -  (K  -  BAGCp)  ]  -  0  (2-11) 

The  latter  equation  is  easily  recognized  as  the  closed  loop  characteristic 
equation  for  the  full-order  structural  model  (2-1),  (2-5)  connected  with  the 
output  feedback  law  u  *  Gy.  (Set  A  -  ito  in  equation  (2-11)  to  obtain  an 
equivalent,  more  familiar,  form.)  Equation  (2-11)  has  2(n  -  N)  more 
solutions  than  equation  (2-10). 

A  simple,  but  useful,  sufficient  condition  for  closed-loop  design 
model  stability  is  based  on  Theorem  2-1.  For  clarity,  the  result  is 
presented  for  the  case  where  rate  sensors  only  are  used  (Cp  *  0). 

Corresponding  results  for  the  general  case  are  easily  obtained. 


( 

i 


0. 


Theorem  2-2.  Assume  a  closed-loop  structure  as  In  Theorem  2-1  with  C  * 

t  * 

Suppose  that  there  exist  6^  >  0,  i  -  such  that  ^B^GCyC^  *  diag  (-6^) . 

Then  the  closed-loop  system  is  stable.  Q 

The  proof  is  given  in  Section  2.5.  This  result  is  useful  for  establishing 
robustness  of  stability. 

Specializing  from  general  output  feedback  systems  to  those  designed 
using  the  extended  Kosut  method,  one  may  ask  under  what  conditions  the 
feedback  gain  matrix  generated  by  that  method  satisfies  the  hypotheses  of 
Theorem  2-2.  This  question  has  not  been  answered  in  general.  It  should 
be  evident,  however,  that  Theorem  2-2  is  by  no  means  necessary  for  design 
model  stability;  it  is  in  fact  a  quite  restrictive  condition.  This  may  be 
seen  by  examining  the  condition  in  the  context  of  the  two-mass  oscillator 

T  2  T  2 

example  worked  out  previously  [5;Sec.  6].  In  that  case,  <J>C  B^GCy^  ■  ($  )  , 

a  scalar  (hence  diagonal),  where 


G  *  G(a;e,6) 


♦2  2 

-a  +  e 


-5 


-e 


-a  + 


♦  l  2 

6  V*2 


The  stability  condition  of  Theorem  2-2  reduces  here  to  the  condition  that 
G  be  negative  definite,  which  requires  [TO]  that: 


and 


y2  2 

o  -  e  /<}>.  >  0 


[- 


2  2 

*22  *12 
e  -  6 


However,  it  was  seen  [5;Sec.  6]  that  this  closed-loop  design  model  Ac  +  BCGCC 

is  independent  of  the  free  parameters  e  and  6,  and  is  in  fact  stable  as 
long  as  a  >  0. 


II 


2.3. 1.2  Robustness  of  Stability  for  the  Extended  Kosut  Method 

The  robustness  question  is  examined  in  two  parts.  First,  the 
sensitivity  of  the  design  model  closed-loop  eigenvalues  is  examined  relative 
to  changes  in  the  feedback  gain  matrix.  This  is  applicable  to  any  output 
feedback  system.  Second,  the  sensitivity  of  the  gains  generated  by  the 
extended  Kosut  method  is  examined  relative  to  changes  in  the  fundamental 
structural  model  parameters.  For  clarity,  it  is  again  assumed  that  rate 
sensing  only  is  used  (Cp  =  0). 

Theorem  2-3.  Assume  a  closed-loop  structure  as  in  Theorem  2-1  with  Cp  =  0. 

The  mapping  which  assigns  to  each  output  feedback  gain  matrix  the  set  of 
closed-loop  eigenvalues  is  continuous.  D 

The  proof  is  given  in  Section  2.5.  The  significance  of  Theorem  2-3  is 
as  follows.  Suppose  that  for  a  particular  gain  Gq  the  system  is  stable. 

Then  the  corresponding  eigenvalues  A^  E  (A*j*, . . . ,  A^j)  lie  in  the  Cartesian 

product  of  the  sets  ReCA^)  <  0,  an  open  set  in  real  Euclidean  2N-space. 

Theorem  2-3  says  that  there  is  an  open  set  containing  Gq  in  the  space  of 

mx£  real  matrices  for  which  the  corresponding  eigenvalues  also  satisfy 
Re(A^)  <  0,  i  =  1,...,2N.  It  should  be  noted  that  Theorem  2-3  is 

essentially  without  significance  if  at  least  one  stable  set  of  eigenvalues 
cannot  be  found  (cf.  Theorem  2-2). 

For  the  following  result,  it  is  assumed  that  structural  model  reduction 
is  by  truncation  (equations  (2-8),  (Z-9)),  and  that  a  specific  reference 
system  and  a  specific  suboptimal  performance  measure  for  the  extended  Kosut 
method  (minimum  error  excitation)  have  been  chosen.  The  mapping  described  is 
then  well  defined. 

Theorem  2-4.  Assume  Cp  =  0.  The  mapping  which  assigns  to  each  mass¬ 
stiffness  matrix  pair  in  Equation  (2-1)  the  extended-Kosut  feedback  gain 
matrix  (Moore-Penrose  inverse  particular  solution  [8])  is  continuous.  0 

The  proof  is  given  in  Section  2.5.  Together  with  Theorem  2-3,  this  result 
shows  that  a  stable  design  can  be  achieved  by  the  extended  Kosut  method, 
then  stability  is  retained  under  appropriately  small  changes  in  the 
fundamental  structural  model  parameters. 

2.3.2  Refinement  of  Previously  Reported  Results 

Extensions  to  the  Kosut  method  reported  previously  [5]  are  based  on 
the  existence  and  properties  of  solutions  G:mx£  to  the  Kosut  gain  equation: 

g(ccPCcT)  -  F*PCcT  (2-12) 
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where  F*:m*2N  is  the  state  feedback  gain  for  the  stable  reference  system  used 
in  the  suboptimal  design,  and  P:2Nx2N  is  the  positive  definite  Lagrange 
multiplier  matrix  obtained  from  the  Liapunov  equation: 

(Ac  +  BCF*)P  +  P(AC  +  BcF*)T  +  I2n  =  0  (2-13) 

It  was  shown  [5;para.  6. 2. 3. 3]  that  Equation  (2-12)  is  solvable  regardless 

T 

of  the  rank  of  C  (equivalently,  rank  of  C  PC  ),  and  that  in  the  rank- 

L  v  L 

deficient  case  (say,  rank  Cc  =  r  <  £  £  2N) ,  the  general  solution  has  the  form: 

G(D  =  Gq  +  TS  (2-14) 

where  Gq  is  a  particular  solution  of  Equation  (2-12),  T  is  a  matrix  whose 

m*(i.-r)  non-zero  elements  are  free  parameters,  and  S  is  an  appropriate  non¬ 
singular  matrix. 

In  a  detailed  design  example,  it  was  observed  that  the  eigenvalues  of 
the  closed-loop  design  model  were  independent  of  the  (two)  free  parameters 
in  the  feedback  gain  matrix.  Since  the  performance  (as  contrasted  with  the 
stability)  of  the  design  model  is  not  of  practical  interest,  it  was  possible 
to  base  the  selection  of  free  parameter  values  solely  on  full-order  system 
performance  considerations.  The  question  arises  as  to  what  conditions 
characterize  this  "independence"  in  general.  To  begin,  a  very  simple 
sufficient  condition  can  be  stated. 

Fact  2-5.  Assume  a  closed-loop  structure  as  in  Theorem  2-1;  refer  to 
Equation  (2-14).  If  B^TSCq  =  0,  then  the  eigenvalues  of  the  closed-loop 

design  model  are  independent  of  the  free  parameter  matrix  T.  0 

The  proof  is  self-evident,  since  both  G(T)  and  Gq  generate  the  same  closed- 

loop  system  matrix.  It  is  worth  observing  that  this  condition  (in  fact, 
even  the  stronger  condition:  TSCq  -  0)  has  been  satisfied  in  all  numerical 

examples  studied  to  date.  The  following  theorem  fully  characterizes  the 
independence  condition. 

Theorem  2-6.  Assume  a  closed-loop  structure  as  in  Theorem  2-1;  refer  to 
Equation  (2-14).  The  eigenvalues  of  the  closed-loop  design  model  are 
independent  of  the  free  parameter  matrix  T  if  and  only  if  A  +  B  G(T)C  and 

Aq  +  BqGqCq  are  similar;  i.e.,  if  and  only  if  there  exists  a  nonsingular 

matrix  Q:2N*2N  such  that: 


[Ac  +  BcG(r)cclQ  -  Q[Ac  +  Wc1  *  0  •  D 
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The  proof  is  given  in  Section  2.5.  The  full  implications  of  this  freedom 
of  choice  in  selecting  a  feedback  gain  matrix  with  the  extended  Kosut 
method  have  not  yet  been  uncovered.  Performance  improvements  for  low-order 
vibrating  systems  have  already  been  demonstrated  [5,6].  It  is  worth 
observing,  however,  that  the  question  of  effective  use  of  available  freedom 
of  choice  in  selections  of  feedback  gains  is  an  important,  and  largely 
unexplored,  question  for  general  multivariable  feedback  systems  [13;p.51]. 

Proof  of  the  general  algebraic  consistency  of  Equation  (2-12)  has 
been  considerably  simplified  (cp.  [6;Theorem  2],  [5;Theorem  6-6]).  Since 
this  is  a  principal  result,  the  streamlined  proof  is  of  interest.  Details 
are  given  in  Section  2.5. 

2.4  New  Results  for  Decentralized  Controllers 
2.4.1  New  Perceptions  Suggested  by  Decentralization 

The  discussion  of  decentralization  in  this  section  of  the  report  is 
restricted  to  the  nature  of  the  information  structure  in  the  sensor-to- 
actuator  feedback  connection.  Totally  decentralized  systems  which  include 
decentralized  plant  structures  are  not  treated.  Specifically,  it  is  assumed 
that  the  information  for  driving  a  given  actuator  is  derived  from  a  (possibly 
proper)  subset  of  all  available  sensors,  and  that  the  collection  of  sensors  in 
this  subset  is  (in  general)  different  from  one  actuator  to  another.  In 
symbols,  the  sensor  combinations  for  transmission  to  each  actuator  are 
generated  from  the  state  vector  as  follows: 

y1  =  C^x  ,  C1:£ix4N  5  i  =  1 . m  (2-15) 


Equation  (2-15)  replaces  Equation  (2-7)  in  a  decentralized  full-order  model. 
The  feedback  connection  is  made  by: 


T  i 

-  8ty 


Wl 


i  ~  1  j  •  •  •  j  in 


(2-16) 


Remark:  In  the  special  case  that  the  components  of  each  y^  are  selections  from 

_  T 

the  components  of  a  fixed  vector  y  =(y  ,...,y^)  ,  Equation  (2-16)  may  be 

_  T 

written  in  the  form  u  =  Gy,  where  u  =  (u,,...,um)  ,  and  G  =  ]  has  the 


r 


m 


property  that  g  =  0  if  the  component  y  of  the  fixed  vector  y  is  not  also 

i  J 

a  component  of  the  output  vector  y  associated  with  u^.  In  general,  however, 
this  situation  will  not  necessarily  exist. 
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The  Kosut  gain  equations  corresponding  to  Equation  (2-12)  for  a 
controller  with  decentralized  feedback  structure  are  [5] : 


^(c*PC*  j  =  F*PC^  ,  i  =  1, . . .  ,t 


(2-17) 


*  *  i 

where  F  :lx2N  are  the  rows  of  F  ,  and  C  :Z  x2N  are  the  matrices  corresponding 

1  L  1 

to  Equation  (2-15)  for  a  reduced-order  model.  The  structure  of  solutions  to 
algebraically  consistent  equations  of  the  form  (2-17)  may  be  expressed  in 
the  following  form. 

Theorem  2-7.  Let  A:A*V,  b:l*v  be  matrices  for  which  the  equation  xA  =  b 
has  a  solution  x^.  Denote  r  =  rank  (A).  The  general  solution  of  this 

equation  has  the  form: 


x 


*o  + 


(2-18) 


where  the  A-r  row  vectors  of  S:(A-r)*A  form  a  basis  for  the  null  space  of 
A  (viewed  as  a  mapping  of  E^  into  EV) ,  and  the  (A-r)-vector  y  is  arbitrary.  Q 

The  proof  is  given  in  Section  2.5.  This  result  immediately  suggests  an 
appropriate  reformulation  of  previously  reported  solution  structures  for 
homogeneous  matrix  equations  [cf.  5,  Theorem  6-3;  6,  Theorem  3;  8,  Theorem  3-9] 
which  sharpens  the  statement  and  identifies  the  true  nature  of  the  "S" 
matrix. 

Theorem  2-8.  Let  A:A*v>  be  a  matrix.  Denote  r  =  rank  (A).  The  general 
solution  of  the  equation  XA  =  0^  has  the  form  X  =  TS,  where  the  A-r  row 

vectors  of  S:(A-r)*A  form  a  basis  for  the  null  space  of  A  (viewed  as  a 

mapping  from  E^  into  EV) ,  and  the  matrix  r:yx(A-r)  is  arbitrary.  D 

This  result  shows  that  the  last  A-r  rows  of  the  S  matrix  in  Theorem  6-3  of 

T 

[5]  and  the  last  H-r  rows  of  the  U  matrix  in  Theorem  3-9  of  [8]  correspond 
simply  to  different  choices  of  a  basis  for  the  null  space  of  A. 

Step  4  of  the  algorithm  for  numerical  solution  of  the  Kosut  gain 
equation  [8,  Sec.  3.5]  should  now  be  modified  so  as  to  be  consistent  with 
Theorem  2-8. 

Step  4  (revised) .  Find  the  general  solution. 

(a)  Select  an  m*(£-r)  matrix  f  of  free  parameters. 
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(b)  Choose  row  ^-vectors  that  form  a  basis  for  the  null 

space  of  A  and  form  the  matrix 


S  & 


£-r 


(One  possible  choice  is  the  set  of  vectors  obtained  by  transposing 
columns  r+l,...,£  of  the  matrix  U  from  Step  1(a)). 


(c)  Compute  the  homogeneous  solution  X  =  TS. 

c 

(d)  Compute  the  general  solution  X  =  X  +  X  ,  using  the  results  of 

Steps  3(b)  and  4(c).  Q  p 


Step  4(b)  makes  explicit  the  fact  that  there  is  considerable  freedom  in 
choosing  the  matrix  S  (although  not  complete  arbitrariness).  Specific 
approaches  to  exploiting  this  freedom  of  choice  await  further  investigation. 


2.4.2  Performance  Implications  of  Decentralization 

There  are  a  number  of  reasons  why  a  decentralized  information  structure 
might  be  incorporated  in  a  physical  system.  For  a  very  large  structure  with 
many  sensors,  the  computational  burden  of  incorporating  information  from 
all  of  the  sensors  in  the  determination  of  control  signals  to  each  actuator 
may  far  outweigh  the  advantages  when  a  centralized  processing  unit  is  used. 

If  numerous  processing  units  are  distributed  throughout  the  structure,  it 
may  not  be  feasible  to  transmit  information  from  all  the  sensors  to  each 
processing  unit.  On  the  other  hand,  there  are  cases  in  which  the  only 
relevant  information  is  that  from  nearby  sensors  (e.g.,  in  the  case  of 
shape  control) ,  in  which  case  information  from  all  the  sensors  at  each 
processing  unit  would  be  largely  unnecessary.  Whatever  the  reasons,  however, 
it  is  important  to  take  these  information  structure  constraints  explicitly 
into  account  in  the  mathematical  model  used  for  controller  design  and  eval¬ 
uation.  In  general,  one  would  expect  some  loss  in  performance  (by  whatever 
measure)  as  the  price  of  operating  with  restrictions  on  the  internal 
information  structure  for  determining  control  signals. 


A  clear  picture  of  such  performance  implications  does  not  emerge  from 
the  design  equations  themselves.  In  fact,  the  gain  equations  (2-17)  for  the 
decentralized  controller  have  exactly  the  same  mathematical  structure  as  the 
gain  equation  (2-12)  for  the  centralized  controller.  Hence,  all  the  results 
obtained  for  the  latter  (e.g.,  algebraic  consistency,  structure  of  general 
solution)  apply  to  the  former.  Therefore,  some  design  calculations  were 
made  on  a  specific  representative  system  (the  two-mass  oscillator  example  [5]) 
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to  get  at  least  a  snapshot  of  the  nature  of  performance  degradation,  if  any, 
associated  with  decentralization  of  the  controller.  The  oscillator  was 
assumed  to  have  one  force  actuator  and  one  rate  sensor  associated  with  each 
mass  [5;  Fig.  2-2].  Each  of  three  possible  types  of  decentralization  in  the 
feedback  structure  was  studied:  (a)  fully  decentralized:  =  u^(y^), 

u2  =  u2^y2^’  (b)  partially  decentralized  with  actuator  1  localized: 

U1  =  ui^yi^’  u2  =  u2^yl’  y2^:  ^ c ^  Partially  decentralized  with  actuator  2 

localized:  u^  =  Uj(y^,  y 2),  u2  =  u2<y2>.  The  reference  system  for  the 

Kosut  design  was  taken  as  determined  in  previous  examples  [5].  Hence  F*, 
and  also  P,  [cf.  Equation  (2-13)]  are  as  previously  computed: 

0.066389774  0.51609641 


-0.028241581  -0.21954252 


P 


1.1373892  -0.5 

-0.5  7.4375655 


2.4. 2.1  Fully  Decentralized  Desis 


Recall  that  the  state  vector  for  the  reduced-order  model 
[cf.  Equation  (2-8)]  is  xc  =  (f]2>  t^)*  that  residual  mode  (r^)  dynamics  are 

being  ignored,  and  that  Equation  (2-4)  relates  the  physical  and  modal 

coordinates.  The  specifications  for  this  configuration  may  be  stated: 

1*2*  i 

y  =  q^,  y  =  q2-  Accordingly,  the  matrices  Cc  corresponding  to 

Equation  (2-15)  for  the  reduced-order  system  are: 


CJ  •  [»  »?]  •  CC  ■  [°  *22]  • 

Since  each  C*  has  maximum  rank,  each  of  the  gain  equations  (2-17)  has  a 
unique  solution.  The  solutions  are: 


where  =  0.59708155  is  the  reference-system-dependent  parameter  occurring 
in  the  solution  of  Equation  (2-12)  for  the  centralized  design  [5,  Eq.(6-43)]. 
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The  closed  loop  design  model  has  the  same  characteristic  polynomial  as  in 
the  centralized  case: 

X2  +  0°(<J)2  <p2)\  +  a)2 

It  is  stable,  with  u)  =  (open-loop  critical  mode  frequency),  and  £  =  0. 1 

(design  objective).  Note  that,  in  contrast  to  the  corresponding  centralized 
design,  there  are  no  free  parameters  in  the  solution  for  the  feedback  gains. 

The  matrices  C1  of  Equation  (2-15)  for  the  full-order  system: 

c1  £  [o  <j>2  o  ♦}]  ,  c2  =  [o  4>2  0  4>l] 

are  used  to  connect  the  reduced-order  controller  to  the  full-order  system 
via  Equation  (2-16).  The  characteristic  polynomial  for  the  closed-loop 
full-order  system  is 


Dynamic  characteristics  of  this  system  are: 

Mode  1  (critical):  =  2.5851902  sec 

=  0.099939961 

Mode  1  (residual):  =  0.54704429  sec 

he 
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0.34631083 


Some  effects  of  decentralization  on  overall  system  performance  can  be  seen 
by  comparing  the  time  response  in  Figure  2-1  and  the  frequency  response  in 
Figure  2-2  with  the  corresponding  plots  for  the  centralized  design 
[5;  Figs.  6-7  and  6-8,  respectively].  Noticeable  degradation  does  occur,  and 
is  traceable  to  the  reduction  in  damping  of  the  residual  mode.  Nonetheless, 
the  degradation  is  not  nearly  as  severe  as  that  occurring  when  the  sensor 
on  the  outer  mass  of  the  two-mass  system  is  deactivated  while  retaining  a 
centralized  information  structure  on  the  remaining  available  sensor 
[5;  Figs.  6-13  and  6-14,  respectively]. 


2. 4. 2. 2  Partially  Decentralized  Design:  Actuator  1  Localized 

1  • 

The  specif ications  for  this  configuration  may  be  stated:  y  =  q, , 

2  •  •  T  i  1 

y  =  (<1|»  q2)  •  The  related  matrices  C^,  corresponding  to  Equation  (2-15) 

for  the  reduced-order  system  are: 

°c  -  [°  ■  °c  •  t°  *2] 


Solutions  to  the  gain  equations  (2-17)  are  as  follows: 


82(«) 


-6  -a0  +  S^1/*2 


] 


where  5  is  a  free  parameter.  An  equivalent  feedback  connection  of  form 
u  =  Gy^,  exists  and  is: 


T 

r  o  .  i 

81  ° 

i 

Q 

a 

G(6)  = 

= 

2 

g2(<S) 

-6  -o°  +  6  1M?\ 

The  characteristic  polynomial  for  the  closed-loop  design  model  is  the  same 
as  that  for  the  preceding  design;  in  particular,  it  is  independent  of  the 
free  parameter  6. 

The  matrices  C1  of  Equation  (2-15)  for  the  full-order  system  are: 

C1  =  [o  <f>J  0  ,  C2  =  [o  <p2  0  (f)1] 
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SYSTEM:  Q101  I.C.  Q2  (0) 


response  to  initial  conditions 
■  decentralized  controller). 


To  assist  in  determining  an  appropriate  value  for  the  free  parameter, 
the  modal  closed-loop  equations  are  developed: 


(2-20) 


Several  possibilities  for  selecting  6  are  as  follows:  (a)  choose  6  so  that 
the  coefficient  of  in  Equation  (2-19)  vanishes;  (b)  choose  6  so  that  the 

coefficient  of  in  Equation  (2-20)  attains  some  desired  value;  or 

(c)  some  combination  of  (a)  and  (b).  Selection  (a)  eliminates  the  effect  on 
the  critical  mode  of  residual  mode  dynamics,  whereas  selection  (b)  approx¬ 
imates  a  specification  for  the  residual  mode  damping.  Past  design  examples 
have  indicated  that  failure  to  damp  the  residual  mode  strongly  and  adversely 
affects  the  full  system  performance,  and  that  the  spillover  effect  admitted 
by  Equation  (2-19)  is  a  lesser  concern.  Since  selection  (a)  will  not  affect 
residual  mode  damping,  whereas  selection  (b)  will  tend  to  reduce  the 
spillover  problem  indicated  by  Equation  (2-19),  selection  (b)  is  made.  Thus: 


6°  & 


2CiD“i  -  °°(*1  4*1) 

(*>2) 


=  0.33580542 


det  $ 


where  is  chosen  to  be  1  /-/l.  The  gain  values  then  become: 


G(6U)  = 


r-o. 
L  -0. 


59708155 

33580542 


-1.3864866 


] 


With  this  connection,  the  dynamic  characteristics  of  the  system  are: 


Mode  2 

(critical) ; 

“2C 

2.5914658  sec' 

C2C 

0.10021320 

Mode  1 

(residual) : 

W1C 

0.54571957  sec 

?1C 

0.70633967 
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Performance  of  this  design  is  so  close  to  that  of  the  corresponding 
centralized  design  that  no  degradation  can  be  observed  from  the  plots. 
Figures  2-3  and  2-4  may  be  compared  with  the  corresponding  plots  for  the 
centralized  system  [5;  Figs.  6-7  and  6-8,  respectively].  Consequently, 
Figures  2-3  and  2-4  may  themselves  be  viewed  as  representing  the 
centralized  system  time  and  frequency  response,  respectively,  for  the 
purposes  of  the  comparisons  in  this  section. 


2. 4. 2. 3  Partially  Decentralized  Desif 


Actuator  2  Localized 


The  specifications  for  this  configuration  may  be  stated: 

1  *  •  T  2  •  i 

y  =  (q^i  q2)  »  y  =  q2-  The  related  matrices  C^,  corresponding  to 

Equation  (2-15)  for  the  reduced-order  system  are: 


cc  ■  [°  *2]'  cc  ’  [°  *2] 


Solutions  to  the  gain  equation  (2-17)  are  as  follows: 


-  [-o°  +  c^  -e]  , 


80  =  -o 


where  e  is  a  free  parameter.  An  equivalent  feedback  connection  of  form 
u  =  Gyc  is: 


G(e)  = 


T 

gj(e) 


0  8 


0  ^2  2 

-<r  +  e  -e 


The  characteristic  polynomial  for  the  closed-loop  design  model  is  the  same 
as  for  the  fully  decentralized  design;  it  is  also  independent  of  the  free 
parameter  e. 

The  matrices  C*  of  Equation  (2-15)  for  the  full-order  system  are: 

c1  =  [0  4>2  0  (j)1]  ,  c2  =  [0  <j>2  0 
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The  modal  closed-loop  equations  (to  help  select  e)  are: 


o2  +  o°(<j>2  4>2)n2  +  ^2  *  M*2  +  e  det  *]  nL 

>1  +  4*1)  -  e  det  4>  (  1/<l>2)]n1  +  WjOj  =  [-o°(^2  (j)1)]  n2 


(2-21) 


Following  the  same  arguments  as  in  the  preceding  section,  the  parameter  e 
is  chosen  so  that  the  coefficient  of  n,  in  Equation  (2-21)  has  the  desired 
value  =  1//5: 


0 

e 


A 


2C1du>l  -  a0^1  t1) 
det  *  (4>}/-4>2) 


0.92785918 


The  gain  values  then  become: 


G(e°) 


-0.99178405 

0 


-0.92785918 

-0.59708155 


With  this  connection,  the  dynamic  characteristics  of  the  system  are: 


Mode  2 

(critical) : 

U)-  _  = 

2C 

2.5683141  sec" 

C2C 

0.098354234 

Mode  1 

(residual) : 

U1C  = 

0.55063889  sec 

C1C 

0.71291349 

The  only  performance  degradation  relative  to  the  centralized  design 
observable  from  the  plots  is  in  the  frequency  domain,  and  is  quite  minor. 
The  frequency  response  shows  slightly  higher  magnitudes  near  u)  * 

(open- loop  critical  mode  frequency),  and  slightly  higher  phase  shifts  for 
w  2  Figures  2-5  and  2-6  may  be  compared  with  the  corresponding  plots 

for  the  centralized  system  [5;  Figs.  6-7  and  6-8,  respectively]. 
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Figure  2-6.  Frequency  response  of  qj  to  periodic  disturbance  fj(t)  “  sin  ut 
(partially  decentralized  controller:  actuator  2  localized) . 


2. 4. 2. 4  Remarks  on  Performance  Comparisons 

A  summary  of  the  design  studies  on  the  two-mass  oscillator  example 
for  all  of  the  centralized  and  decentralized  designs  discussed  in  this 
section  is  shown  in  Table  2-1.  Data  for  the  centralized  designs  was 
previously  reported  [5,7]. 

Broad  conclusions  regarding  the  effects  on  performance  of  decen¬ 
tralizing  the  feedback  information  structure  cannot  be  drawn  from  these 
few  designs  on  such  a  low-order  system.  However,  several  observations  are 
worth  noting.  First,  noticeable  performance  degradation  did  occur,  and  the 
degree  of  degradation  increased  as  the  degree  of  decentralization  increased. 
Second,  the  worst  performance  degradation  (for  the  fully  decentralized 
conf iguration)  was  not  as  bad  as  one  might  have  expected.  Finally,  the 
design  flexibility  due  to  the  free  parameters  occurring  in  the  designs  with 
partial,  but  less  than  total,  decentralization  was  sufficient  to  eliminate 
(nearly)  all  performance  degradation. 

2.5  Proofs  of  Theoretical  Results 


The  proof  of  Theorem  2-1  depends  on  the  following  two  facts. 

USS-l-  *  "c 

Proof.  It  follows  from  a  block  matrix  representation  of  Equation  (2-3) 
that  (^)V  =  and  (<5>i)TK<t>'1  *  wi^ij  ’  ^*3  *  l»...,n,  where  6^  is  the 

Kronecker  indicator  function: 

1  , 

0  , 

In  particular,  these  relations  hold  for 

Fact  2.  (Rosenbrock  [9]).  If  A,  B,  C, 
square  and  A  nonsingular,  then: 

TA  ■  B1  -i 

det  -  -I - -  det  A  det  (-CA  B  +  D) 

LC  i  Dj 


i  =  j 
i  *  3 

i,  j  =  1 . N.  I 

D  are  matrices,  with  A,  D  being 
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Proof  of  Theorem  2-1.  The  closed-loop  system  matrix  for  the  model  (2-8) , 
(2-9)  under  the  connection  u  =  Gy^  is: 


+  “c^c 


-ac +  '!wc 


It  is  easily  seen  from  the  properties  of  determinants  that  the  vanishing  of 

2  T 

det  (A^  +  B^GC^,)  is  equivalent  to  the  condition:  det  [-ft^,  +  ^B^Cp^]  = 

This  proves  result  (2-10)  for  the  case  X  =  0.  For  X  ^  0,  the  matrix  -Xl^ 
is  nonsingular.  Using  Fact  2  and  the  properties  of  determinants  leads 
to: 


det  (Ac  +  BCGCC 


*I2N)  =  (-D 


det 


-x\  +  4c[bac«cv  +  °Pmc 


Result  (2-10)  follows  by  an  application  of  Fact  1.  R 

Proof  of  Theorem  2-2.  Denote  A  =  diag  (-6^).  Using  Theorem  2-1  and  Fact  1, 
the  eigenvalues  of  the  closed-loop  system  are  **een  to  be  the  roots  of: 


0  -  <et(-AV>‘-ac)  -  <-»*  n>2  +  V-i) 

Since  each  6^  >  0,  all  eigenvalues  have  negative  real  parts.  H 


Proof  of  Theorem  2-3.  Denote  by  G  t- 


->  (Xj , . . .  »Y2jj)  the  mapping 


described  in  the  theorem.  The  function  y  is  a  composition  Y^  °  Y2  °  Y^ 

Y  i  ^ 

of  the  following  three  mappings:  (1)  G  I — — ►  $^,B^GC^$^; 

(2)  Z  I — Pz(’)»  where  PZ(X)  =  det  ^X2I^  +  XZ  +  fi2  )  ;  and 


(3)  p(‘)  h 


->  (X . . . .X^^) ,  the  roots  of  p(X)  =  0.  It  can  be  shown  that 


Yp  Y2>  and  Y^  are  continuous;  hence  y  is  continuous.  Continuity  of  y^  is 

based  on  the  well-known  root-coefficient  relations  for  polynomial 
equations  [11] .  I 
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Proof  of  Theorem  2-4.  Denote  by  (M,K)  | - >  G  the  mapping  described  in 

the  theorem.  The  function  8  is  a  composition  f  °  6  °  8,  of  the  following 

ft  J  ^  1 

P1  2  2  /  2  2\ 
three  mappings:  (1)  (M,K)  | — — >  (ft  ,  $) ,  where  ft  =  diag  (id  ,  I, 

2  2  1  ' 
uk  being  the  roots  of  det  (-id  M  +  K)  =0,  and  =  [<#»  j  j  tp ] ,  <[>  being 

the  solutions  of  ^-id^M  +  <f>  =  0  such  that  Equation  (2-3)  holds; 

2  ^2 

(2)  (ft  ,  4>)  | — — >  (A  ,  B  ,  C  ),  the  latter  matrices  being  defined  by 

^  ^  ^  ^ 

Equations  (2-8),  (2-9);  and  (3)  (Ac,  Bc>  Cc)  (— — >  G  =  (F*PCCT) (CCPCCT)+ , 

where  F  is  the  state  feedback  gain  for  the  Kosut  reference  system,  and  P 
is  the  solution  of  Equation  (2-13).  The  superscript  "+"  denotes  the 
Moore-Penrose  inverse  [8].  It  can  be  shown  that  8^»  82*  and  8^  are 
continuous;  hence  8  is  continuous.  I 


result . 


The  proof  of  Theorem  2-6  is  self-evident  from  the  following  general 


Fact  3.  Two  square  matrices  are  similar  if  and  only  if  they  have  the  same 
set  of  eigenvalues. 

Proof.  (Only  if)  This  is  a  well-known  result  [12;  p.144]. 

(If)  Since  both  matrices  have  the  same  set  of  eigenvalues,  they  are  each 
similar  to  the  same  Jordan  block  matrix.  Since  similarity  is  an  equi¬ 
valence  relation,  the  matrices  in  question  are  similar  to  each  other. 

Simplified  proof  of  Theorem  6-6  of  [5;  Sec.  6], 

VI 

Denote  E  =  | - .  Since  the  rank  of  a  matrix  is  equal  to  the  number  of 

,* 

T  T 

its  linearly  independent  rows:  rank  (C  PC  )  5  rank  (EC  ).  Conversely, 

”  ^  ^  ip 

suppose  that  this  inequality  is  strict.  Denote  r  =  rank  (EC„  ).  Using 

T  ^ 

the  fact  that  rank  (C^PCC  )  =  rank  (C^)  when  P  is  positive  definite,  it 

follows  that  rank  (Cr)  <  r.  But  then  it  follows  easily  [12]  that 

*T"i  rri 

rank  (EC  )  -  min  (rank  (E) ,  rank  (CpA)  }  <  r,  which  contradicts  the 
definition  of  r.  Hence,  rank  (C^PC^  )  -  rank  (EC^,  ).  I 

The  proof  of  Theorem  2-7  is  based  on  the  following  well-known  result 
from  linear  algebra  [12], 
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Fact  4.  Let  <(■/  be  a  linear  transformation  defined  on  a  vector  space  X. 
Then: 


dim  [Range  GV)  ]  +  dim  [Null  CV)  ]  =  dim  (X) 


Proof  of  Theorem  2-7. 


Define  the  map  x  xA  on  the  row-space  e\ 


Since  rank  (tr/)  =  rank  (A)  =  r,  it  follows  from  Fact  4  that  Null  (.V),  the 
null  space  of  <V,  has  dimension  A-r.  If  r  =  A,  then  Null  G*/ )  =  0  and  x^ 
is  unique.  Otherwise,  r  <  A,  so  there  exist  row  vectors  s^,...,s^_r 


Let 


A  T  T  1  |  t  ~l  T 

that  form  a  basis  for  Null  (,«/).  Denote  S  =  s,  i  •••  i  s,  . 

x  L  1  •  i  X-rJ 

x  e  E  satisfy  xA  =  b.  Then  (x  -  x  )  e  Null  (,V) ,  so  there  exist  constants 


Y  =  (Y, 


,YA-r) 


such  that 


0 


■  — .  A— r 

>  Y.s.  = 

^i=l  1  1 


T 

=  yS 


Conversely,  any  vector  of  form  (2-18)  satisfies  the  equation  xA  =  b.  I 
2.6  Conclusion 


A  capsule  summary  of  the  results  reported  in  this  section  may  be 
stated  as  follows.  First,  certain  previously  reported  results  relating  to 
the  algebraic  consistency  of,  and  the  structure  of  solutions  to,  the  Kosut 
gain  equation  are  substantially  improved.  Moreover,  conditions  for  design- 
model-  independence  of  free  parameters  occurring  in  such  solutions  are  given. 
Second,  a  clearer  idea  of  the  nature  and  extent  of  performance  degradation 
to  be  expected  due  to  decentralizing  the  feedback  information  structure  is 
obtained.  Finally,  stable  designs  obtained  by  the  extended  Kosut  method 
are  shown  to  be  robust. 

However,  some  substantial  open  questions,  both  in  theory  and  in 
application  methods,  must  be  resolved  before  the  extended  Kosut  approach 
can  become  a  competitive  method  for  large  space  structure  controller  design. 

Major  theoretical  problems: 

1.  Characterize  systems  for  which  the  extended  Kosut  approach  leads 
to  a  stable  closed-loop  design  (at  the  design  model  level). 

2.  Synthesize  a  systematic  procedure  for  assigning  values  to  free 
parameters  (when  available)  so  as  to  improve  performance  of  the  full-order 
system  incorporating  the  reduced-order  controller. 
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Major  applications  problems 

1.  Develop  the  algorithm  for  solving  the  Kosut  gain  equation  [8] 
into  an  effective  design  tool. 

2.  Compare  the  extended  Kosut  method  with  several  other  represen¬ 
tative  design  methods  on  more  sophisticated  structural  models. 
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SECTION  3 


STABILITY  MARGINS  FOR  FIXED-ORDER  COMPENSATORS 


3.0  Abstract 


In  order  that  a  theory  can  be  developed  for  establishing  an  upper 
bound  on  closed-loop  stability  when  evaluating  a  fixed-order  compensator 
with  higher-order  evaluation  models  of  large  space  structures,  the  closed- 
loop  dynamic  equations  were  derived  first  (see  Subsection  3.1).  As  an 
initial  effort,  the  common  observer-based  fixed-order  controllers  were 
considered;  the  design  of  such  compensators  was  assumed  to  be  based  only 
on  a  reduced-order  model  formed  by  truncating  all  but  a  fixed  number  of 
vibration  modes.  The  resultant  closed-loop  dynamic  equations  were  then  put 
into  the  standard  second-order  matrix  form  in  two  different  ways,  all  having 
a  certain  form  of  effective  damping  and  stiffness  matrices. 

The  closed-loop  effective  stiffness  matrices  were  not  symmetric; 
their  symmetrization  was  then  studied,  so  that  standard  stability  conditions 
on  second-order  dynamics  could  be  applied.  Two  steps  were  proposed: 
symmetrizing  transformation  of  the  coordinates  and  orthonormalizing 
transformation  of  the  new  eigenvectors.  (See  Subsection  3.2.)  The  closed- 
loop  dynamic  equations  were  then  put  in  a  modal  form  similar  to  the  open- 
loop  case.  Determining  the  closed-loop  stability  and  its  margin  would 
then  be  a  straightforward  matter  after  carrying  out  this  two-step  procedure, 
say,  by  a  computer.  Analytic  expressions  for  the  closed-loop  stability 
conditions  (with  higher-order  evaluation  models),  however,  were  not  easy 
to  derive. 

The  stability  of  the  closed-loop  system  having  a  fixed-order  compen¬ 
sator  could  be  looked  at  directly  in  the  frequency  domain.  In  other  words, 
the  shortest  distance  of  the  closed-loop  poles  to  the  imaginary  axis  of  the 
complex  plane  could  determine  the  upper  bound  of  closed-loop  stability. 

An  example  was  worked  out  analytically  in  detail  (see  Subsection  3.3)  to 
demonstrate  that  it  would  be  very  difficult  to  obtain  an  analytical 
expression  of  closed-loop  stability  margins  from  this  viewpoint  even  for 
a  structure  as  simple  as  having  only  two  elastic  modes. 
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3. 1  Closed-Loop  Dynamic  Equations  with  Observer-Based  Fixed-Order 
Compensators 

Assume  that  an  open-loop  modal  representation  of  the  large  space  struc¬ 
tures,  in  the  form 


h  +  (G  +  D)*t  +  ft  n  =  $u 
y  =  Pn  +  Vn 


(3-1) 


is  available,  and  that  the  critical  modes  are  the  low-frequency  ones  while  the 
residual  modes  are  of  high  frequencies.  G  and  D  are  gyroscopic  and  damping 
matrices  respectively;  ft  is  a  diagonal  matrix  of  natural  frequencies.  Let  i, 
m,  n,  n  ,  and  denote  the  numbers  of  outputs,  inputs,  modes,  critical  modes, 

and  residual  modes,  respectively.  Consider  the  following  partitioning. 
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We  presume  that  the  parameters  of  the  full  model  are  available,  but  that  the 
design  is  to  be  based  essentially  on  the  critical-mode  part  of  the  model.  The 
feedback  controls  that  will  be  considered  are  of  the  form 


u  =  Kn  +  Ln 
c  c 


(3-4) 
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where 


n  +  (G  +  D  )n  +  fin 
c  cc  cc  c  c  c 


8  u  +  H(y  -  P  n  -  Vn) 
c  c  c  c  c 


(3-5) 


In  (3-5) ,  n  has  the  form  of  an  observer  for  n  •  Note  that  it  is  driven  by 
c  c 

the  output,  y,  which  contains  nr»  and  is  thus  not  actually  an  observer  of  nc 

in  the  formal  sense.  Evidently,  (3-5)  involves  only  the  critical  mode  param¬ 
eters.  The  free  design  parameters  are  the  gain  matrices  K,  L,  and  H,  and  we 
want  to  know  if  they  can  be  chosen  by  a  procedure  so  that  (3-1),  (3-4),  (3-5) 
yield  a  stable  closed-loop  system,  particularly  when  G  and  D  are  very  small. 


The  reason  for  assuming  this  form  is  that  it  is  representative  of  a 
range  of  practical  schemes  for  finding  gains.  Note,  for  instance,  that  K,  L, 
could  be  designed  on  the  basis  of  an  LQ  control  problem,  and  that  H  could  be 
designed  on  the  basis  of  a  Kalman  filtering  problem  or  as  a  model  matched 
filter  for  the  critical  modes.  However,  these  methods  have  not  been  shown  to 
lead  in  general  to  closed- loop  stability,  and  more  refined  conditions  on 
their  application  are  needed. 

The  standard  coordinates  for  the  closed-loop  dynamics  are  n  and 


=  nc  " 


(3-6) 


We  write  these  in  second-order  form  rather  than  state-space  form.  Let 


Then 


x  + 


G+D+HV  0  HV+G+dS 

cc  cc  c  r  cr  cr  I 


Effective  damping 


(3-7) 
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3.2  Symmetrization  of  Closed-Loop  Stiffness  Matrices 

We  note  that  the  closed-loop  system  can  be  put  into  second-order  form, 
but  it  is  not  clear  whether  this  could  be  transformed  to  an  equivalent  de¬ 
coupled  system 


+  [G  +  D]  x  + 


a2  i 


=  0 


(3-11) 


where 

~T  ~T  ~2 

G  =  -G  ,  D  =  D  ,  and  0  is  diagonal. 

If  it  could,  then  a  rather  simple  stability  test  for  second-order  systems 

_2 

could  be  applied.  For  example  [1],  [2],  when  Q  is  positive  definite,  system 
(3-11)  is  asymptotically  stable  if  and  only  if  the  following  associate  system 


z  +  Gz  +  h2  z  =  0 

—  • 

y  =  Dz 


is  completely  observable. 

The  conversion  of  (3-10)  to  (3-11)  proceeds  in  two  steps: 

Step  1 :  A  symmetrizing  transformation,  S,  which  converts  (3-10) 
so  that  its  last  term  is  symmetric. 

Step  2:  An  orthonormal  transformation,  $,  which  puts  the  result  of 
Step  1  into  modal  form. 

Step  1 

Suppose  we  perform  a  change  of  variable  S  as  follows: 

Sv  =  x  (S  nonsingular)  (3-12) 

For  simplicity,  rewrite  system  (3-10)  in  the  form 

x  +  Fx  +  Kx  =  0  (3-13) 

then 

v  +  S_1FSv  +  S_1KSv  -  0  (3-14) 

If  the  last  term  is  to  be  symmetric,  there  must  be  an  S  such  that 
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-1-  T-T  -IT 

S  KS  =  STS 


(3-15) 


T-T  T  -1 

i.e.,  K  =  (SS  )K  (SS  ) 

or  equivalently,  there  must  be  a  nonsingular  symmetric  transformation  T  such 
that 


KT  =  TKT 


(3-16) 


In  other  words,  such  a  transformation  must  be  solved  from  the  following  mat¬ 
rix  equation: 


-Hp  Tn  T12 


2  T 

-6K  Q  T12  T22 


T11  T12  -KB 

T  T  T  2 

T12  T22  -PH  0 


(3-17) 


where 


i  fi  -  e  K  +  HP 
c  c  c 

To  complete  Step  1,  suppose  that  such  a  symmetric  nonsingular  transformation 
T  has  been  obtained.  Then  there  is  a  nonsingular  square  root  of  T,  denoted 
S,  such  that 


SS  =  T 


(3-18) 


This  can  be  found  via  Cholesky  decomposition. 


Step  2 


Define 


K  =  S'^S 


(3-19) 


A  =  S^FS 


(3-20) 


Also  define 


&- 

ON 

H 

V-/ 

(3-21) 

i(A  + 

«T 

A  ) 

(3-22) 

Then  there  exists  an  orthogonal  transformation  4>  such  that 

4>TK4>  =  h2  (3-23) 

Defining 

G  =  <5>TGi>  (3-24) 

D  =  4>TG0  (3-25) 


brings  the  system  into  the  form  (3-11). 

Comment  1 :  If  steps  1  and  2  can  be  carried  out,  it  is  then  possible 
to  interpret  the  closed-loop  system  in  terms  of  modes ,  i.e.,  as  if  it  repre¬ 
sented  a  new  structure.  It  was  not  previously  apparent  that  "closed-loop 
modes"  could  be  defined  in  the  sense  that  open-loop  ones  were. 

Comment  2:  If  any  of  the  steps  can  be  carried  through  analytically, 
then  explicit  stability  conditions  could  be  given.  This  has  not  yet  been 
accomplished . 

Comment  3:  This  procedure  might  represent  a  more  simple  way  to 
determine  closed-loop  stability  numerically  using  existing  computer  programs 
for  high-order  systems. 

Comment  4:  Since  the  design  value  of  the  foregoing  results  often  de¬ 
pend  crucially  on  the  ability  to  obtain  analytical  expressions  involving  the 
gains,  several  examples  were  worked  in  order  to  obtain  a  better  qualitative 
feel  for  how  the  gains  enter  the  exact  analytical  expressions  for  stability. 
In  symmetrization  of  a  two-mode  system  with  control  of  one  critical  mode,  we 
carried  out  step  1  of  the  symmetrization  procedures,  Eq.  (3-17),  assuming 
n^  =  n^_  =  1.  The  expressions  obtained  were  very  complex.  It  should  be  ap¬ 
parent,  however,  that  the  problem  of  finding  a  symmetrizing  transformation 
can  be  expected  to  yield  interesting  stability  results. 

3 . 3  Closed-Loop  Stability  Conditions  for  a  Two-Mode  Example 

The  stability  problem  of  system  (3-10)  can  be  viewed  directly  in  the 
frequency  domain.  In  this  case,  we  seek  the  roots  of  the  polynomial 
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M(s) 
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(C 

cc 


+  D 

cc 


tf  L  +  HV  )s  +  (  2 

c  c  c 


c  K  +  HP  ) 

C  L' 


det 


-H(Vs  +  P) 


-p(Ls  +  K) 


j  -) 

s  •*•((.  +  D)s  + 


Note  that  the  lower  right  block  is  just  the  polynomial  of  the  open-loop 
system,  which  is  stable  by  assumption.  Thus  a  stabilitv  analysis  might  be 
based  on  the  perturbation-theory  expression 


M(s)  =  det  ^  s2  +  (G  +  D  -  p  L  +  HV  ) s 
|  recce  c 

+  (,2  -  8  K  +  HP  ) 
c  c  c 

-  H (Vs  +  P)  [S2  +  (C  +  D)s  +  .2]  t (Ls  +  K) 
•  det  [s2  +  (G  +  D)s  +  ii  ] 


An  effort  has  been  made  to  work  out  M(s)  for  a  two-mode  example,  and  to  apply 
Routh's  criterion  so  that  exact  stability  conditions  might  be  found  for  the 
simplest  case  of  interest.  The  results  follow 
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cr  cr 


where  lower-case  letters  have  been  used  to  denote  scalar  quantities.  Note 
that  g  =  =  0  by  the  requirement  that  gyroscopic  coupling  be  antisym¬ 

metric.  Also  d  =  d  and  g  =  -g  have  been  used.  Subscript  c  denotes 
rc  cr  °rc  cr  r 

"critical"  and  subscript  r  denotes  "residual"  as  above.  Lengthy  algebra 
yields 

6  . 

M(s)  =  £  Mhs 

i-0 


A3 


where 
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Forming  the  Routh  array  for  a  general  6th-order  polynominal, 
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The  most  simple  special  case  assumes 

1.  No  gyroscopic  effects  g  =  0 

2.  No  off-diagonal  damping,  d  =  d  =0 

cr  rc 

3.  Velocity  sensing  only  p^  =  p  =  0 

4.  Feedback  only  on  velocity  estimate,  k  =  0 
In  this  case 
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It  is  still  impractical  to  work  out  the  entire  Routh  array  analytically,  but 
we  note  that 

M,  =  l  -  0  (first  column  must  be  >0  for  stability) 

6 

Mq  •  0  (last  element  of  first  column  is  positive) 

Mr  0  implies  v  h  -  g  ?  >  -  (2d  +  d  ) 

5  r  c  c  cc  rr 
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This  shows  that  the  third  entry  in  the  first  column  of  the  Routh  array  will 
be  positive  under  the  more  restrictive  assumptions  that 

h  >  0 

l  <  0 

6  v  -  3  v  >0  (more  "control  authority"  on  critical  than 
c  c  r  r 

residual  modes) 

At  this  point,  our  algebraic  capacity  was  exhausted,  although  it  is  clear  that 
additional  conditions  can  be  derived. 

Conclusion:  Since  in  practice  d  ,  d  are  very  small,  we  have  in  ef- 

-  cc  rr 

feet  shown  that  all  but  three  conditions  of  Routh's  test  can  be  satisfied 
under  conditions  of  a  typical  example,  by  choosing  l  <  0,  h  >  0.  It  is  un¬ 
likely  that  the  remaining  conditions  will  allow  any  less-restrictive  condi¬ 
tion  than  this. 
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SECTION  4 


AGGREGATE  MODEL  DEVELOPMENT 


4.0  Abstract 


A  common  practice  in  obtaining  a  reduced-order  model  of  a  large  space 
structure  for  designing  its  vibration  controllers  is  to  truncate  the 
vibration  modes  directly,  according  to,  say,  the  line-of-sight  motions. 

Such  design  models  do  not  account  for  the  influences  of  actuators  on  the 
truncated  modes  nor  the  influences  of  the  truncated  modes  on  the  sensors. 

The  attendant  control  spillover  and  observation  spillover  can  cause  the 
vibration  controllers  to  destabilize  the  structures,  even  when  these 
controllers  have  been  carefully  and  satisfactorily  designed  to  ensure  the 
closed-loop  stability  with  respect  to  the  reduced-order  design  models. 

The  development  of  aggregate  design  models  was  thus  intended  to  provide 
better  reduced-order  models  of  the  structures  for  use  in  the  design  of 
vibration  controllers.  Vibration  modes  can  be  aggregated  with  respect  to 
actuator  and  sensor  influences  by  a  nonsingular  transformation  before 
making  any  truncation;  the  truncated  portion  can  then  be  independent  of  the 
actuators  and  sensors.  Such  an  aggregation  transformation  in  the  general 
form  was  carefully  studied  (see  Subsection  4.1).  A  useful  set  of  necessary 
and  sufficient  conditions  and  two  separate  sets  of  sufficient  conditions 
for  the  aggregation  transformation  to  be  nonsingular  were  established; 
such  conditions  are  useful  as  a  guide  in  choosing  specific  values  of  a 
diagonal  block  of  free  parameters  in  the  aggregation  transformation. 


The  aggregation  transformation  is  not  unique.  Several  interesting 
ideas  of  utilizing  the  uniqueness  (i.e.,  of  selecting  the  specific  values 
of  various  free  parameters  in  the  transformation)  for  imposing  additional 
properties  on  the  aggregate  design  models  were  examined  (see  Subsection  4.3); 
some  of  their  potential  difficulties  in  obtaining  closed-form  analytical 
results  were  analyzed. 

Vibration  modes  can  be  aggregated  also  with  respect  to  their 
influences  on  performance  and  the  influences  of  disturbances  in  addition  to 
actuator  and  sensor  influences.  For  example,  the  control  systems  for  a 
large  space  structure  might  be  designed  to  minimize  the  transmission  from 
the  disturbance  sources  to  the  critical  performance  variables  besides 
providing  the  structures  with  closed-loop  stability.  By  including  the 
influences  of  the  disturbances  as  a  "part"  of  actuator  influences  and, 
similarly,  the  influences  on  the  performance  as  a  "part"  of  sensor  influences, 
this  "generalized"  aggregation  problem  was  able  to  be  recast  into  the  same 
initial  form  (see  Subsection  4.3),  so  that  the  same  arguments  and  the  same 
techniques  could  be  applied  and  extended  directly. 
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4.1 


Aggregation  of  Vibration  Modes  with  Respect  to  Actuator  and 
Sensor  Influences 


A  standard  approach  to  model  reduction  is  to  truncate  the  open-loop 
modal  representation  directly,  by  neglecting  (generally  high-frequency)  non- 
critical  natural  modes.  This  generally  results  in  control  energy  being  coupled 
into  the  neglected  portion  in  a  way  which  can  be  destabilizing.  Here  we  explore 
a  different  choice  of  basis.  Rather  than  using  the  open-loop  modes  directly, 
we  aggregate  them  first  based  on  the  influence  of  the  actuators  and  sensors. 

The  advantage  is  that  the  coupling  with  the  neglected  portion  becomes  indepen¬ 
dent  of  feedback  control  gains,  and  is  thus  much  more  readily  predicted. 
Furthermore,  this  choice  of  basis  is  also  not  affected  by  feedback. 

Consider  the  following  usual  finite-element  modal  representation  of  LSS. 

( q  +  ft2q  =  $TB„u 

I  F 

|  y  «  cDtn  +  cv»n  (4-1) 


where  n  denotes  an  L-vector  of  normal  coordinates;  ft  =  {to^,  ....  ui^}  and  $  - 
[<(>.,  ...»  4>  ]  are  L  x  L  matrices;  oi.  and  <j> .  denote  the  natural  frequency  and 

1  iw  11 

the  natural  shape  of  mode  i,  respectively.  L  is  usually  a  very  large  integer, 
n  *  . ..,  n^)  denotes  the  corresponding  normal  coordinates,  u  =  (u^,  .... 

u  )  denotes  the  control  inputs  to  m  force  actuators,  and  B„  the  L  x  m  influence 
m  r 

matrix,  y  =  (y^,  ...,  y^)  denotes  the  observation  outputs  from  displacement 

sensors  and  velocity  sensors  (0  £  i,  0  <_  <_  &);  Cq  and  denote  the 

corresponding  l  *  L  influence  matrices.  Superscript  T  denotes  transpose.  Nor¬ 
mally,  m  <  L  and  l  <  2L. 

A  corresponding  state-space  form  is 


c  =  rcD<i>  cv<t] 
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We  assume  that  both  B  and  C  are  of  full  rank.  The  rows  of  matrix  B  and  the 
columns  of  matrix  C  can  be  rearranged  such  that 


B 


(4-3) 


with  B^  being  an  n  *  m  matrix  of  rank  m  and  an  l  x  n  matrix  of  rank  l  for 

the  smallest  integer  n.  We  assume  the  corresponding  rearrangements  have  been 
made  on  the  rows  and  columns  of  matrix  A. 


Now,  to  aggregate  the  state  variables  of  (4-2),  namely,  the  vibration 
modes  of  (4-1),  with  respect  to  actuator  and  sensor  influences,  we  aggregate 
the  rows  of  matrix  B  and  the  columns  of  C  by  a  similarity  transformation  T 
such  that 


tc1 ;  o] 


(4-4) 


where  B^  is  also  an  n  *  m  matrix  of  rank  m  and  also  an  l  *  n  matrix  of  rank 
l.  By  defining 

x  =  Tx ,  A  =  TAT_1,  B  -  TB,  C  =  CT_1  (4-5) 

we  now  have  the  following  special  but  equivalent  representation  of  the  LSS 

*  Ax  +  Bu 


I1  ' 

|y  -  cx 


(4-6) 


Partitioning  these  matrices  conformably  with  (4-4)  yields 


A  = 


l~All 

!  A12 

B  = 

V 

! - - 

A21 

a 

22 

_  0  . 

[c1  0] 


where  ,  in  particular,  is  an  n  *  n  matrix.  The  corresponding  state  equation 
can  thus  be  rewritten  as 


AnXl  +  A12X2  +  B1U 


A21xl  +  A22X2 


C1X1 


(4-7) 


Using  any  such  aggregation  transformation  (and^the  orthogonal  trans¬ 
formation  representing  the  rearrangement  of  matrices  A,  B,  and  C  before  ag¬ 
gregation),  we  have  aggregated  the  vibration  modes  with  respect  to  the  influ¬ 
ence  of  actuators  and  sensors.  Observe,  specifically,  that  the  control  of  (4-7) 
does  not  directly  affect  nor  is  the  observation  directly  affected  by  it. 

Neglecting  we  8et  the  following  aggregate  design  model  of  order  only  n: 


x  =  A^x  +  BjU 


C1X 


(4-8) 


Roughly  speaking,  such  a  reduced-order  model  represents  an  equivalent 
dynamics  of  those  linear  combinations  of  natural  modes  which  can  be  directly 
measured  and  controlled.  It  is  not  an  exact  aggregation  [1]  -  [3]  of  system 
(4-2)  in  the  sense  of  Aoki.  By  proper  selection  of  the  aggregation  transforma¬ 
tion  T,  this  reduced-model  may  be  made  to  approximate  open-loop  performance 
of  system  (4-1),  much  closer  than  any  one  obtained  by  direct  truncation  of  the 
vibration  modes. 


Comment  1 :  In  general,  it  is  sufficient  to  have  £  +  m  rows  in  the  block  B^  and 

£  +  m  columns  in  the  block  C ^ ;  hence  n  <_  £.  +  m.  Depending  on  the  specific 

form  of  matrices  B  and  C  at  hand  have,  or  due  to  the  existence  of  both  velocity 
and  displacement  measurements,  the  blocks  B^  and  need  only  max{£,m}  rows 

and  columns  respectively;  hence  n  £max{£,m}.  Note  that  the  possibility  of 
having  an  aggregate  design  model  (4-8)  of  order  less  than  £  +  m  is  an  improve¬ 
ment  over  the  preliminary  results  [7]. 


Comment  2:  The  rearrangement  of  matrices  B  and  C  to  have  them  bear  the  form 
described  by  (4-3)  can  be  carried  out  using  common  computational  procedures  such 
as  Gaussian  elimination  or  Gram-Schmidt  orthogonalization  on  rows  of  B  and 
columns  of  C.  Any  prior  effort  in  the  placement  of  actuators  (and  sensors), 
or  in  the  synthesis  of  actuator  (and  sensor)  influences  for  preventing  control 
spillover  to  (and  observation  spillover  from)  secondary  modes,  as  discussed  in 
[4-6],  will  facilitate  the  rearrangement  of  these  matrices,  since  there  will 
be  fewer  nonzero  rows  in  B  (and  fewer  nonzero  columns  in  C)  to  work  with. 

Since  the  number  of  actuator-influenced  (and  sensor-influenced)  modes,  partic¬ 
ularly  those  which  are  secondary  modes,  will  be  decreased,  such  spillover 
prevention  efforts  will  also  facilitate  the  aggregation  and  reduce  the  coupling 
Aj £  and  A2i‘ 
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Comment  3:  The  aggregation  of,  say,  matrix  B  as  described  in  (4-4)  is  equiv¬ 
alent  to  making  proper  linear  combinations  of  the  rows  of  B  to  produce  the 
specified  form.  In  conjunction  with  representation  (4-7),  this  can  be  inter¬ 
preted  as  synthesizing  the  actuator  influences  so  that  control  inputs  will  not 
spill  over  part  of  the  aggregated  state  x.  Compare  this  interpretation  with 

the  idea  [4-6]  of  preventing  control  spillover  to  secondary  modes  by  synthesis 
of  actuator  influences.  To  prevent  control  spillover  by  synthesis  is  to  make 
proper  linear  combinations  of  columns  of  matrix  B  to  achieve  a  partitioning 
similar  to  what  is  described  by  (4-4),  namely,  to  determine  a  transformation  T 
such  that 


Br 


It  follows  from  this  comparison  that  the  technique  developed  for  computing 
synthesis  transformation  T  can  be  modified  to  help  compute  the  aggregation 
transformation  T.  This  comment  is  equally  applicable  to  the  aggregation  of 

matrix  C  with  respect  to  an  aggregation  transformation  S  =  T  observation 
spillover,  and  synthesis  of  sensor  influences. 

Let  the  transformation  be  partitioned  with  respect  to  (4-3)  and  (4-4) 
as  follows 


The  dimension  of  the  blocks  are  as  follows: 

Tu>  n  x  n;  T^,  n  x  (2L  -  n);  T21>  (2L  -  n)  *  n;  T^,  (2L  -  n)  x  (2L  -  n) . 
Then  Eq.  (4-4)  can  be  written  as 


T11B1  +  T12B2 

—* 

B1 

(4-9) 

T21B1  +  T22B2 

= 

0 

(4-10) 

C1T11  = 

(4-11) 

C1T12  = 

S 

(4-12) 

is  of  the  full  row  rank,  its  right 

inverse  CB  exists. 

Therefore,  from 

(4-12),  we  have 
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(4-13) 


where 


T12  '  ClS 


R  T  T  —  1 

C1  =  C1<C!C1>  +N12 


and  is  an  arbitrary  matrix  such  that 


C1N12  =  0 


Similarly,  has  its  left  inverse  and,  from  (4-10),  we  have 


where 


T21  ~  T22B2B1 


B1  5  +  S21 


and  N21  is  an  arbitrary  matrix  such  that 


N21B1  “  ° 


It  follows  from  (4-14)  and  (4-16)  that 


T  = 


r  t.  .  cRc„' 

*1  o' 

p 

11  12 

-t„„b„b^  T„„ 

0  T„„ 

22  2  1  22 

22 

11 

*2B1 


Therefore,  the  aggregation  transformation  T  is  nonsingular  if 
T22  and  |tli  +  are  nonsingular. 


Comment 


_4:  The  specific  value  of  matrices  and  in  (4-4) 


should  be  consequences  of  the  aggregation  transformation  T. 
contrary  that  C.  has  been  prescribed  before  the  aggregation. 
(4-11) 


11 


cRc 

cri 


(4-14a) 


(4-14b) 


(4-15) 


(4-16a) 


(4-16b) 


and  only  if  both 

,  in  general. 

Suppose  on  the 
Then,  from 
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Consequently,  the  matrix  sum 


11 


Cld2®2=l 


:i(£l  +  e2626l) 


cannot  be  nonsingular  unless  for  such  a  very  special  case  that  n  =  i  >_  m  and 
+  C2B2B^  is  nonsingular.  Therefore,  it  is  better  not  to  prescribe  the 
value  of  before  aggregation. 

Matrix  should  not  be  prescribed,  either.  Premultiplying  (4-9)  by  Cj^ 
and  using  (4-11)  -  (4-12)  and  (4-3),  we  have 

C1B1  =  +  C2B2  H  CB  (4-17) 


This  relationship  in  general  will  not  be  satisfied  if  B  is  assigned  a  specific 
value  before  aggregation. 

In  the  sequel,  we  assume  that  both  T^  and  T22  are  nonsingular,  and  that 

both  0^  and  B^  are  computed  after  aggregation.  Now,  from  (4-13)  and  (4-15),  we 
have 


lll 

0 


01 


22j 


-®2S1 


T~ 

1irr2 


Theorem  1;  The  aggregation  matrix  T  is  nonsingular  if  and  only  if  the 


—  1  R"  ~  -  L  \ 

n  x  n  matrix  (I  +  T^CjC2B2Bj)  or  the  (2L  -  n)  x  (2L  -  n)  matrix 


|l  +  B2B^Tj  |c^C2  j  is  nonsingular,  where 


■  eiTu 


(4-18) 


Theorem  2:  The  following  are  two  separate  sets  of  sufficient  conditions 
for  the  matrix  T  to  be  nonsingular. 

(a)  B2  -  0  or  C2  =  0 

T  ~T 

(b)  TiiTn  =  <1  =  ra,  and  C  *  EB  for  some  nonsingular  matrix  E 
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Comment  5:  Condition  (a)  appears  to  be  trivial,  since  either  the  actuator  in¬ 
fluences  or  the  sensor  influences  have  been  aggregated  properly.  This  condi¬ 
tion,  however,  clearly  indicates  the  advantage  of  using  Gram-Schmidt  or 
Gaussian  elimination  procedure  for  rearranging  the  rows  of  B  and  the  columns 
of  C,  and  the  advantage  of  prior  effort  made  in  preventing  control  and  obser¬ 
vation  spillover  by  placement  or  synthesis;  see  Comment  2. 

Comment  6:  First  part  of  condition  (b)  means  that  T  is  a  special  kind  of  non¬ 
singular  matrix  that  is  orthogonal;  the  identity  matrix  obviously  is  such  an 
orthogonal  matrix.  The  second  and  third  parts  mean  that  there  are  as  many 
sensors  as  there  are  actuators  and  that  the  sensor  influences  as  a  whole  are 
similar  to  the  actuator  influences  as  a  whole.  For  LSS  having  only  force 
actuators,  this  implies  that  no  displacement  sensors  are  considered  for  aggre¬ 
gation;  namely,  C^  =  0  in  (4-1).  The  special  case  in  which  E  =  kl  represents 

the  commonly  considered  case  where  each  force  actuator  is  collocated  with  a 
velocity  sensor,  and  vice  versa.  Notice  that  vibration  modes  with  such  a 
special  case  can  always  be  aggregated  to  achieve  the  desired  representation  (4-7) . 

Comment  7 :  The  expressions  (4-13)  -  (4-16)  for  T^  and  T2i  are  slightly  dif¬ 
ferent  from  those  given  in  [7].  These  slight  improvements,  however,  have  made 
it  possible  to  obtain  useful  new  results,  such  as  Theorems  1  and  2. 

4 . 2  Ideas  for  Specifying  the  Nonunique  Aggregation  Transformation 

The  similarity  transformation  T  for  accomplishing  the  aggregation  as  des¬ 
cribed  by  (4-4)  is  not  unique;  subject  only  to  some  rather  minor  constraints, 
submatrices  T^,  Nj2>  anc*  N21  are  virtuaHy  free  parameters.  Additional 

properties  may  thus  be  imposed  on  the  aggregate  design  model  (4-8)  by  specifying 
these  parameters.  The  following  are  some  interesting  ideas  we  made  a  quick 
look  at. 


4.2.1  Minimization  of  Coupling  Coefficients  ,,  and  A^ ^ 

Unlike  the  reduced-order  design  model  usually  obtained  from  direct 
truncation  of  vibration  modes,  the  aggregate  design  model  (4-8)  has  the  desirable 
property  that  the  coupling  between  the  retained  part  x^  and  the  truncated  part 

x^  remains  unchanged  even  in  the  closed  loop  with  feedback  control  systems. 

This  constancy  does  not  depend  on  how  the  control  systems  are  designed,  be  they 
designed  as  LQ  regulators  and  observers,  or  as  dynamic  output  feedback  com¬ 
pensators  [7].  Therefore,  it  would  be  very  desirable  if  the  open-loop  coupling 
coefficients  A^  and  A£^  could  be  minimized  by  a  suitable  choice  of  the  non¬ 
unique  aggregation  transformation. 

Partition  the  aggregated  matrix  A  defined  in  (4-5)  and  the  original 
matrix  A  conformably  with  the  similarity  transformation  T.  Then  we  have 
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■’I 


A11  A12  T11 

T12  T11  T12  A11  A12 

a21  A22  T21 

T22_  _T21  T22_  A21  A22_ 

Solving 

for  A12 

and  A^  yields 

> 

ro 

If 

[T11A12 

+  T12A22  "  (TllAll 

+  *12*21  >Ti}T12](T22  -  Iji^Jlij)'1 

(4-19) 

A21  " 

[T21A1 1 

+  t22a21  -  <t21a12 

+  T22A22)T22T2l]<Tll  '  T12T22T21>'‘ 

(4-20) 

In  terms  of  and  A^,  we  have 


4l  (T1LA11  +  T12A21  "  A12T2I)T1I 


A22  =  (T21A22  +  T22A22  ~  A21T12)T 


-1 

22 


Substituting  (4-13)  and  (4-15)  in  (4-19)  and  (4-20), 
A12  =  Tll(A12  +  Tll( 


-  a  T-^r^r 
'12  11  11  1^2 


-  t11c1c2a21t11c1c2)(i  +  b2b1th  c1c2)  t22 


21 


+  A2i  - 


!2  (~^2^1A1 1 
A22B2B1 ) (r  +  T1 1 C 1 C  2B  2B2 )  T11 


(4-21) 

(4-22) 


(4-23) 


(4-24) 


With  the  expressions  (4-14a)  and  (4-18),  the  terms  within  the  parentheses  of 
(4-23)  and  (4-24)  are  virtually  independent  of  T ^  and  T22.  Thus,  to  minimize 

A^2  is  to  minimize  and  maximize  T22>  whereas  it  is  the  opposite  for  mini¬ 
mizing  A,^.  Therefore,  the  coupling  coefficients  k^^  and  A2j  cannot  be  mini¬ 
mized  simultaneously  by  any  choice  of  T^  and  T22< 
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The  only  possibility  lies  in  the  selection  of  the  null-space  matrices 
and  involved  respectively  in  the  right  inverse  and  the  left  inverse 

B^.  The  terms  in  the  first  parenthesis  of  (4-23)  can  be  minimized  by  a  separate 
choice  of  N^«  while  those  of  (4-24)  can  be  minimized  simultaneously  by  a  separ¬ 
ate  choice  of  N,^.  The  matrix  products  in  the  second  parenthesis  of  (4-23)  and 
(4-24)  can  be  maximized  simultaneously  by  a  joint  choice  of  and 


The  constancy  of  the  terms  in  the  parentheses  of  both  (4-23)  and  (4-24) 
can  be  seen  more  easily  if  we  consider  the  class  of  orthogonal  matrices  for 
T^,  namely 


For  this  class 


T~VRr 

illLlL2 


=  CRC 

C1C2 


(4-25) 


which  is  independent  of  T^.  All  the  terms  within  these  parentheses  are 
completely  independent  of  T^! 

Nevertheless,  the  numbers  of  nonzero  columns  in  and  nonzero  rows  in 

A^^  can  be  reduced.  Yeung  Yam  of  MIT,  while  expanding  the  preliminary  results 

in  [7],  has  recently  derived  a  specific  choice  of  the  free  parameters  in  T 
such  that,  for  the  case  n  =  i  +  m,  only  n  columns  of  A^  and  n  rows  of  A 
are  nonzero. 

4.2.2  Minimization  of  Differences  in  Closed-Loop  Transfer  Functions 

Feedback  control  systems  are  to  be  designed  based  only  on  the  reduced- 
order  aggregate  design  model  (4-8).  The  resultant  closed-loop  transfer  functions 
for  the  design  model  alone  generally  will  be  different  from  the  resultant  ones 
for  the  overall  system  representing  the  LSS.  It  would  be  very  desirable  if 
the  differences  could  be  minimized  by  suitably  specifying  the  free  parameters 
of  the  aggregation  transformation. 

As  in  [7],  consider  fixed-order  dynamic  compensators  of  the  form 

|  z  =  Fz  +  Gy 

| u  =  Hz  +  Ky  +  u  (4-26) 

where  the  compensator  state,  z  has  a  fixed  dimension  which  is  usually  much 
less  than  2L;  u  denotes  the  command  inputs;  F,  G,  H,  and  K  are  design 
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parameters  to  be  determined  based  on  the  reduced-order  model  (4-8).  Substitu¬ 
ting  (4-26)  in  (4-8)  yields  the  following  closed-loop  system  of  the  design  model. 


where 


f-  •  -1 

z 

z~ 

'  0  ■ 

x . 

II 

>1 

L*1. 

+ 

B, 

1 

1 

(4-27) 


11 


GC, 


BjH  An  +  B  KCr 


(4-28) 


The  transfer  function  matrix  of  the  closed  loop  system  (4-27)  between  u  and  y 
is  found  to  be 


Q(s) 


[0  Cj]  (si  -  An) 


lr0n 
B1 


=  C, 


si  -  A u  -  B1KC1  -  BjlKsI  -  F)"1  GC1  1  Bx  (4-29) 


Now  substituting  (4-26)  in  the  overall  system  (4-7)  yields  the  following 


closed-loop  system. 


(4-30) 


where 


A  = 

A 

An 

_A'l2_' 

ss 

F 

B1H 

gc1  ; 

A11  +  BlKcli 

0 

A12 

A21 

a 

22 

0 

A21  ! 

A22  - 
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with  the  same  A^ 
closed  loop  system 


Q(s) 


as  defined  by  (4-28). 
(4-30)  between  u  and 


1  1 


B 


KC. 


B  ll(sl  - 


The  i  rails!  er  !  unction  matrix  of  the 
v  is  found  to  be* 


A12(hI 


A  t 


B 


1 


Define  the  difference  by 


(4-32) 


E(s)  =  Q(s)  -  Q(s) 


Then,  one  could  attempt 
that 


to  choose  the  parameters  T^,  T77,  N  ,  and  ^21 


(4-33) 

such 


OO 


— oo 


is  minimal,  where  s  =  o  +  jj  and  s  denotes  the  complex  conjugate  of  s.  To  ob¬ 
tain  a  solution  to  this  minimization  problem  is  impra ttical ,  if  not  impossible, 
judging  from  the  algebraic  complexity  involved  when  the  expressions  (4-19)  - 
(4-22)  for  A.,  and  (4-13)  -  (4-16)  for  T.~  and  T~.  are  substituted  in  (4-29) 
and  (4-32).  1J 

Intuitively,  the  main  difference  between  Q(s)  and  Q(s)  is  the  term 
A12(sl  -  A22^  A21”  T*ie  minimizat3i-on  °f  this  term  will  also  minimize  the  over¬ 

all  difference  as  defined  by  (4-33)  or  (4-34).  See  the  preceding  subsection 
(Subsection  4.2.1)  for  discussion  on  minimizing  the  coupling  coefficients  A^7 
and  A21. 

Another  intuitive  alternative  is  to  maximize  the  terms  within  the  brack¬ 
et  of  (4-29)  so  that  the  corresponding  terms  in  (4-32)  dominate  the  term 

A ! 2  (si  “  a22^  A21 ‘ 

These  two  intuitive  approaches  may  be  much  easier  to  take  than  direct 
minimization  of  the  overall  difference  defined  by  (4-34)  and  any  other  way  using 
(4-33).  The  necessary  algebra  still  can  be  very  complicated. 

4.2.3  Assurance  of  Closed-Loop  Stability  of  LSS 

Again,  feedback  controllers  are  to  be  designed  using  only  the  aggregate 
design  model  (4-8),  and  at  least  to  assure  closed  loop  stability  with  respect 
to  the  design  model.  It  would  be  desirable  if  closed  loop  stability  could 
still  be  assured  when  such  controllers  were  implemented  on  LSS,  which  contains 
not  only  the  design  model  but  also  the  truncated  part.  Such  an  assurance  is 
generally  impossible  for  the  case  of  direct  truncation  of  vibration  modes.  Re¬ 
call  that  the  coupling  A^  and  A21  of  the  aggregate  design  model  are  independent 

of  any  controllers  designed  by  any  method.  Thus,  if  closed-loop  stability  of  the 
LSS  having  a  controller  designed  for  the  aggregate  design  model  using  some 
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model  using  some  method  can  be  guaranteed,  then  closed-loop  stability  of  the 
LSS  having  any  such  controllers  but  designed  by  other  methods  can  also  be 
guaranteed . 


Consider  again  the  feedback  controllers  of  the  form  (4-26).  Assume  the 

design  parameters  F,  G,  H,  and  K  are  computed  based  on  (4-8)  so  that  the  closed- 

loop  system  (4-27)  is  stable.  The  question  is  whether  it  is  possible  to  find 

(and  specify)  among  the  nonunique  aggregation  transformation  one  such  that  the 

closed-loop  system  (4-30)  is  also  stable.  This  requires  the  use  of  the  stability 

theory  for  interconnected  large  scale  subsystems;  but  the  development  of  such 

stability  theory  has  just  begun  [8]  and  no  results  are  directly  applicable  to 

(or  general  enough  for)  our  problem  at  hand  yet.  Moreover,  with  the  expressions 

for  A..,  and  T  .  substituted  in  (4-31),  the  resultant  matrix  A  is  a  very  com- 
ij 

plicated  nonlinear  matrix  function  of  parameters  T^ ,  T22’  ^12’  anc*  ^21'  ®e~ 

riving  an  analytical  stability  condition  for  system  (4-30)  is  very  difficult, 
if  not  impossible,  at  present.  An  attempt  was  made,  however,  in  applying 
Willems'  generalization  [9]  of  Lyapunov  theory  to  system  (4-30). 

4 . 3  Aggregation  of  Vibration  Modes  also  with  Respect  to  Tracking  and 
Disturbance  Rejection 

In  the  previous  derivations,  it  was  implicitly  assumed  that  the  pur¬ 
pose  of  the  fixed-order  compensator  (4-26)  was  merely  to  enhance  the  stability 
of  the  closed-loop  system.  The  aggregation  was  thus  based  only  on  the  actu¬ 
ator  and  sensor  influences.  Now  assume  that  there  are  a  finite  number  of 
disturbance  sources  on  the  structure  and  a  finite  number  of  critical  response 
variables  (such  as  LOS  motion  in  x  and  y  directions).  Then  the  compensator 
might  be  designed  for  minimizing  the  transmission  from  the  disturbance  sources 
to  the  critical  response  variables  for  all  critical  frequencies.  This  implies 
that  the  compensator  places  particular  emphasis  on  damping  out  those  vibration 
modes  which  can  influence  the  critical  response  variables  and  can  be  excited 
by  the  disturbance  sources.  For  such  a  generalized  purpose,  aggregation  of 
the  vibration  modes  should  also  include  the  influences  of  disturbance  sources 
and  critical  response  variables. 

We  extend  the  preceding  studies  in  the  following  way.  Suppose  the 
original  representation  is  now  expanded  as  follows: 

Ax  +  Bu  +  Dv 
Cx  +  Ew 
Mx  +  Nu 


(4-35) 


fc>0 


where 


v  is  a  vector  of  disturbance  inputs 
w  is  a  vector  of  observation  noises 

r  is  a  vector  of  response  (not  necessarily  measurable) 
u,  x,  y  are  as  before. 

Now  define 


B  =  [B  D  0],  D 


~o 

0 

E  ' 

u 

II 

<  U 
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_N 
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-  w_ 

C  = 
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,  y  = 

_  M  _ 

_ r  _ 

,  A  =  A,  x  =  x 


x  = 


Xy.  +  Bu 


y  = 


-A  A  /V  A 

Cx  +  Du 


(4-36) 


Now  by  arguing  as  before,  there  is  a  transformation  which  brings  the  repre- 

a  —l  -  -  -i 

sentation  (4-35)  into  the  form  (A,B,C,D)  =  (TAT  ,  TB,  CT  ,  D) ,  where 


A  = 


A  ,  A,  „ 

B, 

11  12 

1 

,  B  = 

A  ,  A„  „ 

0 

21  22 

_ 

C  = 


(4-37) 


here  Aj^  has  the  dimension  between  max  {dim  u,  dim  y}  and  (dim  u  +  dim  y) . 

The  representation  (4-36)  is  now  transformed  into  a  similar  form  to  (4-7). 
By  making  the  same  truncation,  we  obtain  an  aggregate  design  model  which 
is  also  in  a  similar  form  to  (4-8). 


Now  suppose  we  have  a  compensator 


z  =  Fz  +  Gy 
u  =  Hz  +  Ky 
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Then  the  closed-loop  dynamics  is  given  by 


GC 


*  *  * 
SB  A  +  B  KC 


r  =  [NH  M*  +  NK.C* ] 


0  GE 


*  * 

D  B  KE 
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z 

V 

+  [0  NKE] 

X 

-  ~  /V  *  *  *  * 

where  N  =  N,  E  =  E  (since  D  =  D),  and  B  ,  C  ,  D  ,  M  are  identified  by  parti¬ 
tioning  B  and  C,  i.e.. 
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SECTION  5 


PARETO-OPTIMAL  CONTROL 


5.0  Abstract 


Canavin's  "modal-dashpot"  concept  and  Aubrun's  "low-authority 
control"  theory  were  two  off-the-shelf  examples  used  to  illustrate  the 
multiple-objective  optimization  problems  underlying  the  design  of  vibration 
controllers  for  large  space  structures  (see  Subsection  5.1).  Canavin  used 
pseudo-inverses  to  compute  the  output  feedback  gains  required  for  augmenting 
a  certain  level  of  damping  to  primary  modes  (or,  specifically,  critical 
modes) ;  the  computed  gains  were  of  dismally  high  values.  A  careful  exam¬ 
ination  showed  that  such  a  high-gain  problem  need  not  exist  because  there 
were  free  parameters  that  could  be  adjusted  to  minimize  the  magnitude  of  the 
required  gains.  It  also  showed  that  the  desire  for  large  possible  achievable 
damping  and  for  small  possible  required  gains  could  easily  be  formulated  as 
a  two-objective  optimization  problem  and  solved  directly  and  efficiently  by 
multiple-objective  methods  without  computing  the  pseudo  inverses  and  the 
appropriate  values  for  the  free  parameters.  In  Aubrun's  theory,  the  feed¬ 
back  gains  were  to  minimize  a  cost  functional  that  was  the  sum  of  two 
disparate  design  objectives:  the  magnitude  of  the  errors  between  predicted 
and  desired  closed-loop  damping,  and  the  magnitude  of  feedback  gains.  The 
tradeoff  between  the  closed-loop  damping  errors  and  the  feedback  gains,  and 
the  variation  of  the  former  versus  the  latter,  could  be  systematically  and 
efficiently  made  by  formulating  the  underlying  two-objective  optimization 
problem  and  then  applying  a  multiple-objective  optimization  method. 

The  applicability  and  application  of  two  multiple-objective  opti¬ 
mization  methods  were  demonstrated  by  Pareto-optimal  designs  of  modal 
dashpots  for  the  12-mode  Draper  model  (see  Subsection  5.2).  Not  only  a 
sequence  of  the  two-objective  Pareto  optimal  designs  but  also  the  maximum 
damping  achievable  for  each  level  of  feedback  gains  were  computed  system¬ 
atically.  Insights  into  the  movement  of  the  closed-loop  poles  as  the 
feedback  gains  increase  their  magnitude  were  easy  to  obtain  from  the  resulting 
root  loci.  Systematic  evaluation  for  robustness  of  the  resulting  sequence 
of  designs  to  parameter  variations  and  to  spillover  became  easy.  The 
variations  of  the  sensitivity  (i.e.,  of  the  robustness)  were  also  obtained. 
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5.1 


Applicability  of  Multiple-Objective  Optimization 
Methods  to  the  Design  of  Vibration  Controllers 

The  design  of  a  control  system  normally  has  more  than  one  performance 
objective  to  achieve:  e.g.,  large  modal  damping,  low  control  response,  little 
small  feedback  gains,  small  errors  or  overshoots,  fast  response,  little 
spillover.  As  a  standard  practice,  a  designer  usually  combines  the  multiple 
objectives  into  one  single  performance  index  using  some  rather  arbitrary 
weights.  The  designer,  usually  not  satisfied  with  the  resultant  control 
systems,  then  repeatedly,  aimlessly,  and  sometimes  endlessly,  changes  the 
weighting  coefficients  to  obtain  new  designs.  Such  frustration  is  not 
uncommon  in  the  design  of  "modern  control  systems"  using  the  linear- 
quadratic  regulator  theory;  the  weighting  matrices  Q  and  R  frequently  are 
changed  over  and  over  again.  Experience  with  artificial  single-objective 
optimization  used  on  multiple-objective  design  problems  has  often  left  the 
designer  dissatisfied  with  the  lack  of  information  (or  systematic  information) 
on  the  available  tradeoffs  between  his  multiple  objectives  and  his  "multiple- 
objective  optimal"  designs.  Realistically,  incompatible  objectives  cannot 
be  unified;  important  multiple  objectives  should  be  addressed  directly  and 
their  trade-off  be  made  efficiently  and  systematically. 

Pareto  optimality  is  a  practical  notion  of  "multiple-objective" 
optimality  that  is  useful  to  control  systems  designers  faced  with  multiple 
disparate  performance  objectives.  Stated  simply,  a  design  is  Pareto  optimal 
if  none  of  the  multiple  objectives  can  be  further  improved  without  degrading 
any  other.  Multiple-objective  optimization  by  the  method  of  proper  equality 
constraints  (PEC)  or  the  method  of  proper  inequality  constraints  (PIC)  help 
the  designers  obtain  both  Pareto  optimal  designs  of  control  systems  and 
efficient  tradeoffs  of  the  multiple  objectives  directly  and  systematically. 

Consider  the  following  two  different  approaches  to  the  design  of 
vibration  controllers  for  large  space  structures:  Canavin's  "modal- 
dashpot"  concept  [1-4]  and  Aubrun's  "low  authority  control"  theory  [5,6], 

Let  us  point  out  their  underlying  multiple-objective  design  problems  and 
briefly  demonstrate  the  applicability  of  the  above-mentioned  multiple- 
objective  methods  [7-10]. 

Assume  that  large  space  structures  are  represented  in  modal  coordinates 
as  follows: 


••2  T 

n  +  n  n  =  $  b_,u 

r 


y  =  cv<J>n 


(5-1) 

(5-2) 
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where  ft  =  diag  , . . .  ,gj^  | ,  and  $  *  [$^,...,$^1  are  L*L  matrices:  and 

4> ^  denote  respectively  the  natural  frequency  and  shape  of  the  ith  normal 

mode  of  vibration,  n  =  denotes  the  corresponding  normal 

coordinates,  u  =  (u,,...,u  )  denotes  the  control  inputs  to  m  force  actu- 

1  m 

ators,  and  B  denotes  the  Lxm  influence  matrix,  y  =  (y ,,..., y .)  denotes 
F  IX- 

the  observation  outputs  from  l  velocity  sensors:  denotes  the  correspond¬ 

ing  £xL  influence  matrix.  Superscript  T  denotes  transpose. 


In  general,  a  very  large  number  of  vibration  modes  are  included  in 
this  model.  Let  {u)pi>  4>p^}t  *  =  1»***»P»  denote  a  selected  subset  for 

suppression;  we  call  them  primary  modes  for  a  self-explanatory  reason. 
Completely  neglecting  all  nonprimary  modes  yields  the  following  truncated 
form 


np  +  ^np  =  opbfu 


y  =  Vpnp 


(5-3) 

(5-4) 


where  ftp  and  <tp  are  similarly  defined  in  terms  of  the  p  primary  modes, 

and  np  =  (lpl , , , . .hpp)  denotes  the  corresponding  normal  coordinates. 

Vibration  control  systems  are  then  designed  for  this  model  as  if  it  had 
exactly  modeled  the  LSS  in  question. 

Both  Canavin's  and  Aubrun's  approaches  are  for  designing  constant-gain 
velocity-output  feedback  control  in  the  form 

u  =  -Gvy  (5-5) 


Substituting  (5-5)  in  (5-3)  and  (5-4)  yields  the  following  closed-loop 
system  dynamics 

np  +  $pbfgvcvVp  +  =  0  (5"6) 


Consider  Canavin's  modal-dashpot  concept  first.  In  order  that  each 
primary  mode  may  be  damped  independently  and  arbitrarily,  the  feedback  gain 
matrix  Gv  is  computed  using  pseudo- inverses  as  follows: 


(*X)T[*X(*>F)Tr  (2W  [(Vp)T  CV*p]'‘  (VpJ 


(5-7) 
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where  Zp  denotes  a  diagonal  matrix  of  desired  damping  ratios  £  for  the  p 

primary  modes.  It  is  assumed  that  S,  >  p  and  m  i  p.  Each  actuator  was 
assumed  to  be  collocated  with  a  velocity  sensor  in  his  numerical 
example  [1,2].  For  10%  damping  for  each  primary  mode,  the  computation  using 
(5-7)  resulted  in  a  gain  matrix  of  shockingly  large  values.  He  felt  that 
"modal  dashpots  may  be  of  limited  utility  due  to  the  high  gains  produced 
by  this  approach"  [1,2]. 

Such  a  high-gain  problem  actually  need  not  exist.  For  exactly  the 
same  purpose  of  independently  and  arbitrarily  damping  the  primary  modes. 


the  gain  matrix  G^  can  be 

T 

computed  using  a  right  inverse  of  $  B_  and  a 

r  r 

left 

inverse 

of  Cylp,  instead. 

as  follows: 

cv 

’  (4pBr)R  (2zA)(Vp)L 

(5-8) 

where 

(*WR  ■- 

(*pbf)T  [♦Jbf(*p»f)T]"1  +  »f 

(5-9) 

(cvV)L  5 

[(cv*p)T  Vp]"‘  (Vp)1  +  Nv 

(5-10) 

with  Np 

and  Nv  being  arbitrary  null-space  matrices  such  that 

=  0 

(5-11) 

=  0 

(5-12) 

(See  Section  6.2,  particularly  Comments  7  and  8,  concerning  a  generalization 
and  an  extension  of  the  modal  dashpot  concept.)  The  similarity  and  the 
difference  between  (5-7)  and  (5-8)  are  obvious.  Whenever  the  rank  of 
T 

<I>pBp  is  less  than  the  number  m  of  actuators,  matrix  Np  need  not  be  zero. 
Similarly,  whenever  the  rank  of  Cy^p  is  less  than  the  number  l  of  sensors, 
matrix  N„  need  not  be  zero,  either.  Nonzero  N  and  N  serve  as  free  para- 
meters  and  can  be  used  to  reduce  the  magnitude  of  gain  matrix  G^.  The  unde¬ 
sirable  consequence  of  neglecting  these  adjustable  free  parameters  is  thus  clear. 

Now  that  we  have  recognized  the  existence  and  advantage  of  free 
parameters  Nr  and  N  ,  an  obvious  question  is  how  to  choose  their  values 

to  minimize  the  magnitude  of  the  required  feedback  gains  while  using  the 
feedback  gains  to  augment  as  much  damping  to  the  primary  modes  as  possible. 
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The  underlying  two  performance  objectives  now  become  evident,  and  are 
formulated  as  a  bi-objective  optimization  problem  as  follows.  Note  that 
in  the  interest  of  improving  model  fidelity,  we  also  retain  in  the  control 
design  model  some  other  important  modes  which  we  call  secondary  modes  and 
denote  by  subscript  S.  Active  damping  of  secondary  modes  is  not  required, 
but  prevention  of  control  spillover  to  and  observation  spillover  from  them 
is  desired.  Both  primary  and  secondary  modes  are  altogether  referred  to 
as  modeled  modes  and  denoted  by  subscript  H.  Therefore,  the  bi-objective 
modal-dashpot  design  problem  is  to  find  the  feedback  gains  g^  and  the 


damping  ratios 
maximize 


Ci 


that 


and  minimize 


fl  =  £  «  C. 

1  -1- 


i=l 


Pi 


m  i 
i=l  j=l 


(5-13) 


(5-14) 


subject  to 


VfWh  ■  2ZA 


c,  ;  c. 


for  primary  modes 


’i  "  Pmin 

L.  1  C„  .  for  secondary  modes 
i  bmin 


(5-15) 

(5-16) 

(5-17) 


where  weighting  factors  w^  are  useful  for  expressing  desired  mode  movement 

(see  Subsection  5.2).  Observe  that  no  computation  of  the  pseudo- inverses , 
the  right  or  the  left  inverses,  or  the  free  parameters  is  required.  They 
are  all  embedded  in  the  equality  constraint  (5-16)  of  closed-loop  modal 
decoupling.  The  method  of  proper  equality  constraints  (PEC)  is  to  set 
the  first  objective  to  a  parametric  level,  i.e.. 


i=l 


Wi^Pi 


=  a 


(5-18) 


whereas  the  method  of  proper  inequality  constraints  (PIC)  is  to  bound  the 
first  objective  by  a  parametric  level,  i.e., 

P 

L  *  a  (5-19) 

i-1  1  Pi 
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Whether  the  method  of  PEC  or  PIC  is  used,  the  bi-objective  problem  (5-13)  - 
(5-17)  is  all  converted  to  a  parametric  single-objective  quadratic 
programming  problem  of  the  following  standard  form:  minimize 

J  =  xTQx  (5-20) 


subject  to 


Ax  ^  b 
x  >  0 


Many  algorithms  and  computer  programs  have  been  developed  for  solving  the 
quadratic  programming  problem  (5-20);  see,  for  example,  References  [11-15]. 

For  each  selected  value  of  the  parameter  ot,  the  resulting  quadratic 
programming  problem  is  solved  for  the  single-objective  optimal  solutions. 

Using  the  test  for  Pareto  optimality  as  given  in  References  [7-10],  the 
Pareto  optimal  values  of  feedback  gain  g  and  damping  ratios  are 

obtained  systematically  in  the  increasing  order  of  a.  The  resulting 
minimized  values  of  the  objective  (i.e.,  the  magnitude  of  gain  matrix  G^) , 

plotted  versus  parameter  a,  will  help  the  designer  make  optimal  tradeoffs 
efficiently  and  systematically  between  achievable  damping  and  minimum 
required  gains.  See  Subsection  5.2  or  Refs.  [16-17]  for  more  details. 

It  is  worth  mentioning  that  the  free  parameters  N  and  N  are 

r  V 

implicitly  used  to  minimize  the  magnitude  of  required  feedback  gains  and  to 
decouple  secondary  modes  from  primary  modes  in  the  closed  loop.  Such  was 
not  possible  with  Canavin's  origin  modal-dasnpot  concept.  On  the  other  hand, 
there  is  no  question  at  all  about  the  applicability  of  multiple-objective 
optimization  methods  to  design  of  modal  dashpots. 

Now,  consider  Aubrun’s  low-authority-control  (LAC)  theory.  According 
to  Aubrun  [5,6],  for  sufficiently  "small"  feedback  gains  g„  ,  the  closed- 

loop  damping  of  mode  i  is  predicted  using  the  following  formula  [5,6]: 

dA.  *  2^.  *  ^Wi  (5"21) 

using  the  open-loop  frequency  oo^  and  mode  shape  4^.  To  compute  the  gains  g„  , 
let  (dA^)p  denote  the  predicted  closed-loop  damping  given  by  the  above 
formula,  and  let  (dX^p  denote  the  desired  closed-loop  damping  imposed  by 
the  LAC  controller  design.  Then  the  gains  g„  are  chosen  so  as  to  minimize 
the  quadratic  cost  functional 

J(CV)  .  E  [(dA±)D  -  (d*.)p]  2  +  E  g2j  (5-22) 
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in  which  the  modal  weights  w^  help  specify  pole  locations  and  the 

term  £  .  improve  robustness  of  the  controller  [6].  It  is  not  difficult 

to  see  from  (5-22)  that  two  disparate  objectives,  one  concerning  damping 


and  the  other  concerning  gains,  are  "unified", 
optimization  problem  is  to 


The  underlying  bi-objective 


minimize 


«Vd  -  «Vp 


(5-23) 


and  minimize 


=  £ 


i,  j 


’ij 


(5-23) 


Pareto  optimal  solutions  of  this  problem  by  the  method  of  PEC  or  PIC  will 
provide  a  curve  of  J~  versus  for  systematically  making  efficient  tradeoff 

between  achievable  damping  errors  and  the  corresponding  feedback  gains 
required.  The  first  objective,  J^,  is  converted  to  a  parametric  equality 

constraint  if  the  method  of  PEC  is  used,  while  it  is  converted  to  a 
parametric  inequality  constraint  if  the  method  of  PIC  is  used  instead.  All 
converts  the  bi-objective  problem  to  a  simplest  form  of  nonlinear  programming 
problem,  where  the  single  objective  function  and  the  parametric  constraint 
are  quadratic.  Many  nonlinear  programming  algorithms  and  computer  programs 
already  developed  [14,15]  can  be  used  in  solving  the  resultant  constrained 
optimization  problem.  Applicability  of  the  multiple-objective  optimization 
methods  is,  again,  not  a  problem  at  all. 

Comment  1 :  Both  underlying  bi-objective  design  problems  (5-13)  - 
(5-17)  and  (5-23)  -  (5-24),  have  the  minimization  of  feedback  gains  as  an 
objective;  such  is  a  very  practical  engineering  objective  in  addition  to 
the  design  objective  on  damping. 


Comment  2 :  The  closed- loop  damping  expressed  by  formula  (5-21)  may 
not  be  valid  since  the  closed-loop  modes  may  not  remain  decoupled  with  the 
same  open-loop  frequencies,  even  with  respect  to  the  reduced  design  model. 
Only  when  the  gains  g^  are  extremely  small  will  the  closed-loop  modes  be 

virtually  uncoupled.  On  the  other  hand,  the  desired  independence  can  be 
explicitly  imposed,  such  as  (5-15);  the  resultant  Pareto  optimal  designs 
then  will  certainly  ensure  that  the  open-loop  modes  remained  decoupled  at 
least  in  the  closed  loop  of  the  design  model. 


Comment  3:  The  minimization  of  feedback  gains,  such  as  expressed  by 
(5-14)  or  (5-24),  does  not  actually  provide  robustness  of  the  controllers 
designed.  It  will  reduce  control  and  observation  spillover,  however, 
because  the  residual  modes  will  be  less  excited  and  less  observed. 
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5 . 2  Application  to  Draper  Tetrahedral  Model  of  LSS 

The  12-mode  Draper  tetrahedral  model  of  LSS  was  used  as  an 
example  [16,17].  Pareto  optimal  designs  of  modal  dashpots  and  the  tradeoff 
curve  were  sequentially  generated  by  a  bi-objective  optimization  program 
developed.  In  search  of  more  or  less  uniform  movement  of  the  primary  modes 
under  efficient  tradeoffs  of  damping  and  gain  objectives,  we  found  the 
form  of  inverse-sum-of-squared-inf luence  (ISS1)  for  the  weighting  factors  w^ 
in  (5-13)  more  satisfactory  than  others. 

Each  of  the  resulting  designs  was  then  examined  for  its  robustness 
to  model  parameter  variation  and  spillover  (resulting  from  modal  truncation). 
For  all  designs  with  velocity  sensors  collocated  with  force  actuators, 
robustness  of  stability  to  parameter  variation  and  to  spillover  was  clearly 
demonstrated.  The  damping  performance  is  relatively  insensitive  to 
parameter  variation  and  to  spillover;  however,  the  sensitivity  rapidly 
increases  as  the  desired  level  of  critical-mode  damping  increases.  In  other 
words,  modal  dashpots  can  theoretically  be  designed  to  provide  any  desired 
level  of  damping  to  the  critical  vibration  modes.  In  practice,  however, 
only  low  designed  levels  of  modal  damping  can  be  closely  achieved  due  to 
this  sensitivity  to  parameter  variation  and  spillover. 

For  those  designs  where  the  sensors  and  actuators  were  not  all 
collocated,  stability  and  damping  performance  are  also  relatively  insen¬ 
sitive  to  parameter  variation.  However,  spillover  causes  all  designs  t< 
destabilize  the  originally  stable  structure;  at  least  half  of  the  resultant 
closed-loop  poles  have  positive  real  parts.  Destabilization  by  modal 
dashpots,  which  are  also  referred  to  as  robust  controllers,  has  not  been 
demonstrated  before.  Consequently,  collocation  of  sensors  with  actuators 
is  essential  to  the  robustness  of  modal  dashpots  when  used  as  vibration 
dampers  for  realistic  large  space  structures. 
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SECTION  6 


GENERAL  REQUIREMENTS  FOR  STABILIZATION 


6.0  Abstract 


Feedback  from  velocity  sensor  outputs  can  augment  damping  and 
feedback  from  displacement  sensor  outputs  can  augment  stiffness  to  large 
space  structures.  The  design  of  the  feedback  gains  are  normally  based  only 
on  modeled  modes;  the  exisi_ence  of  unmodeled  modes  is  completely  ignored. 

Due  to  the  attendant  control  and  observation  spillover,  even  stability  of 
the  closed-loop  system  having  feedback  controllers  thus  designed  may  not  be 
ensured,  let  alone  •'he  additional  damping  and  stiffness  desired.  Useful 
general  conditions  in  various  useful  forms  were  recently  established  (see 
Subsection  6.1)  for  guaranteeing  closed-loop  asymptotic  stability.  Such 
conditions  are  useful  as  a  guide  to  the  design  of  closed-loop  stabilizing 
feedback  gains.  Large  space  structures  considered  were  not  limited  to 
those  having  only  stable  elastic  modes;  rigid-body  modes  and  unstable  modes 
were  not  excluded.  For  those  special  structures  having  only  stable  elastic 
modes  and  rigid-body  modes,  useful  simple  conditions  were  also  established. 
The  conditions  for  robustness  to  parameter  variation  and  to  modal  truncation 
were  also  established. 

Canavin's  modal-dashpot  concept  was  generalized  and  extended  (see 
Subsection  6.2).  The  similar  concept  of  "modal  springs"  was  therefore 
proposed.  The  general  formulas  for  computing  the  feedback  gains  of 
generalized  modal  dashpots  and  modal  springs  were  derived.  These  formulas 
contain,  various  adjustable  free  parameters.  The  closed-loop  asymptotic 
stability  ensurance  and  the  robustness  to  parameter  variation  and  to  modal 
truncation  of  modal  dashpots  and  modal  springs  were  established,  thus  at 
least  providing  a  theoretic  foundation  for  Canavin's  original  concept. 
Although  no  analytic  proofs  for  the  closed-loop  stability  ensurance  and 
robustness  of  modal  dashpots  were  given  by  Cnnavin  when  he  proDosed  the 
intuitive  concept,  a  numerical  demonstration  was  made  with  a  typical 
example  of  large  space  structures. 

Control  spillover  to  and  observation  spillover  from  secondary  modes 
can  be  prevented  before  the  beginning  of  the  controller  design  process  by- 
proper  placement  of  actuators  and  sensors  or  by  proper  synthesis  of 
actuator  and  sensor  iniiuences.  A  general-purpose  computer  program  is  being 
developed  for  synthesizing  actuator  influences  so  that  control  spillover 
to  secondary  modes  can  be  prevented  and  the  desired  control  influences  on 
primary  modes  can  be  obtained  simultaneously.  A  recent  theoretical  analysis 
has  led  to  the  conclusion  that  such  a  computer  program  can  be  used  to 
synthesize  sensor  influences  by  duality  and  even  to  determine  the  proper 
location  of  actuators  and  sensors  (see  Subsection  6.3). 
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6  .  1  General  Conditions  for  Closed-Loop  St ab  1 1 iy.at:  ion  by 
Output  Feedback  Control 

Consider  the  active  control  of  a  generic  large  space  structure  (LSS) 
which  may  contain  rigid-body  modes,  unstable  and  stable  elastic  modes.  As 
usual,  let  the  modal  equations  generated  by  a  finite-element  computer  program 
(such  as  NASTRAN)  be 


r 


T 

v 


(6-1 ) 


where  r,  denotes  the  L  (normal)  modal  coordinates  and  1  the  natural  mode  shape.-.; 
f  denotes  the  m  generalized  forces;  A  and  Z  are  diagonal  matrices  denoting, 
respectively,  the  modal  damping  and  the  modal  stiffness  coefficients.  In  gen¬ 
eral,  the  diagonal  elements  of  A  and  1  may  be  positive,  zero,  or  even  negative, 
and  can  be  extremely  small  if  not  zero.  Note  that  if  q  denotes  the  generalized 
physical  coordinates,  then 


q  =  ;  i  (6-2  ,i 

Consider  specifically  the  case  of  constant-gain  output  feedback  control,  namely, 
i  velocity  measurements  and  displacement  measurements  are  fed  back  to  m 

force  actuators  with  constant  gains  as  follows 


f 

=  BFu 

(6-3) 

u 

=  -“G  y  ~  G  v_. 

(fc-4 ) 

\rV  DM) 

>’v 

=  cvq 

(6-5a) 

VD 

(6-5b) 

In  other  words 


f 


BFGvcvq 


BF°DCDq 


(6-6) 


Substituting  (6-6)  in  (6-1)  yields 


T0  :)•  +  ( 

..  +  : tk  ; ] 

\ 

A  /  \ 

A  / 
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n  + 


0 


(6-  7 ) 


where 


da  j  bfgvcv 


ka  5  bfgdcd 


(6-8a) 

(6-8b) 


denote  the  additional  damping  and  stiffness  matrices  respectively.  The  design 
of  gain  matrices  and  should  make  both  the  augmented  damping  matrix 


(a  +  $TDA$j 


and  the  augmented  stiffness  mat 


rix  iy ! 


£  +  <t  K.^J  symmetric  and  posi¬ 


tive  definite,  in  addition  to  achieving  certain  design  objectives  (e.g.,  in¬ 
creasing  the  damping  of  critical  modes  to  desired  levels).  The  reason  is  sim¬ 
ple,  since  it  is  well  known  that  a  closed-loop  system  as  described  by  (6-7)  is 


asymptotically  stable  if  both  matrices 


(a  +  !TDa^)  and  |  L  +  j 


are  sym¬ 


metric  and  positive  definite,  and  closed-loop  stability  is  a  basic  requirement. 

Note,  however,  that  the  matrices  A  and  7  are  not  positive  definite  un¬ 
less  all  their  diagonal  elements  are  positive.  Also  note  that  both  matrices 
and  are  singular,  since  there  are  more  modes  represented  by  (6-1)  than 


there  are  actuators  or  sensors  placed  on  the  structure  (i.e.,  L 


V 


J- 


Moreover,  the  design  of  gains  G^  and  G^,  in  general,  can  be  based  only  on  a 

reduced-order  model  of  the  structure;  in  other  words,  onlv  a  truncated  version 
of  (6-1)  can  be  used.  Therefore,  the  asymptotic  stability  of  the  closed-loop 
system  (6-7)  does  not  follow  automatically  with  an  arbitrary  design  of  G^  and  G^. 

Let  subscript  M  denote  the  part  of  (6-1)  that  is  included  in  the  re¬ 
duced-order  model  for  control  system  design  and  subscript  U  the  unmodeled 
(i.e.,  truncated)  part.  Then  Eqs.  (6-1)  -  (6-5)  can  be  rewritten  as  follows 


nu  + 


Ju  'u 


i  +  Vvi 

= 

*mV 

(6-9a) 

[  +  V'u 

= 

<uV 

(6-9b) 

~Vv  ■  ' 

"DyD 

(6-10) 

CVr';MnM  + 

V 

unu 

(6-1  la) 

cpV'm  + 

V 

unu 

(6-1  lb) 

The  design  of  gains  G^  and  G^  ignores  the  existence  of  the  unmodeled  part  com¬ 
pletely,  and  is  to  make  only  the  modeled  part  achieve  desired  damping  and 
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stiffness  objectives.  Due  to  control  spillover  to  and  observation  spillover 
from  the  unmodeled  part,  as  can  be  seen  from  Eq .  (6-9b)  and  (6-11),  the  overall 
closed-loop  system  (6-7)  may  not  achieve  the  damping  and  stiffness  objectives. 
In  particular,  the  symmetry  and  positive  definiteness  of  the  overall  matrices 

|  A  +  4>  j  and  j  do  not  follow  automatically  from  only  making  their 

modeled  parts,  j  and  j ,  symmetric  and  positive  defi¬ 

nite.  We  have  recently  obtained  the  following  general  conditions  (Theorems  1 
through  5)  for  guaranteeing  the  symmetry  and  positive  definiteness  of  such 
overall  matrices  and  hence  for  guaranteeing  the  asymptotic  stability  of  closed- 
loop  system  (6-7).  Such  conditions  are  also  useful  as  a  guide  in  the  design  of 
gain  matrices  Gy  and  G^.  A  useful  special  case  is  discussed  in  Subsection  6.2. 


Theorem  1:  Assume  that  is  symmetric  and  |  is  posit 

definite.  Then,  for  (a  +  j  to  be  positive  definite,  it  is  necessary 

and  sufficient  that  ' 


ive 


Au  +  *uVu 


WmK  +  vme>a’m)  Wu 


(6-12) 


is  positive  definite. 

Theorem  2:  Assume  that  KA  is  symmetric  and  j is  positi 
definite.  Then,  for  | 1  +  to  be  positive  definite,  it  necessary  and 


ive 


and  sufficient  that 
L1 


T  T  /  T  \-l  T 

v  +  if  R  t  -  4.  K  it  '  +  I  K  ;  }  I  } 
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(6-13) 


is  positive  definite. 

Comment _ 1 :  For  closed-loop  asymptotic  stability  of  the  overall  system,  it  is 
not  enough  just  to  guarantee  the  positive  definiteness  of  the  modeled  parts 

( AM  +  Wm)  an<^  (  £'m  +  WM  )'  According  to  Theorems  1  and  2,  the  matrices 
Aj  and  2^  defined  by  (6-12)  and  (6-13),  respectively  must  be  made  positive 
definite  as  well  in  the  design  of  Gy  and  G^.  This  explains  why  current  ap¬ 
proaches  to  the  design  of  LSS  controllers  (e.g.,  direct  applications  of  state- 
of-the-art  control  techniques  only  to  the  modeled  part  without  any  considera¬ 
tion  of  the  unmodeled  part)  cannot  guarantee  closed-loop  asymptotic  stabilitv. 

Comment  2:  It  is  convenient  as  well  as  intuitively  desirable  to  be  able  to 

find  a  design  of,  say,  G  such  that  the  additional  damping  matrix  for  the 

V  T 

modeled  part  is  diagonal.  For  example,  using  the  pseudo- inverses  of  and 

Cf^M  a  diagonal  matrix  will  make  diagonal.  Notice,  however  that 

T 

those  Gy  designed  using  such  pseudo- inverses  will  make  only  symmetric;  to 


make  the  whole  matrix  4>  DA1>  symmetric  requires  that  each  velocity  sensor  be 

A  T 

collocated  with  a  force  actuator  such  that  C„  =  k,,B  for  some  nonzero  constant 

V  V  F 

ky.  Note,  however,  that  collocation  is  not  assumed  in  either  Theorem  1  or  2. 

Comment  3:  These  theorems  are  applicable  to  the  general  case  of  large  space 
structures  that  contain  not  only  stable  elastic  modes  but  also  rigid-body 
modes  and  unstable  modes.  To  provide  a  better  insight  into  the  requirements 
in  stabilization  of  general  large  space  structures,  refine  the  partitioning 
of  matrices  A  and  2  as  follows.  Denote  the  group  by  subscript  G  if  their 
diagonal  elements  of  A  are  greater  than  zero,  and  by  double  subscript  LE,  if 
otherwise.  Similarly,  denote  the  group 'by  G'  if  their  diagonal  elements  of 
2  are  greater  than  zero,  and  by  L'E',  if  otherwise.  Then  we  get  the  following 
similar  general  conditions 

Theorem  3:  Assume  that  D,  is  symmetric  and  (  A,  „  +  „  )  is  positive 

A  \  LE  LE  A  LL  J 

definite.  Then,  for  ( A  +  j  to  be  positive  definite,  it  is  necessary  and 

sufficient  that 


T 

<J>  D  i 

GAG 


T 

:  d  ,*> 

G  A  LE 


(* 


LE 


<p  D  t' 

LE  A  LE 


)-1 


T 

i>  D  4 

LE  A  G 


(6-14) 


is  positive  definite. 

Theorem  4:  Assume  that  Ks  is  symmetric  and  (  2, 

A  \LE  L.LALE/ 


positive  definite.  Then,  for  /  2  +  J 

necessary  and  sufficient  that  ' 


+  4>  K.$  )  to  be  positive  definite,  it  is 
A  1 


is 


ZG’  +  ^G,KA^G' 


*J'Vl 


’  E  ’  (‘ 


L'E' 


Ve'Vl 


♦l-E'Vc 


(6-15) 


is  positive  definite. 

Comment  4:  The  unstable  or  undamped  modes,  having  either  nonpositive  modal 
damping  or  nonpositive  modal  stiffness,  naturally  should  be  included  in  the 
model  for  designing  gain  matrices  G^7  and  G^,  and  at  least  both 

(“LE  +  ';lEBFGvVle)  and  \  L'E'  +  "l'E'Wd^L'E*  )  Sh°uld  be  made  Positive 
definite.  Theorems  3  and  4,  however,  show  that  this  is  not  sufficient  unless 
both  ^2  and  as  defined  by  (6-14)  and  (6-15),  are  guaranteed  to  be  positive 

definite.  In  general,  flexible  structures  having  small  positive  modal  damping 
(i.e.,  small  diagonal  elements  of  A  ) ,  and  lower  natural  vibration  frequencies 

(i.e.,  small  diagonal  elements  of  Z  ,)  will  be  difficult  in  providing  such 

Lj 

guarantees . 

Consider  the  following  common  special  cases:  only  rigid-body  modes,  and 
elastic  modes,  undamped  and  lightly  damped,  are  present.  In  the  above  partition¬ 
ing,  the  groups  denoted  by  LE  and  L'E'  now  correspond  to  only  zero. 
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diagonal  elements  of  A  and  I,  respectively.  Denoting  these  groups  now  by  E 
and  E',  we  get  the  following  two  very  useful  consequences  of  Theorems  3  and  4, 
concerning  closed-loop  stability  and  its  robustness. 

Theorem  5:  Assume  that  each  unmodeled  elastic  mode  has  some  positive 
damping,  no  matter  how  small  it  is.  Then  the  closed-loop  overall  system  (6-7) 

is  asymptotically  stable,  i.e.,  both  (a  +  C  fD  < )  and  (z  +  <fTK  C>)  are  symmetric 
and  positive  definite  if  '  /  V  A  / 

(a)  and  are  symmetric  and  at  least  positive  semidef inite,  and 

T  T 

(b)  and  are  positive  definite. 

Comment  5:  This  theorem  implies  that  the  closed-loop  asymptotic  stability  is 
robust  to  modal  truncation,  so  long  as  conditions  (a)  and  (b)  are  satisfied 
and  all  the  unmodeled  (i.e.,  the  truncated)  modes  have  some  positive  damping, 
finite  or  infinitesimal.  First,  unmodeled  modes  do  not  appear  in  conditions 
(a)  and  (b)  at  all,  control  spillover  and  observation  spillover  thus  have  no 
effect  on  closed-loop  asymptotic  stability  so  long  as  the  assumption  and  con¬ 
ditions  of  Theorem  5  satisfied.  Second,  closed-loop  asymptotic  stability  does 
not  depend  on  how  many  damped  elastic  modes  are  included  in  the  design  model 
nor  on  which  of  them  are,  since  neither  conditions  (a)  nor  (b)  depend  on  any 
damped  elastic  mode. 

Theorem  6:  Suppose  the  gain  matrices  C;  and  Gy  have  been  designed  to 

satisfy  conditions  (a)  and  (b)  of  Theorem  5.  LeL  there  be  parameter  variations 
(or  errors)  in  the  mass  matrix,  the  stiffness  matrix,  the  natural  frequencies, 
or  the  mode  shapes  of  the  structure  concerned.  Assume  that 

(i)  the  number  of  rigid-body  modes  and  that  of  undamped  modes 
remains  unchanged,  and 

(li)  4>e  =  $ECE»  *£«  =  ''e’CE' 

Where  the  overbars  denote  the  presence  of  parameter  variations  (errors);  C  and 

E, 

C  ,  are  nonsingular  matrices.  Then  the  closed-lo-.p  overall  system  remains 

t, 

asymptotically  stable. 

Comment  6;  Theorem  6  implies  that  closed-loop  asymptotic  stability  is  robust 
to  changes  in  the  shape  of  any  damped  elastic  modes,  modeled  or  not. 

The  following  subsection  (Subsection  6.2)  continues  the  discussion  of 
this  special  case  with  a  specific  design  of  gain  matrices  and  Gy  that  can 

provide  not  only  the  closed-loop  asymptotic  stability  but  also  its  robustness 
to  modal  truncation  and  parameter  variation.  Before  concluding  this  subsection, 
we  remark  that  the  controllability  and  observability  of  all  the  rigid-body 
modes,  all  the  unstable  modes  (with  either  negative  modal  damping  or  negative 
modal  stiffness  coefficients),  and  all  the  critical  elastic  modes,  are  assumed. 
The  actuators  and  sensors  should  therefore  be  properly  placed  or  their  existent 
influences  be  properly  synthesized,  according  to  the  techniques  and  conditions 
previously  derived;  see  Refs.  1-3;  Subsections  6.3  and  6.4. 


6.2 


Closed-Loop  Stabilization  with  "Modal  Dashpots"  and  "Modal  Springs" 


A  special  output  feedback  control  design  technique  is  presented  in  this 
section  for  integrated  control  of  elastic  and  rigid-body  modes,  for  increasing 
both  the  damping  and  frequency  of  selected  modes,  and  for  providing  closed- 
loop  asymptotic  stability  and  its  robustness  to  modal  truncation  and  to  param¬ 
eter  variation.  This  technique  represents  an  improvement  over  previous  work 
[4-8]  concerning  output  feedback  control  of  LSS  or  control  of  rigid-body  modes, 
and  is  an  advancement  along  Canavin's  modal-dashpot  design  philosophy  [4], 

Consider  again  the  common  special  case  of  large  space  structures: 
rigid-body  modes,  undamped  and  lightly  damped  elastic  modes  may  all  be  present. 
Assume  that  the  rigid-body  modes  are  not  to  be  ignored  as  was  in  [4]  but  to 
be  controlled  along  with  the  undamped  elastic  modes  and  those  lightly  damped 
elastic  modes  which  are  critical  to  the  damping  performance  of  the  LSS.  These 
modes  constitute  the  primary  modes.  Assume  that  the  design  model  is  formed  by 
these  primary  modes.  Therefore,  all  other  modes  are  unmodeled  modes  and  all 
are  assumed  to  have  some  positive  damping,  despite  how  small  they  are.  Using 
the  same  partitioning  as  described  at  the  end  of  the  foregoing  subsection,  we 
have  the  following  direct  consequences  of  Theorem  5. 


Theorem  7 :  Denote  the  number  of  modeled  (i.e.,  the  primary)  modes  by 
p.  Assume  that  th^re  are  at  least  as  many  force  actuators  as  there  are  modeled 
modes  (i.e.,  m  ^  p)  ,  and  that  each  force  actuator  is  collocated  with  a  velocity 

T  T 

sensor  and  a  displacement  sensor,  (i.e.,  Cy  =  kyBp,  CQ  =  kpBj,,  for  some  non¬ 
zero  constants  ky  and  k^) .  Assume  further  that  all  modeled  modes  are  com¬ 
pletely  controllable.  Then  the  closed-loop  system  is  asymptot ical ly  stable  if 
the  gain  matrices  Gy  and  GQ  are  computed  as  follows 

S  ■  K»f)R  <(Vm)L  (6-16> 


(6-17) 

(6- 18a) 

(6-1 9a) 

(6-20a) 


^  & 

where  A  and  £  denote  two  arbitrary  positive  definite  p  x  p  matrices  desired 
M  M 

as  the  additional  damping  and  stiffness  for  the  modeled  modes;  Np,  Ny,  and  NQ 
are  any  null-space  matrices  such  that 
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$mbfnf 

=  0 

(6—1 8b ) 

Vv  = 

(6- 19b) 

(6-20b) 

Comment  7 :  It  is  technically  convenient,  though  not  necessary,  to  consider 

*  -k 

A„  and  Z  as  diagonal  matrices.  Canavin  [4]  originally  called  the  resulting 

modally  decoupled  velocity  output  feedback  controllers  "modal  dashpots",  in 
which  the  rigid-body  modes  were  ignored,  no  displacement  feedback  was  consider- 

* 

ed,  both  Np  and  were  zero,  and  A^  was  diagonal.  Here,  we  extend  the  concept 

in  two  ways.  First,  any  velocity  feedback  controllers  designed  according  to 
(6-16),  (6-18)  and  (6-19)  are  called  (generalized)  modal  dashpots,  so  long  as 

Apj  is  diagonal.  Second,  any  displacement  feedback  controllers  designed  accord¬ 
ing  to  (6-17),  (6-18)  and  (6-20)  are  similarly  called  "modal  springs",  so  long  as 

•k 

the  matrix  Z^  is  diagonal.  It  should  be  noted  that  these  terms  (including  the 

one  initiated  by  Canavin)  and,  particularly,  the  implied  modal  independence  in 
the  closed  loop  make  sense  only  with  respect  to  the  design  model.  In  the 
presence  of  unmodeled  modes,  the  modeled  modes  may  not  maintain  their  designed 
independence  due  to  control  and  observation  spillover. 

Comment  8:  Since  the  gain  matrices  < >v  and  given  by  (6-16)  -  (6-20)  represent 

a  special  form  that  satisfies  conditions  (a)  and  (b)  of  Theorem  5,  the  asser¬ 
tions  on  robustness  of  closed-loop  asymptotic  stability  to  modal  truncation 
and  to  parameter  variation  remain  valid;  see  Comments  5  and  6,  and  Theorem  6. 

Comment  9:  Stability  and  robustness  results  for  the  special  subcase  [9]  in 

which  there  are  no  undamped  elastic  modes,  matrices  N  ,  N., ,  and  N  in  the  design 

r  V  Jl) 

*  ★ 

equations  (6-18)- (6-20)  are  all  zero,  matrices  A^  and  in  design  equations 

(6-16)- (6-17)  are  all  diagonal,  and  matrix  C  in  Theorem  6  is  diagonal  were  cb- 

tained  jointly  with  Dr.  Y.H.  bin  of  the  Jet  Propulsion  Laboratory,  California 
Institute  of  Technology  (who  was  formerly  with  the  C.S.  Draper  Laboratory). 
Various  numerical  demonstrations  were  also  made,  using  a  10-meter  free-free 
beam. 
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6 . 3  Prevention  of  Control  and  Observation  Spillover  by  Synthesis  of 
Actuator  and  Sensor  Influences 


Theorem  5  implies  that  for  LSS  having  no  unstable  modes  whatsoever  and 
no  undamped  unmodeled  modes,  control  and  observation  spillover  will  not  cause 
closed- loop  instability  so  far  as  output  feedback  control  is  concerned.  The 
design  of  the  output  feedback  controllers  can  be  based  entirely  on  the  modeled 
modes  which  include  all  rigid-body  modes,  all  undamped  elastic  modes,  as  well 
as  critical  elastic  modes.  Nevertheless,  the  closed-loop  performance  of  the 
vibration  controllers  (e.g.,  the  damping  and  frequency  shifting  of  primary 
modes)  will  be  degraded  if  control  spillover  and  observation  spillover  exist. 
Theorems  4  and  5  imply  that  control  and  observation  spillover  problems  may 
make  it  difficult  to  guarantee  closed-loop  stability  of  LSS  having  some  unsta¬ 


ble  modes. 


unless  the  control  spillover  matrices 


T 

*gbf 


and 


observation- 


spillover  matrices  C  <t  and  C  4>  ,  are  all  of  small  magnitude. 

V  0  1-)  u 


It  is  now  a 


standard  understanding  that  control  spillover  and  observation  spillover  will 
cause  closed-loop  instability  for  LSS  using  state-feedback  controllers  with 
state  estimators. 


Control  spillover  to  and  observation  spillover  from  secondary  modes  can 
be  prevented  by  proper  placement  of  actuators  and  sensors  on  the  structures 
and  by  proper  synthesis  of  actuator  and  sensor  influences  [10-12],  We  have  re¬ 
cently  concluded  that  the  technique  proposed  in  [10-12]  for  synthesizing  actu¬ 
ator  influences  can  be  extended  not  only  to  the  synthesis  of  sensor  influences 
(because  of  duality)  but  also  to  the  placement  of  actuators  and  sensors.  The 
formulas  derived  [10-12]  for  proper  placement  of  actuators  (and  sensors)  can 
be  recasted  into  similar  mathematical  forms  for  proper  synthesis  of  actuator 
(sensor)  influences.  A  general-purpose  computer  program  is  being  developed  for 
synthesizing  actuator  influences  so  that  control  spillover  to  secondary  modes 
is  prevented  and,  moreover,  desired  control  influences  on  primary  modes  is  ob¬ 
tained  . 


For  preventing  control  spillover  to  secondary  modes  by  synthesis  of 

T 

actuator  influences,  the  control  spillover  matrix  v  gB^  is  to  be  transformed 
using  column  operations  Q  and  row  operations  P  such  that 


W 

w 


11 

W12 

21 

W22 

P(*sVQ 


(6-21) 


T 

with  being  nonsingular  and  having  the  same  rank  at  'J’gBp.  Such  transforma¬ 

tion  (using  Q  and  P)  can  be  accomplished  by  Gaussian  elimination,  Gram— Schmidt 
orthogonalization,  or  singular-value  decomposition.  Preferably,  W^  is  made 

to  be  zero.  Then  the  synthesizer  F  for  prevention  of  control  spillover  to 
secondary  modes,  i.e.,  as  a  solution  of 
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0 


(6-22) 


T 

»sV 


is  given  by 


=  ( 


VuV’ 


(6-23) 


where  Q  =  [Qj  »  Q,,  ] .  The  matrix  F  represents  free  parameters  and  can  be  adjust¬ 
ed  for  achieving^desired  control  influences  on  primary  modes.  Let  B  denote 
a  desired  control-influence  matrix  for  the  primary  modes.  Then  V  is  to  be 
solved  from  the  following  equations 

(-VLW11W12  +  V  2) B  =  B  (6-24) 

where 

LV1’V2J  =  ^?Bf[Q1*Q2^  (6"25) 

The  program  under  development  is  first  to  carry  out  the  transformation 
(6-21)  using  Gaussian  elimination,  so  that  W  =  0.  The  transformation  (6-25) 

T  * 2 

of  actuator  influence  matrix  '1  B„  on  primary  modes  is  carried  out  simultan- 

r  r 

eously,  since  the  same  operations  Q  are  used.  This  program  is  then  to  compute 
f  from  (6-24).  A  desired  synthesizer  F  is  thus  obtained  from  (6-23).  Notice 
that  Wj2  in  both  (6-23)  and  (6-24)  is  zero.  A  very  simple  and  specific  version 
was  developed  for  application  to  the  Draper  tetrahedral  model  of  LSS.  Modifica¬ 
tion,  generalization,  and  addition  are  underway. 

For  preventing  control  spillover  to  secondary  modes  by  placement  of 

actuators  the  corresponding  formulas  for  determining  the  proper  actuator- 

influence  matrix  B„  are  as  lollows 
F 


rW,,W.7'  =  lip  (6-26) 

k.  I  Z  j  o 

T 

with  being  nonsingular  and  having  the  rank  of  <J>g,  which  is  equal  to  the 
number  of  secondary  modes.  B  is  to  be  a  solution  of 

r 


T 

<P  B_ 
S  F 


=  0 


(6-27) 
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and  is  given  by 


bf  •  K“iS +  ‘>2)b 


(6-28) 


where  Q  =  [Qj,Q9]  and  B  represents  free  parameters.  Let  B  denote  a  desired 
control-influence  matrix  for  the  primary  modes.  Then  B  is  to  be  solved  from 

(~vlvl~lv2  +  V 2 )8  =  8  (6-29) 

where 

[vrv2]  -  <6-30> 

The  one-to-one  mathematical  similarity  between  (6-26)  -  (6-30)  and  (6-21)  - 
(6-25)  is  obvious. 

To  place  actuators  to  prevent  control  spillover  to  all  nonprimgrv  modes, 
if  possible,  is  to  find  a  matrix  B  such  that  the  composite  matrix  Mlp§  is  re- 

alizable  as  actuator  influences,  where  M  deontes  the  mass  (inertia)  matrix  and 
the  orthogonality  of  mode  shapes  is  used.  The  influences  of  a  real  actuator  on 
LSS  are  zero  at  some  modes,  and  nonzero  at  others.  So,  the  realization  of 
the  above  ideal  actuator  influences  by  placement  is  to  find  B  such  that  "con¬ 
trol  spillover"  to  some  modes  is  prevented  while  desired  control  influences  at 
some  other  modes  are  obtained.  The  applicability  and  usefulness  of  the  above- 
mentioned  general-purpose  computer  program  initially  only  for  synthesis  of 
actuator  influences  is  obvious. 
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SECTION  7 


STOCHASTIC  OUTPUT  FEEDBACK  CONTROL 


7.1  Introduction 


The  minimum  variance  fixed  form  output  feedback  compensator  has  been 
recently  proposed  [1,2]  as  a  possible  solution  to  the  problem  of  designing 
reduced  order  controllers  for  plants  of  large  dimension.  The  zero  order 
limit  of  this  compensator,  the  stochastic  output  feedback  controller  [3], 
is  of  particular  interest  because  it  is  the  easiest  such  compensator  to 
implement  and  simultaneously  results  in  the  simplest  design  problem.  The 
price  of  this  simplicity  is  that  this  compensator  generally  achieves  poorer 
performance  than  optimal  compensators  of  higher  order.  As  a  result,  the 
stochastic  output  feedback  controller  should  prove  useful  in  assessing 
complexity/performance  tradeoffs  in  the  design  of  simple  compensators  for 
complex  plants. 

Like  the  reduced  order  compensators  considered  earlier  by  Levine, 

Athans  and  Johnson  [4,5,6],  however,  even  in  the  zero  order  limit  the 
minimum  variance  fixed  form  output  feedback  compensator  gains  and  cost  are 
difficult  to  compute,  especially  for  plants  of  large  dimension.  Consequently, 
upper  and  lower  bounds  on  the  optimal  constant  gain  cost  that  are  relatively 
easy  to  compute  have  been  developed  and  are  described  in  the  following 
sections.  Specifically,  Section  7.2  presents  a  more  detailed  discussion  of 
the  stochastic  output  feedback  controller  and  its  relation  to  several  other 
compensator  formulations  such  as  the  minimum  order  dual  observer  based 
compensator  [1,7]  and  the  Levine-Athans-Johnson  reduced  order  compensator 
[4,5,6].  Section  7.3  derives  an  easily  computed  lower  bound  on  the  optimal 
stochastic  output  feedback  controller  cost  while  Section  7.4  derives  three 
easily  computed  upper  bounds  on  the  cost.  Finally,  Section  7.5  discusses  the 
principal  results  of  the  previous  sections. 

7. 2  Minimum  Various  Fixed  Form  Output  Feedback  Compensator 

Consider  a  plant  of  the  form 

ic  ( t )  =  Ax(t)  +  Bju(t)  +  v(t)  (7-1) 

with  observations 

Z(t)  =  Cx(t)  +  w(t)  (7-2) 

Here  x(t)  £  Rn  is  the  plant  state  vector,  jj(t)  £  Rm  is  the  control  input 
vector,  y(t)  e  is  the  observation  vector  and  A,  B,  C  are  constant  matrices 
of  dimension  nxn,  nxm  and  pxn,  respectively.  The  vectors  v(t)  £  Rn  and 
w(t)  £  R^  are  white  noise  processes  representing  plant  driving  noise  and 


observation  noise,  respectively.  These  are  independent  stationary  zero-mean 
processes  with  covariances 


and 


E{v(t)  v'(T)} 
E{w(t)  w'(T)} 


=  V<S(t  -  x) 
=  W6(t  -  t) 


(7-3) 


As  defined  in  [1,2],  the  minimum  variance  fixed  form  output  feedback  compen¬ 
sator  is  of  the  form 


u(t)  =  Hz(t)  +  Mjr(t)  (7-4) 

where  z^(t)  £  is  the  compensator  state  vector  whose  dynamics  are  described 
by 


£(t)  =  Fz(t)  +  Gjr  (t)  (7-5) 

Here  F,  G,  P  and  M  are  constant  matrices  of  dimension  qxq,  q*p,  mxq  and  mxp, 
respectively.  If  mean-square  stabilizing  compensators  of  this  form  exist 
for  given  dynamic  order  q,  the  optimal  compensator  is  defined  as  the  one 
that  minimizes  the  asymptotic  state  covariance  weighted  by  the  n*n  matrix  Q 
of  design  parameters  with  respect  to  the  matrices  F,  G,  H,  and  M,  i.e., 

J  =  lira  E{x' (t)  Q  x(t)}  (7-6) 

t-x» 


Necessary  conditions  for  the  optimality  of  this  compensator  are  presented 
in  [1,2]  for  compensators  of  arbitrary  dynamic  order  q.  Two  cases  of 
particular  interest  are  the  stochastic  output  feedback  compensator  (SOFC) , 
corresponding  to  q  =  0,  and  the  minimum  order  dual  observer  based  compen¬ 
sator  (MODOBC) ,  corresponding  to  q  =  n  -  m. 

For  the  SOFC  problem,  the  necessary  conditions  reduce  to 

(A  +  BMC) X  +  X(A  +  BMC)'  +  V  +  BMWK ' B '  =  0  (7-7) 

M  =  -(B’PB)_1B’PXC'W_1  (7-8) 

P  (A  +  BMC)  +  (A  +  BMC) '  P  +  Q  =  0  (7-9) 

assuming  the  required  inverse  exists.  The  n*n  matrix  X  represents  the 
asymptotic  state  covariance  matrix  so  the  cost  may  be  expressed  as 


J  =  tr [QX] 


(7-10) 


<4 
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As  In  Che  case  of  the  compensators  considered  In  [4,5,6]  It  is  to  be  expected 
that  the  cost  functional  will  exhibit  local  minimum  [8]  so  th  t  equations 
(7-7),  (7-8)  and  (7-9)  will  not  possess  a  unique  solution. 

Specific  necessary  conditions  can  also  be  developed  for  the  MODOBC 
problem  [1],  Here  the  process  is  more  difficult,  but  the  solution  can  be 
shown  to  be  unique  [7]  implying  the  conditions  are  also  sufficient  for 
optimality.  These  equations  are  derived  by  first  transforming  the  problem 
into  standard  controllable  form  and  then  applying  the  necessary  conditions 
for  the  minimum  variance  fixed  form  output  feedback  compensator.  Determination 
of  the  optimal  MODOBC  cost  then  reduces  to  the  solution  of  a  single  n*n 
Riccati  equation  while  determination  of  the  gains  requires  the  solution  of  an 
additional  (n  -  m)  *  (n  -  m)  Ricatti  equation.  This  compensator  is  of  general 
interest  because  it  is  of  sufficient  dynamic  order  to  allow  arbitrary 
placement  of  the  closed  loop  poles  [7]  and  is  of  specific  interest  here 
because  the  optimal  MODOBC  cost  represents  a  lower  bound  on  the  optimal  SOFC 
cost.  (This  follows  from  the  fact  that  any  compensator  of  order  zero  is  a 
degenerate  case  of  a  compensator  of:  order  n  -  m,  assuming  n  >  m). 

Before  proceeding  with  derivations  of  other  bounds  on  the  optimal  SOFC 
cost,  it  is  useful  to  put  this  reduced  order  compensator  problem  in  perspective 
relative  to  two  other  reduced  order  formulations.  First,  note  that  no 
explicit  control  penalty  is  Included  in  the  cost  functional  considered  here. 

As  noted  in  [1],  no  such  direct  penalty  is  required  because  as  the  compensator 
gains  increase,  so  does  the  contribution  of  the  observation  noise  to  the  state 
covariance.  In  particular,  lack  of  an  explicit  control  penalty  does  not 
result  in  a  singular  problem  in  this  case  as  it  does  in  the  standard  linear 
regulator  problem.  If  an  explicit  control  penalty  is  included  in  the  cost, 
the  minimum  variance  compensator  formulation  of  Martin  and  Bryson  [9,10]  is 
obtained.  Specifically,  the  cost  functional  considered  in  [9,10]  is 

J  =  lim  E{x'(t)  Q.  x(t)  +  u'(t)  Q~  u(t)}  (7-11) 

t-x»  1  c 

providing  the  designer  greater  freedom  to  control  the  magnitude  of  the  compen¬ 
sator  gains  at  the  expense  of  specifying  additional  design  parameters.  The 
principal  disadvantage  of  this  compensator  is  that  it  becomes  ill-posed  in 
the  constant  gain  limit  since  the  cost  becomes  infinite  if  the  control  is 
allowed  to  contain  direct  feedthrough  from  a  plant  output  contaminated  by 
white  noise.  Furthermore,  since  the  SOFC  problem  considered  here  does  not 
penalize  the  control,  the  resulting  optimal  cost  represents  the  best 
performance  achievable  with  any  constant  gain  compensator.  Also,  since  the 
optimal  gains  are  finite,  they  provide  a  physically  useful  indication  of  the 
difficulty  of  achieving  such  performance. 

The  other  reduced  order  compensator  formulation  with  which  the  present 
formulation  should  be  compared  is  the  Levine-Athans- Johnson  formulation 
[4,5,6].  The  compensator  structure  considered  in  this  case  is  the  same,  but 
the  plant  driving  noise  and  observation  noise  are  replaced  by  Initial  state 
uncertainties  and  the  cost  is  the  same  as  that  for  the  linear  regulator 
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problem.  Consequently,  the  resulting  optimal  compensator  depends  on  the 
Initial  state  statistics  which  may  be  more  difficult  to  specify  than 
observation  and  driving  noise  statistics.  Also,  at  least  partially  as  a 
result  of  the  greater  number  of  design  parameters  involved,  the  necessary 
conditions  for  this  compensator  are  more  difficult  to  solve  than  the  SOFC 
necessary  conditions. 

7.3  A  Lower  Bound  on  the  Optimal  SOFC  Cost 

As  noted  in  the  last  section,  the  optimal  MODOBC  cost  is  a  lower 
bound  on  the  optimal  SOFC  cost.  However,  it  is  possible  to  derive  a  more 
easily  computed  lower  bound  by  noting  that  the  necessary  conditions  can  be 
solved  in  the  special  case  where  B  is  a  nonsingular  nxn  matrix.  This 
approach  also  has  the  advantage  of  providing  additional  insight  into  the 
nature  of  the  SOFC  problem  and,  in  fact,  motivates  the  upper  bounds 
developed  in  the  next  section.  Specifically,  if  B-  exists  and  P-^  exists, 
determination  of  the  optimal  SOFC  gain  reduces  to  the  solution  of 


M  =  -B_1X  C’W-1 

o  o 

(7-12) 

+  X  A'  +  V  -  X  C’W-1CX  =  0 

(7-13) 

o  o  o  o 


Note  that  in  this  case,  the  optimal  gain  Mq  is  independent  of  the  cost 

weighting  matrix  Q  provided  only  that  it  is  sufficiently  well-behaved  for  the 
solution  of  equation  (7-9)  to  be  positive  definite.  This  would  appear  to 
imply  that  the  effect  of  the  optimal  gain  M  is  to  minimize  all  of  the 
asymptotic  state  covariances  independently. 

An  even  more  significant  observation  is  that  the  asymptotic  state 
covariance  X  is  independent  of  B  itself  and  only  depends  on  the  existence  of 

B-*.  Because  of  this,  for  any  given  matrix  Q  the  same  optimal  cost  Jq  is 

achieved  by  applying  the  optimal  gain  Mq  (which  does  depend  on  B) ,  around 

any  open  loop  plant  described  by  the  same  A  and  C  matrices  and  subject  to  the 

same  disturbances  V  and  W,  provided  only  that  B  exists.  In  addition,  if 
Jg(M)  represents  the  cost  associated  with  applying  the  gain  M  around  the  open 

loop  system  (A,  B,  C,  V,  W),  then  for  any  B  and  § 

BM  -  BM  ->  Jg(M)  *  Jg(M)  (7-14) 

A. 

provided  A,  C,  V,  W  and  Q  are  the  same  in  both  cases.  Now  suppose  M*  is  the 
optimal  SOFC  gain  for^the  plant  (A,  B,  C,  V,  W)  and  that  B  is  not  necessarily 
invertible,  but  that  B  is.  Thus,  we  can  define  an  "equivalent  gain"  MT  for 
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It  • 


mam 


the  plant  (A,  B,  C,  V,  W)  by 


M 


B  BM 


(7-15) 


so  that 


<y«vT  A  /V  ”v4-  AJ, 

BM  =  BM  =>  JS(M  )  =  Jc(M*) 

D  D 


(7-16) 


However,  since  B  *  exists,  the  optimal  gain  M*  associated  with  (A,  B,  C,  V,  W) 
can  be  computed  from  (7-12)  with  the  optimal  cost  computed  from  (7-13)  as 

Jq  =  tr[QXQ)  (7-17) 

Clearly,  then,  we  must  have 

Jq  <  Jg(M+)  =  Jg(M*)  (7-18) 


so  that  Jq  represents  a  lower  bound  on  the  optimal  SOFC  cost  for  the  plant 
(A,  B,  C,  V,  W)  for  any  actuator  influence  matrix  B. 

Note  that  if  B  1  does  not  exist,  the  computation  of  the  lower  bound 
cost  Jq  is  easier  than  the  computation  of  the  MODOBC  cost  since  it  is  not 

necessary  to  first  transform  the  plant  into  standard  controllable  form.  If 

B  *  does  exist,  it  follows  that  m  =  n  so  the  MODOBC  problem  and  the  SOFC 
problem  are  identical.  In  this  case,  the  problem  Is  effectively  in  standard 

controllable  form  (i.e.,  the  required  transformation  is  z^(t)  =  B  *x(0>  so 
no  additional  work  is  necessary.  Finally,  note  that  the  covariance  equation 
(7-13)  is  the  same  algebraic  Riccati  equation  that  is  used  to  compute  the 
Kalman  filter  gains  required  to  estimate  the  state  of  the  open  loop  plant. 
Consequently,  many  theoretical  and  numerical  results  are  available  pertaining 
to  its  solution.  For  example,  it  is  known  [11]  that  if  the  open  loop  plant 
is  stabilizable  and  detectable  with  respect  to  the  noise  sources,  there 
exists  a  unique  positive  semidef inite  solution  Xq  that  is  positive  definite 

if  and  only  if  the  system  is  completely  controllable  from  the  plant  driving 
noise. 


7.4  Suboptimal  SOFC  Compensators 

The  simplifications  of  the  optimal  SOFC  problem  resulting  when  B  * 
exists  suggest  that  it  may  be  profitable  to  cast  the  problem  in  terms  of  a 

generalized  inverse  [12,13]  of  B  when  B  *  does  not  exist.  To  do  this,  note 
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that  the  gain  expression  may  be  written  as 


M  =  -B+XC'W-1 


(7-19) 


where  the  generalized  inverse  B  is  given  by 


B+  =  (B'PB)_1B'P 


(7-20) 


and  is  in  fact  a  weighted  least  squares  inverse  [13]  of  B.  That  is,  defining 
S  as  the  positive  definite  square  root  of  P  (assuming  P  is  positive  definite),, 
i.e. , 


S'S  =  P  ,  S  >  0 


(7-21) 


the  vector  x^  =  is  the  minimum  S-norm  approximate  solution  of  the 

inconsistent  equation  Bx  =  Consequently,  an  alternate  way  of  viewing  the 
optimal  SOFC  design  problem  is  as  that  of  minimizing  the  error 


e  -  | |Bu  +  XC'W-1^! |s  -  | |SBu  +  SXC'W"1^!  j  (7-22) 

subject  to  the  constraints 

u  =  (7-23) 

(A  +  BMC)X  +  X(A  +  BMC)'  +  V  +  BMWM'B'  =  0  (7-24) 

P(A  +  BMC)  +  (A  +  BMC) 'P  +  Q  =  0  (7-25) 


where  S'S  =  P.  Physically,  equation  (7-23)  just  specifies  the  form  of  the 
control  law,  equation  (7-24)  defines  X  as  the  closed  loop  state  covariance 
and  equation  (7-25)  relates  the  norm  ||*|L  to  the  cost  weighting  matrix  Q. 

Note  also  that  if  B  exists,  it  is  imnediately  apparent  that  the  minimum 
error  e  =  0  is  achievable  independent  of  the  norm  ||*||g-  This  formulation 

of  the  problem  thus  makes  immediately  clear  the  independence  of  the  optimal 
gains  from  the  design  parameters  noted  in  Section  7.3  when  B  is  non¬ 
singular. 

This  formulation  of  the  optimal  SOFC  design  problem  also  suggests 
several  suboptimal  control  problems  that  may  be  useful  in  obtaining  either 
upper  bounds  on  the  optimal  SOFC  cost  or  initial  gain  estimates  for 
iteratively  computing  the  optimal  SOFC  gains.  For  example,  since  the 
principal  difficulty  in  computing  the  optimal  gains  is  the  coupling  between 


equations  (7-24)  and  (7-25),  one  possible  approach  would  be  to  assume  a 
least  squares  weighting  S  and  minimize  the  error  defined  in  equation  (7-22) 
subject  to  constraints  (7-23)  and  (7-24).  Defining  P  =  S'S,  the  solution 
to  this  problem  is 

Ml  =  -(B'PB)-1B'PX1C,W-1  (7-26) 

where  the  closed  loop  covariance  satisfies  the  equation 

AXX  +  XXA'  +  V  -  X1C'W-1CX1  =  -(I  -  B(B'PB)-1B,P)X1C'W-1CX1(I  -  PB(B,PB)"1B' ) 

(7-27) 

A  particularly  convenient  choice  of  weighting  matrix  S  would  be  the  n*n 
identity  matrix,  reducing  the  problem  to  a  standard  least  squares  minimi¬ 
zation. 

Another  possible  approach  to  decoupling  equations  (7-24)  and  (7-25) 
would  be  to  assume  a  value  X  for  the  closed  loop  state  covariance  and  minimize 
the  error  defined  in  equation  (7-22)  subject  to  constraints  (7-23)  and 
(7-25).  The  solution  to  this  problem  is 

M2  =  - (B ' P2B) _  XB ' P2XC  *  W~ 1  (7-28) 

where  P2  satisfies  the  equation 

p2a  +  a'p2  +  q  -  p2b(b,p2b)_1b'p2xc,w_1c  -  c,w_1cxp2b(b,p2b)_1b,p2  *  0 

(7-29) 

with  the  least  squares  weighting  matrix  S2  defined  as  the  positive  definite 
square  root  of  Pj.  A  particularly  convenient  choice  of  covariance  matrix 
would  be  the  lower  bound  covariance  Xq  defined  by  equation  (7-13)  since  the 
suboptlmal  solution  then  becomes  exact  when  B  is  nonsingular. 

Finally,  a  third  suboptimal  SOFC  problem  results  if  both  the  least 
squares  weighting  matrix  S  and  the  closed  loop  covariance  matrix  X  are  fixed. 
In  particular,  assuming  S  -  S  and  X  ■  X,  the  suboptimal  gain  is  given  by 

M3  -  -(B'PB)_1B,PXC,W~1  (7-30) 
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where  P  *  S'S  as  before.  Note  that  if  S  is  taken  as  the  nxn  identity  matrix 
this  approach  ia  very  similar  to  Kosut's  minimum  norm  suboptimal  compensator 
for  the  Levine-Athans-Johnson  reduced  order  compensator  [14]. 

7.5  Conclusions 


Motivated  primarily  by  the  exact  solution  of  the  optimal  SOFC  problem 
when  the  actuator  influence  matrix  B  is  nonsingular,  bounds  have  been  derived 
for  the  optimal  SOFC  cost  for  plants  with  singular  B  matrices.  These  bounds 
satisfy  the  inequality 

min  (J  ,J  )  *  max  (J  ,J  )  <  J  <  T  <  J  <  J, , .  (7-31) 

O  —  O  ~  T  TT  fir 

where  J_  is  the  optimal  MODOBC  cost  for  the  plant  (A,  B,  C,  V,  W) ,  Jq  is  the 

optimal  SOFC  cost  for  the  same  plant  but  with  a  nonsingular  B  matrix,  J  is 
the  optimal  SOFC  cost  and  J  ,  J  ,  J , , ,  are  the  costs  associated  with  the 

i  TT  TTT 

suboptimal  gains  M^,  M2,  and  M^  ranked  in  order  of  increasing  cost.  It 

should  be  noted  that  if  B  is  nonsingular,  it  follows  immediately  that 

J  =  J  =  J  =  J,  and  that  if  X  is  taken  as  the  solution  X  of  equation  (7-13) 
o—l  o 

in  computing  M2  and  M^,  ^  and  are  also  equal  to  the  optimal  cost. 

In  terms  of  computational  difficulty,  the  cost  Jq  is  clearly  the 

easiest  to  evaluate,  requiring  the  solution  of  a  single  Riccati  equation.  If 
the  plant  is  in  standard  controllable  form,  the  optimal  MODOBC  cost  J_  is 
equally  easy  to  compute,  but  if  it  is  not  (as  is  normally  the  case),  the 
evaluation  of  J_  can  be  considerably  more  difficult.  Of  the  upper  bounds, 
the  cost  J  is  the  easiest  to  evaluate,  requiring  the  solution  of  a  single 

J  A 

Riccati  equation  if  X  is  available.  Evaluation  of  is  somewhat  more 

difficult  since  it  requires  the  solution  of  equation  (7-27)  which  is  not  a 
standard  Riccati  equation.  Similarly,  the  cost  should  be  still  more 

difficult  to  evaluate  because  of  the  more  complex  nonlinearity  of 
equation  (7-29).  Note,  however,  that  all  of  these  bounds  on  the  optimal 
SOFC  cost  are  more  easily  evaluated,  assuming  they  exist,  than  the  optimal 
cost  itself,  since  they  involve  the  solution  of  uncoupled  equations. 

In  terms  of  applications,  three  possibilities  are  apparent.  First, 
the  lower  bound  Jq  and  the  upper  bound  should  be  useful  in  preliminary 

evaluation  of  optimal  SOFC  performance  for  particular  applications.  This 
would  be  particularly  convenient  if  the  lower  bound  covariance  Xq  were  used 

in  the  computation  of  the  suboptimal  cost  J^.  Secondly,  an  examination  of 

the  suboptimal  compensators  M^  and  M2  suggests  a  possible  computational 

algorithm  for  evaluating  the  optimal  SOFC  gains  M.  Specifically,  if 
reasonably  efficient  algorithms  could  be  developed  for  solving 
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equations  (7-27)  and  (7—29) , Ait  might  be  reasonable  to  approach  the  problem 
by  first  fixing  either  P  or  X,  solving  the  appropriate  equation  to^obtain 
X  or  P,  respectively,  and  using  these  to  derive  updated  values  of  P  or  X 
from  the  other  equation.  Finally,  any  of  the  suboptimal  gains  M^,  M2  or  M^ 

could  be  used  as  a  starting  estimate  for  any  iterative  calculation  of  the 
optimal  SOFC  gains,  as  has  been  done  with  Kosut's  suboptimal  compensator  in 
evaluating  the  Levine-Athans- Johnson  reduced  order  compensator  gains  [15] . 
Because  of  its  relative  ease  of  evaluation,  M^  appears  to  be  a  particularly 
attractive  candidate  for  this  application. 
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SECTION  8 


CALCULATION  OF  OPTICAL  LINE-OF-SIGHT 


8.1  Introduction 


For  the  ACOSS  test  problem  It 2  we  need  a  line-of-sight  model  giving 
the  law  of  displacement  of  the  image  of  a  target  when  the  structure  deforms. 
The  purpose  of  the  model  is  to  give  line-of-sight  changes  of  the  same  order 
of  magnitude  as  those  produced  by  current  DARPA  large  wide-angle  optical 
systems. 

We  assume  that  the  mirror  surfaces  are  simply  displaced  and  maintain 
their  nominal  shape.  We  also  neglect  any  effect  of  light  redistribution  in 
the  image. 

The  structural  model  given  is  shown  in  Figures  8-1  and  8-2.  It  was 
assumed  that  the  base  optical  design  would  be  close  to  an  existing  DARPA 
design.  The  mirror  shapes  are  off  axis  sections  of  rotationally  symmetric 
surfaces  which  are  coaxial.  As  we  want  only  a  first  order  model  and  neglect 
the  influence  of  light  distribution  in  the  image  of  a  point,  the  asphericity 
of  each  mirror  was  disregarded. 

8.2  Results 


The  line-of-sight  change  is  the  sum  of  the  changes  due  to  the 
displacement  of  each  of  the  components  of  the  optical  system  taken  indi¬ 
vidually.  For  a  given  mirror  or  focal  plane  displacement,  the  location  of 
the  intermediate  and/or  final  images  are  calculated  using  the  well-known 
(see  Ref.  1)  first  order  optics  equations 


1 


(8-1) 


m 


(8-2) 


relating  the  distance  d^  between  the  image  and  the  mirror,  the  distance 

between  the  object  and  the  mirror  (both  distances  being  taken  along  the 
optical  axis),  the  radius  of  curvature  r  of  the  mirror  and  the 
magnification  m.  The  equations  (8-1)  and  (8-2)  utilize  the  vertex  of  the 
mirror  as  origin;  therefore,  when  a  mirror  is  displaced,  the  results  are 
given  in  relation  to  the  displaced  positions,  and  this  is  taken  into 
consideration  in  the  transfer  to  the  next  mirror. 
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Using  the  paraxial  optics  equations  (8-1)  and  (8-2),  expressions 
relating  optical  element  movements  and  the  resultant  longitudinal  defocus 
and  transverse  defocus  (x„,  y_)  are  given  for  the  three-mirror  ACOSS 

r  t 

system.  Translations  of  the  mirrors  are  represented  by  the  displacements 
of  the  vertices,  and  rotations  are  defined  for  each  mirror  with  the  pivot 
point  located  at  its  own  vertex  (on  the  optical  axis).  Movements  of  the 
focal  surface  are  also  considered;  x^.  and  y^  are  the  coordinates  of  the 

image  point,  and  the  pivot  point  for  rotation  is  again  at  the  intersection 

of  the  focal  plane  with  the  optical  axis  (see  Figure  8-3).  Line-of-sight  (LOS) 

error  is  determined  by  the  transverse  defocus  (x  ,  y  ). 

F  F 

The  pertinent  variables  are  defined  as  follows: 


Rp  =  radius  of  curvature  of  the  primary 

R^_  =  radius  of  curvature  of  the  tertiary 

t^  =  axial  distance  between  the  primary  and  the  secondary 


t 2  =  axial  distance  between  the  secondary  and  the  tertiary 

We  will  use  the  letter  a  as  a  running  index  which  stands  for  the  letter  p 
for  primary,  s  for  secondary,  t  for  tertiary  and  F  for  focal  plane,  x. 
will  represent  the  displacement  of  the  element  a  along  the  ith  axis  1 
(see  Figure  8-3)  and  0  will  represent  the  counter-clockwise  rotation 

of  the  element  a  about  the  ith  axis  (looking  towards  the  origin) .  With 
these  notations  the  defocus  components  (x^,,  y^,,  z^)  will  be  given  by 


■x.  +  x  -  R 

lp  It  p 


+  Xlt  -  X1F  +  93Fyl 


62p  +  2e2s(/+t.)  -  262tT]2T^ 


[-x0  +  xt 

2p  t 


+  R  0.  -  20. 

p  lp  Is 


{?  +  0  +  29itT]  Trht; 


(8-3) 


+  X2t  "  X2F  "  93FXI 


(8-4) 


(X3P  "  2x3s  +  X3t} 


“  yl6lF  +  XI02F 


+  x„.  -  x. 


[R  -  R  +  2 (t.  +  t,)] 
p  t  12 


2  3t  3F 
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(8-5) 


where 


T 


R 

_2. 

2 


+  t, 


+  t 


2 


The  LOS  error  is  given  by 


where  f  is  the  focal  length  of  the  system  and  is  expressed  as 


(8-6) 


(8-7) 


f 


R  R 
1  E 


2(2T 


V 


The  magnitude  of  the  LOS  error  is  given  by 


(8-8) 


(8-9) 


The  location  of  the  image  plane  is  a  function  of  the  parameters  of  the 
three-mirror  optical  system.  The  distance  BF  between  the  tertiary  mirror 
and  the  image  plane  (see  Figure  8-3)  is  given  by 


BF 


RfcT 


2T  -  R 


(8-10) 


and  for  all  cases  the  focal  surface  must  be  placed  at  this  location,  or, 
alternately  the  radius  of  curvature  of  the  tertiary  should  be  calculated 
using 


2T  »  BF 
Rt  '  T  +  BF 


(8-11) 


so  that  the  final  image  coincides  with  the  focal  plane. 

In  the  above  expressions,  all  radii  of  curvature  and  axial  distances 
should  have  positive  values. 

It  should  also  be  reiterated  that  all  rotations  of  the  individual 
elements  are  defined  with  the  pivot  points  on  the  optical  axis.  The  justi¬ 
fication  for  this  definition  lies  in  the  fact  that  all  three  elements  of 
the  optical  system  are  coaxial,  i.e.,  the  centers  of  curvature  of  these 
elements  lie  along  a  common  axis. 
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As  a  final  note,  the  (x  ,  y  ,  z  )  expressions  were  derived  with  first 

r  r  r 

order  paraxial  optics.  As  a  consequence  of  this  approximation,  the  angle  of 
incidence  of  the  object  radiation  and  all  displacements  and  rotations  of 
the  individual  elements  must  be  small.  If  this  condition  is  not  met,  not 
only  will  higher  order  terms  become  appreciable,  but  the  focal  length  and 
consequently  the  backfocus  BF  will  change  and  aberrations  of  the  optical 
system  will  no  longer  be  negligible.  Therefore  those  equations  cannot  be 
realistically  used  in  an  open  attitude  control  loop. 

8. 3  Application  to  the  Test  Case  #2 

For  the  test  case  til  corresponding  to  the  structural  model  given  on 
Figures  8-1  and  8-2,  the  primary  and  tertiary  mirror  vertices  are  coincident 
and  the  focal  plane  and  tertiary  mirror  (flat)  are  coplanar. 

The  numerical  values  of  the  basic  variables  which  have  been  used  are: 

R 

P 


BF  = 

T,  Rt  and  f  are  calculated  using  Equations  (8-6),  (8-11),  and  (8-8). 
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SECTION  9 


DRAPER  MODEL  #2 


9. 1  Introduction 


In  order  to  assess  the  performance,  sensitivity,  and  hardware 
requirements  of  the  various  active  structural  control  methods  which  are 
being  developed,  a  universal  system  model  is  required.  This  section  contains 
a  complete  description  of  Draper  Model  #2  which  is  a  simple  but  realistic 
evaluation  model.  The  design  of  this  system  was  driven  by  the  desire  to 
incorporate  certain  attributes  into  the  overall  system  characteristics. 

The  desired  features  were 

•  structural  design  based  on  realistic  sizes  and  weights 

•  a  simple  unclassified  optical  system  with  associated 
performance  measures  and  tolerances 

•  a  set  of  disturbances  typical  of  equipment  vibration  and 
attitude  control  (slew). 

The  resulting  model  is  described  in  detail  in  Section  9.2  along  with  the 
resulting  mode  shapes  and  natural  frequencies  of  the  system.  Section  9. A 
contains  a  description  of  the  line  of  sight  (LOS)  performance  measure 
including  theory  and  implementation.  The  tolerances  on  the  optical  system 
and  a  set  of  disturbances  are  given  in  Section  9.5.  A  perturbed  model 
which  can  be  used  to  evaluate  the  sensitivity  of  the  control  system  is 
given  in  Section  9.6. 

9. 2  Structural  Design 

The  example  problem  is  shown  in  Figure  9-1.  It  consists  of  a 
flexible  optical  support  structure  and  an  isolated  equipment  section  which 
contains  the  solar  panels. 

9.2.1  Optical  Support  Structure 

The  optical  support  structure  consists  of  the  upper  mirror  support 
truss,  lower  mirror  support  truss,  and  the  metering  truss  which  maintains 
mirror  separation.  The  finite  element  model  is  shown  in  Figure  9-2.  It 
contains  35  node  points  and  117  beam  elements.  The  structure  is  configured 
as  a  truss,  but  it  is  assumed  that  all  joints  allow  a  full  moment  connection. 
Thus,  both  bending  and  axial  stiffness  are  included  in  the  model  for  all 
structural  members.  All  members  are  assumed  to  be  graphite  epoxy  hollow 
tubes.  The  sizes  and  section  properties  are  given  in  Table  9-1.  The  mass 
of  this  subsystem  including  the  structural  mass,  mirrors,  and  equipment  has 
been  lumped  at  18  node  points  as  shown  in  Figure  9-3.  The  total  mass  of 
the  optical  support  structure  is  508A  kg. 
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model  No.  2— Finite  element  model. 


Table 


099  x  10~  m  J  =  3.749 


EQUIPMENT  SECTION  &  SOLAR  PANELS 


EOUIPMENT  SECTION  -  MEMBER  CONNECTIVITIES 


LOWER  SUPPORT  TRUSS  -  NODE  LOCATIONS 


DIAGONAL  MEMBERS  OMITTED  FOR  CLARITY 


! 


i 
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METERING  TRUSS  -  NODE  LOCATIONS 


UPPER  SUPPORT  TRUSS  -  NODE  LOCATIONS 


DIAGONAL  MEMBERS  OMITTED  FOR  CLARITY 


UPPER  SUPPORT  TRUSS  -  MEMBER  CONNECTIVITIES 


9.2.2  Equipment  Section 


The  equipment  section  consists  of  a  central  rigid  body  and  two 
symmetric  solar  panels.  The  support  structure  for  the  solar  panels  is  a 
hollow  graphite  epoxy  tube.  The  structural  properties  of  the  solar  panels 
are  not  included  in  the  model.  The  inertia  properties  are  lumped  at  6  node 
points  on  the  support  structure.  Details  of  the  finite  element  model  and 
the  inertia  properties  of  the  equipment  section  are  given  below  in 
Figure  9-4. 


Rigid  Body  Properties 

Solar  Panel  Properties 

Node  44 

Nodes  48,  52,  53,  57 

Node 

50,  55 

m 

=  3500  kg 

m  =  90  kg 

m  = 

180  kg 

I 

XX 

=  2100  kg-m2 

I  =270  kg-m2 

XX  & 

I 

XX 

540  kg-m2 

I 

yy 

55  2100  kg-m2 

i 

zz 

=  4200  kg-m2 
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9.2.3  Isolation  System 


The  optical  support  structure  and  the  equipment  section  are  connected 
by  a  passive  isolator  at  three  points  acting  in  all  three  translational 
directions.  The  isolator  is  designed  to  reduce  the  amount  of  force  trans¬ 
mitted  from  the  equipment  section  to  the  optics.  This  is  accomplished  by 
designing  an  isolator  with  the  natural  frequency  much  lower  than  expected 
disturbance  frequency.  In  addition,  the  isolator  frequency  should  be  lower 
than  the  optical  structure  bending  frequency.  A  passive  isolator  consists 
of  a  spring  and  a  viscous  dashpot  in  parallel  as  shown  in  Figure  9-5. 


PRECISION  SECTION 


EQUIPMENT  SECTION 


_1  1*1*1 
211  V  M1M2 

Figure  9-5. 

If  the  isolator  frequency  is  lower  than  the  lowest  natural  frequencies  of 
the  two  sections,  they  can  be  modelled  by  their  inertia  properties  only. 

For  this  system,  the  isolator  was  sized  to  have  a  frequency  of  0.5  Hz.  The 
dashpot  size  has  been  specified  as  it  is  assumed  to  be  part  of  the  overall 
control  design.  The  spring  constants  are  given  in  Table  9-2. 

9. 3  Structural  Analysis 

A  modal  analysis  was  performed  using  the  NASTRAN  finite  element 
program.  The  input  deck  is  listed  in  Table  9-3.  This  model  contains 
84  dynamic  degrees  of  freedom,  thus  84  frequencies  and  mode  shapes  were 
extracted.  Table  9-4  gives  a  list  of  frequencies  and  descriptions  of  the 
first  fifty  mode  shapes. 

9.4  System  Performance 

The  performance  of  this  system  is  measured  by  the  ability  of  maintain 
line-of-sight  (LOS)  rotation  (jitter)  and  defocus  within  specified  tolerances. 
These  quantities  can  be  computed  using  simple  linear  equations  with  the 
movement  of  the  mirrors  as  inputs.  The  motion  of  the  mirrors  is  a  function 
of  the  movement  of  the  support  points.  This  section  presents  the  equations 
necessary  to  compute  LOS  error  from  support  point  displacements. 
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isol 


ISOLATOR  SPRING  CONNECTIVITIES 
TABLE  9-2 


CBftfl  *5  190 
CB»f»  *6  300 
CBftR  47  700 
CM»0  *4  300 
CBAil  55  700 
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9.4.1  LOS  Error  Algorithm 

The  equations  (9-2),  (9-3),  and  (9-4),  relating  optical  surface 
motion  and  the  resultant  LOS  rotations  about  the  X  and  Y  axes 
below  with  notations  and  orientation  consistent  with  the  structural 
analysis.  The  movement  of  each  mirror  is  defined  to  be  the  displacements 
and  rotations  in  the  global  X,  Y,  Z  directions  at  a  point  on  the  sphere 
(primary,  tertiary)  or  plane  (focal  plane)  which  intersects  the  optical 
(global  Z)  axis.  The  LOS  rotation  and  defocus  are  given  by 


LOSX  =  Y/F 

L0SY  =  X/F 

Defocus  =  Z 

(9-1) 

where 

F 

_ 

8.051 

X 

= 

A. [-X  +  X 

1  p  t 

-  R*0Y  +  A  *0Y 
p  2  s 

-  2T*0Y  ]  +  X  -  Xf 

(9-2) 

Y 

= 

A.  f-Y  +  Y 

1  p  t 

+  R*0X  -  A  *0X 
p  2  s 

+  2T*0X  ]  +  Y  -  Y, 
t  t  f 

(9-3) 

Z 

= 

A.[Z  -  2Z 

3  p  s 

+  zt]  +  zt  -  Zf 

(9-4) 

where 

A1 

= 

0.2987 

A2 

= 

93.90 

A3 

= 

0.0892 

R 

= 

53.9 

T 

— 

66.95 

The  terms  X, ,  Y.,  Z.,  9X.,  0Y.,  0Z.  refer  to  the  translations  and  rotations 
1  l  i  l  i  i 

in  the  global  X,  Y,  and  Z  directions  of  the  primary  (p)  and  secondary  (s) , 
tertiary  (t) ,  and  focal  plane  (f). 

9.4.2  LOS  Error  Algorithm  Implementation 

The  motion  of  each  mirror  is  based  on  the  displacements  of  the 
support  points  and  its  position  relative  to  the  optical  axis.  It  is  assumed 
that  the  mirrors  are  kinematically  mounted  in  order  to  decouple  their  dynamic 
motion  from  the  structure.  A  typical  kinematic  mount  is  shown  in  Figure  9-6. 
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p 

Y 


Figure  9-6.  Kinematic  Mount 


Point  A  is  supported  in  X,  Y,  and  Z  directions.  Point  B  in  the  Y  and  Z 
directions,  and  Point  C  in  the  Z  direction.  Using  small  displacement  theory, 
the  displacement  and  rotation  of  any  point  P  can  be  found. 


where 


> 

ii 

Axa  +  zpA6y  - 

yPA6z 

(9-5) 

Ayp  = 

AyA  -  ZpAex  - 

V0Z 

(9-6) 

Azp  = 

AzA  +  yPA9x  " 

xPA0y 

(9-7) 

A0 

X 

=  <Axc 

-  AzA  +  xc(Aza 

-  AzB)/xB)/yC 

(9-8) 

A0 

y 

=  <Aza 

-  Azb)/xb 

(9-9) 

A0 

z 

=  <AyB 

-  AyA)/xB 

(9-10) 

The  locations  of  the  support  points  and  point  P  for  each  surface  in  this 
system  are  shown  in  Figure  9-7  and  are  summarized  in  Table  9-5. 
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Surface 

Node  A 

Node  B 

Node  C 

X 

P 

Y 

p 

Z 

P 

Primary 

34 

35 

30-28 

10.0 

0.0 

Secondary 

40 

- 

- 

- 

- 

Tertiary 

27 

29 

33-32 

-3.0 

0.0 

F.P . 

11 

9 

40 

-4.0 

-5.0 

0.0 

Table  9-5 


It  should  be  noted  that  since  node  C  of  the  primary  and  tertiary  is  not 
located  at  a  structural  node  point,  the  displacement  can  be  found  by  inter¬ 
polating  between  the  two  nodes  given  in  the  table.  Also,  all  six  degrees  of 
freedom  of  the  secondary  correspond  to  the  six  degrees  of  freedom  at  node  40. 
The  location  of  point  P  is  in  the  local  coordinate  system  of  each  mirror 
which  are  all  parallel  to  the  global  coordinate  system.  Using  equations 
(9-5)  -  (9-10)  and  the  information  in  Table  5-5,  a  displacement  transformation 
matrix  (DTM)  which  relates  the  motion  of  each  mirror  to  the  displacements  of 
the  support  points  in  the  global  coordinate  system. 


If  modal  analysis  is  employed  to  compute  the  displacements,  the  LOS 
error  and  defocus  for  each  modal  vector  can  be  computed  beforehand  and  the 
total  errors  found  by  modal  superposition. 

t  LOSX  | 

|  LOSY  I  =  [$L0S]  {q}  (9-11) 

where  [<J>LQg]  is  a  matrix  which  relates  LOS  error  and  defocus  to  modal 

amplitude.  It  is  computed  by  first  computing  the  mirror  motions  due  to  each 
mode  shape  using  the  DTM  and  then  applying  equations  (9-1)  through  (9-3). 

T 

For  Draper  model  if 2  [<t>LQgJ  is  given  in  Table  9-6  for  modes  1-50.  It  was 
calculated  using  mode  shapes  which  were  normalized  to  unit  generalized  mass. 

9.5  Parameter  Variations 


In  order  to  assess  the  sensitivity  of  the  control  system  to  changes  in 
the  natural  frequencies  and  mode  shapes,  a  set  of  perturbed  frequencies  and 
mode  shapes  can  be  used.  This  perturbed  model  can  be  obtained  by  either 
varying  the  frequencies  and  mode  shapes  directly  or  by  changing  selected 
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Table  9-6 


MOOE 

LOSX 

LOSY 

DE FOCUS 

LOSMAG 

FREQ 

**** 

««»# 

**** 

******* 

****** 

«*** 

1 

0.27050-04 

0.37730-03 

0.93130-09 

0.37670-03 

0.0001 

2 

-.15300-04 

0.84260-06 

-.13630-07 

0.15820-04 

0.0001 

3 

0.49270-03 

-.31330-04 

-.96040-09 

0.49370-03 

0.0001 

4 

0.11660-03 

0.1155D-03 

0.43660-10 

0.16410-03 

0.0001 

5 

-.35020-04 

-.81100-03 

2721D-08 

0.6117D-03 

0.0001 

6 

-.73350-03 

0.46760-04 

0.69650-09 

0.78490-03 

0.0001 

7 

-.26410-06 

0.32630-03 

-.4904D-07 

0.32630-03 

0.1455 

S 

0.34390-06 

-.21410-03 

-.  1099D-05 

0.21410-03 

0.2632 

9 

-.15050-05 

-.34520-06 

-.5238D-05 

0 . 15440-05 

0.3173 

10 

-.22930-03 

0. 44070-06 

-. 1467D-06 

0. 2293D-03 

0.3329 

11 

0.60410-05 

0.10660-05 

0.15590-04 

0.6135D-05 

0.4432 

12 

0.73590-03 

-.67630-05 

0.19060-05 

0 . 7360D-03 

0.5779 

13 

0.23300-05 

0. 87310-03 

-.42760-06 

0.87310-03 

0.5814 

14 

-.10270-05 

-.23620-03 

-. 4831D-07 

0.23620-03 

1.2238 

15 

0.37590-05 

-.15190-03 

-. 7763D-05 

0. 1519D-03 

1.3002 

16 

-.15390-03 

-.14130-05 

0.93640-05 

0.15890-03 

1.3475 

17 

0.27330-04 

-.27770-02 

-.4955D-04 

0.2777D-02 

1.7209 

1C 

0 . 2550D-05 

-.49760-06 

0 . 2634D-05 

0.26370-05 

1.8187 

19 

0. 11510-04 

-.23640-05 

-.54930-06 

0. 11750-04 

1.8187 

20 

-.35970-13 

-.35160-14 

0 . 3030D-14 

0 . 3615D-13 

1.8892 

21 

0.22370-04 

0.13760-01 

0. 90280-05 

0.1376D-01 

2.3635 

22 

0.13450-03 

0.27240-05 

-.15090-04 

0 . 1346D-03 

2.9395 

23 

-.43550-05 

-.26020-03 

0.25000-05 

0.26020-03 

3.1759 

24 

0 . 5573D-06 

-.45360-03 

- .42560-06 

0. 45360-03 

3.3873 

25 

-.18140-13 

-.43980-15 

0.26010-14 

0.18150-13 

5.1617 

26 

-.33950-04 

-.70130-06 

0.47470-05 

0.33960-04 

5.2603 

27 

0.10330-14 

0.72640-16 

- . 9614D-16 

0.10400-14 

7.8769 

23 

- .  4275D-03 

-.17010-02 

-.19500-03 

0. 1753D-02 

8.1168 

29 

0.29160-03 

0. 47750-06 

-.3636D-05 

0 . 2916D-03 

8.3600 

30 

0. 13660-01 

-.11030-04 

0.1072D-02 

0.13680-01 

8.5706 

31  • 

-.33220-04 

- . 2915D-06 

- . 3755D-05 

0.38220-04 

8.8135 

32 

-.44940-04 

0 . 7780D-07 

0 . 9546D-06 

0.44940-04 

8.0135 

33 

0.37620-05 

0.69950-05 

-.52400-06 

0.9750D-05 

11.3462 

34 

-.12250-02 

0.1070D-02 

0. 39820-03 

0.2236D-02 

11.4978 

35 

0 . 10Q30-Q1 

-.56620-03 

- .  1593D-01 

0.10990-01 

12.7258 

36 

0.20970-04 

0. 74750-03 

-  .6C65D-04 

0.7478D-03 

13.5832 

37 

-.71430-03 

-.10110-01 

0.11420-02 

0.10140-01 

13.7141 

33 

0.39530-02 

0.59350-02 

-. 1084D-02 

0.71310-02 

14.1604 

39 

-.14140-01 

0.32130-02 

-.49660-02 

0 . 1451D-01 

15.6523 

40 

0.19070-02 

0.10440-01 

- . C294D-03 

0.10620-01 

16.0724 

41 

-.67600-02 

-.22200-02 

0. 1413D-01 

0 . 7134D-02 

16.5248 

42 

0.33530-02 

-.95400-03 

-.13750-01 

0.39690-02 

16.7453 

43 

0.11320-03 

-.64050-02 

-. 1739D-02 

0.8406D-02 

17.1555 

44 

0.3244D-02 

-.1166D-02 

0.71740-02 

0. 34480-02 

17.8203 

45 

-.64390-03 

-.14460-01 

- . 1545D-02 

0 . 1447D-01 

19.0713 

46 

- .  2266D-01 

0.16690-01 

0. 52170-01 

0. 26140-01 

23.7716 

47 

0.90530-02 

0.37110-01 

-.22930-01 

0.33200-01 

24.4140 

43 

0 . 266S0-03 

0.12050-01 

- . 3317D-Q2 

0 . 1206D-01 

25.908  > 

49 

-.14440-01 

-.10420-01 

0.23620-02 

0.17810-01 

26.3625 

50 

0.94160-14 

-.62570-14 

-. 1485D-14 

0.11310-13 

26.4292 
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physical  parameters  which  will  cause  a  variation  in  the  modal  characteristics. 
Since  the  effects  of  modelling  and  manufacturing  errors  are  of  primary 
concern,  the  latter  method  was  chosen.  Two  sets  of  physical  parameter 
variations  are  presented  here.  Modal  analyses  including  the  modal  LOS 
error  were  conducted  using  NASTRAN.  The  resulting  natural  frequencies  and 
mode  shapes  are  given  and  compared  with  the  original  system. 

Two  sets  of  parameter  variations  are  summarized  in  Table  9r7.  The 
changes  are  confined  to  the  optical  support  structure  and  involve  mass 
distribution  changes  on  the  upper  optical  support  truss  and  stiffness 
changes  in  the  metering  truss.  These  parameters  were  chosen  based  on  their 
impact  on  LOS  error  in  several  critical  modes.  The  amount  of  variation  is 
driven  by  the  desire  to  vary  the  natural  frequency  by  10%. 

Dynamic  analyses  were  conducted  on  both  perturbed  models  using  NASTRAN. 
The  resulting  natural  frequencies  are  given  in  Table  9-8  along  with  those 
of  the  original  model.  The  modal  LOS  errors  about  the  X  and  Y  axes  and 
the  overall  magnitude  are  given  for  all  three  cases  in  Table  9-9.  The 
results  show  that  both  sets  of  parameter  variations  had  more  impact  on  the 
higher  frequency  modes  because  the  upper  support  truss  and  the  metering  truss 
contain  a  higher  percentage  of  the  total  energy  in  these  modes  than  in  the 
lower  modes. 

The  perturbed  models  presented  here  are,  of  course,  two  point  designs 
from  an  infinite  number  of  possible  parameter  variations.  In  order  to  fully 
evaluate  the  stability  limits  of  a  particular  control  system  design,  it  may 
be  necessary  to  evaluate  several  perturbed  models.  Another  type  of  variation 
which  must  also  be  considered  are  the  changes  in  amplitude  and  frequency 
of  the  system  disturbances. 

9.6  System  Disturbances  and  Tolerances 

9.6.1  Description 

The  disturbance  model  consists  of  two  sinusoidal  forces  which  are 
applied  to  a  point  on  the  optical  support  structure  and  a  point  on  the 
equipment  section.  The  location,  direction  magnitude  and  frequency  are 
given  in  Figure  9-8.  These  forces  simulate  on-board  vibrating  equipment. 

It  is  assumed  that  both  forces  will  act  simultaneously. 

The  tolerances  for  system  line  of  sight  error  and  defocus  are 


LOS  = 

1.0  x  10"6 

radians 

Defocus  = 

.50  x  10"3 

meters 
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Table  9-7 


•  Mass  changes 


Node 

Original 

Mass 

Case  If 2  -  2A 

Case  //4  -  4A 

27 

375.0 

350.0 

375.0 

28 

375.0 

350.0 

350.0 

29 

375.0 

350.0 

350.0 

30 

375.0 

350.0 

350.0 

32 

500.0 

500.0 

550.0 

33 

500.0 

550.0 

450.0 

34 

250.0 

300.0 

300.0 

35 

250.0 

250.0 

200.0 

•  Stiffness  changes 


Table  9-9 


Figure  9-8 


9.6.2  Open  Loop  Performance 

To  demonstrate  that  structural  control  is  required,  the  open  loop 
performance  Is  shown  In  Figure  9-9.  In  this  analysis  each  mode  is  assumed 
to  have  0.1%  critical  damping  and  there  is  no  additional  isolator  damping. 
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SECTION  10 


PRELIMINARY  ASSESSMENT  OF  DRAPER  MODEL 


10. 1  Introduction 

The  Draper  model  described  in  Section  9  was  analyzed  with  respect  to 
the  two  given  sinusoidal  disturbance  functions.  The  steady  state  line-of- 
slght,  LOS,  amplitude  error  is  determined  and  then  compared  to  the  LOS 
performance  requirements.  The  amount  of  modal  damping  required  to  meet  the 
LOS  performance  is  determined.  The  preliminary  results  show  that  modal 
damping  alone  is  not  sufficient  to  suppress  the  sinusoidal  disturbance 
impact  on  the  LOS  error. 


10. 2  Steady  State  LOS  Error 

The  amplitude  of  the  steady  state  LOS  error  is  determined  by  eval¬ 
uating  the  modal  response  at  t  =  «>.  The  modal  equation  is 

Hi  +  2^iwir*i  +  = 

where 

r)^  is  the  modal  amplitude  for  the  ith  mode 

is  the  modal  frequency  of  the  ith  mode 

is  the  damping  ratio  of  the  ith  mode 

<p  is  the  disturbance  influence  at  location 
’  ' d'  to  the  ith  mode 

f  is  the  disturbance  forcing  function 


For  this  particular  case,  the  disturbance  forcing  function  is  assumed 
sinusoidal . 


A  sin  lit 


t  *  0 


Evaluation  of  the  modal  response  at  t  “  00  yields  the  steady-state  solution 

A4>, 


\/K“ia)2  ♦  K  -  °2)! 


sin  (Qt  +  6) 


where 


3  -  tan 


■it  2W_\ 

R-W 
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The  steady  state  amplitude  of  the  modal  response  is  defined  as 


The  performance  requirement  for  the  LOS  error  is 

|  LOS  error  about  X-axis 

I  LOS  error  about  Y-axis 
defocus  error  in  Z-direction 

The  asterisk  in  Table  10-1  shows  that  19  of  50  modes  have  a  steady-state  LOS 
amplitude  greater  than  the  performance  requirement.  Note  that  modes  1  to  6 
are  the  rigid  body  modes.  The  defocus  error  always  meets  the  performance. 

The  definitions  for  the  headings  in  Table  10-1  are  as  follows: 

LOSXj  =  LOS  error  about  X-axis  due  to  disturbance  f.  ; 

j  =  1,2  (see  equation  (10-1,2))  ^ 

LOSYj  =  LOS  error  about  Y-axis  due  to  disturbance  f .  ; 

j  =  1,2  (see  equation  (10-1,2))  J 


Table  10-1.  Steady-state  line-of-slght  amplitude  as  functions  of 
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DEFOCUS j 


=  LOS  defocus  error  in  Z-direction  due  to 


LOS. 

J 


LOS 


disturbance  f 


j  ’ 


~\j  (LOSXj  ) 2  +  (LOSYj  )  2 


j  =  1,2  (see  equation  (10-1,2)) 

j  =  1,2 


V(L0S1)2  +  (LOS  2)  ^ 


10.4  Modal  Damping  Requirements 

The  damping  requirements  for  each  mode  are  determined  by  solving 
equations  (10-1),  (10-2)  for  and  specifying  the  desired  performance 
amplitude,  e  Qc,. 


__1 _ 

2w.fi 

l 


i,d^i,L0S 

£L0S 


Table  10-2  shows  the  results  for  the  disturbance  f^  at  node  #37  while 

Table  10-3  is  for  disturbance  f^  at  node  #46.  Ignoring  the  rigid  body, 

modes  1  to  6,  it  can  be  seen  that  there  are  numerous  damping  ratios  »  1. 

In  order  to  obtain  some  insight,  a  frequency  response  was  determined  (Fig.  10-1) 
between  disturbance  f  and  the  L0SX.  Note  the  solid  line  at  62.8  rad/sec 
(10  Hz).  Since  this  line  intersects  the  frequency  response  at  a  location 
removed  from  the  peaks,  just  adding  damping  will  only  affect  the  magnitude 
of  the  response  peaks  and  not  the  general  shape  of  the  frequency  response 
in  the  neighborhood  of  the  disturbance  frequency. 

10. 5  Summary 

Preliminary  results  show  that  the  defocus  error  is  not  degraded  by 
the  two  sinusoidal  disturbances.  Ignoring  the  rigid  body,  17  modes  degrade 
the  LOS  error  about  the  X  and  Y  axes.  Furthermore,  the  addition  of  damping 
to  individual  modes  is  not  sufficient  to  suppress  a  sinusoidal  disturbance 
impact  on  LOS  error. 
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Table  10-2 


Damping  ratio  required  to  meet  LOS  performance  (f^) 


DAMPING  RATIO 


^?0E 

FREQ 

AMPl 

L05X 

LOSY 

OEFOCUS 

X 

Y 

(r/s) 

(r> 

<r) 

(ml 

1 

0.000382 

-.67080-04 

- . 13150-03 

-.25340-07 

-.62470-13 

c 

0.CC04C9 

.52710-03 

-.33230-03 

.4442D-C9 

-.93190-11 

3 

0.000420 

-.55620-04 

-.29320-07 

.18680-08 

.57260-13 

4 

0.000436 

-.41970-04 

-.439^0-08 

-.43480-08 

-.18320-14 

5 

0.000460 

.17750-03 

-.62173-03 

-.14403-06 

-.43310-12 

6 

0.0C0472 

.13493-03 

- . 1C570-C6 

.63090-03 

. *54240-13 

7 

0.914426 

-.76720-04 

.20260-10 

-.25030-07 

.37*30-11 

8 

1.653811 

.40860-04 

.14260-10 

- .  e75CD-03 

-.44910-10 

9 

1.593596 

-.26220-03 

.39470-09 

.90560-10 

.13740-03 

10 

2.091*05 

. 52410-04 

-.12020-07 

.23100-10 

-.76680-11 

11 

2.7S4732 

.39970-03 

. 2415D-03 

.42*10-09 

.62300-03 

T  *> 

3.630322 

-.15740-03 

-.11520-06 

.10653-03 

-.30000-09 

13 

3.653243 

-. 17993-03 

-.41910-09 

-.15710-06 

. 76920-10 

14 

7.639052 

.71230-05 

-.73120-11 

-.16330-03 

- . 3441D-12 

15 

6.165607 

-.1O020-C4 

- .60240-10 

.24340-08 

.  12s4D-09 

16 

3.466574 

.67303-04 

-.10770-07 

-.95730-10 

.63490-09 

17 

10.S13C30 

.43660-04 

. 11510-03 

-.12133-06 

.21630-03 

IS 

11.427218 

.40210-05 

.10410-10 

-.20010-11 

.  10590-10 

19 

11.427431 

-.33220-05 

- . 3326D-10 

.73563-11 

.18250-11 

CO 

11 . 669540 

.11603-13 

-.41730-27 

-.40510-23 

.35140-23 

Cl 

14.850258 

-.72330-03 

-.16130-07 

- . 994  *0-05 

-.65300-03 

19.77 

-  ty 

18.733630 

—  •  «.  4  300  —  :  ^ 

-.32700-03 

-.66213-10 

.36680-09 

C3 

19.555C63 

.18600-04 

-.81020-10 

- .45410-03 

.465i0-10 

C4 

21.232741 

.25500-04 

.14420-10 

-.11753-07 

-.11020-10 

C5 

32.431671 

.15933-14 

- .25910-23 

-.70073-30 

.4144D-29 

Co 

33.051529 

.28600-05 

-.97100-10 

-.200*0-11 

.13570-10 

27 

49.4C2C33 

-.63520-16 

-.71100-31 

-.49720-32 

.65380-32 

C8 

50.995269 

. 7652D-03 

- . 3271D -06 

-.13010-05 

-.14020-06 

.1750 

C9 

52.527542 

-.13510-04 

-.39393-03 

-.64500-11 

.49130-10 

30 

53.e:csi5 

-.65330-03 

-.12310-04 

.73010-08 

-.70620-06 

1.899 

31 

55.3766C2 

. 74620-C6 

-.26523-10 

-.21720-12 

-.28020-11 

32 

55.376617 

.28410-05 

-.12770-09 

.22110-12 

.27130-11 

33 

71.2*0451 

-.16970-05 

-.63260-11 

-.15270-10 

.38540-12 

34 

72 . 242cu6 

.40493-03 

-.49600-06 

.75720-06 

.16120-06 

35 

79.555420 

-.69410-03 

- . 7615D-C5 

.39300-06 

.11C6D-04 

1.838 

3  3 

85.345627 

-.46400-04 

-.97330-09 

- . 34*90-07 

.31850-03 

37 

85.16S513 

.62470-03 

-.44630-06 

-.63170-05 

.71360-06 

2.003 

33 

83. 972637 

- . 9951D-04 

- . 393hD-C6 

- . 53050  -  06 

.107^0-06 

39 

93.346039 

-.80270-04 

.11350-05 

-.25790-06 

.39860-06 

.2990 

40 

100.535977 

- . 65340-04 

-.12460-06 

- . 63243- Co 

.14990-07 

41 

103.320430 

. 1307D-03 

-.82603-06 

-.29000-06 

.18460-05 

42 

1C5.214C50 

-.24333-03 

-.93753-06 

.23210-06 

. 33450-05 

*43 

107.791077 

-.20660-03 

-.23390-07 

.17360-05 

. 3592D-06 

.8038 

s4 

112.016337 

.17540-03 

.56900-06 

-.2C49D-06 

.12520-05 

45 

119.828273 

.12330-03. 

-.626CD-07 

-.16550-05 

-.1OS3D-06 

1.080 

46 

149.3-51542 

- .23160-05 

.52493-07 

- . 38S5D-07 

- . 1208D-C6 

47 

153.397614 

.20210-04 

. 1331D-G6 

.75010-06 

-.46350-06 

43 

162. 75G466 

-.28870-04 

- . 7762D-08 

- . 34300-06 

.95790-07 

49 

165.640686 

-.22523-03 

.32520-05 

.23480-05 

-.53680-06 

3.992 

2.398 

50 

166.055479 

.10660-15 

. 1CC40-29 

-.66720-30 

-.15840-30 

Table  10-3.  Damping  ratio  required  to  meet  LOS  performance  (fj). 


DAMPING  RATIO 


MODE 

FREQ 

AttPL 

10SX 

LOST 

DEFOCUS 

X 

Y 

!r  /«) 

Ir) 

(r) 

(m) 

1 

0.000332 

- . 5566D-03 

-.15110-07 

-.21100-06 

-.52020-12 

2 

0.000409 

.42300-02 

-.66830-07 

. 35640-08 

-.7C79D-10 

3 

0 . 00C420 

-.87630-03 

-.43180-06 

. 2  746D- 07 

.  64170-12 

4 

0.000436 

-.43030-03 

-.50170-07 

-.49700-07 

-.18780-13 

5 

0.000460 

.14490-02 

-.50730-07 

-.11750-05 

39420-11 

.21050 

6 

0.000472 

.17140-02 

-.13430-05 

. 6016D-07 

.11980-11 

.29330.05 

7 

0.914426 

.25060-02 

-.66180-09 

.81750-06 

- .  1229D-09 

8 

1.653211 

-.33470-04 

-.13430-10 

.32390-03 

.42290-10 

9 

1.993546 

-.30250-03 

.45580-09 

.10460-09 

. 15670-03 

10 

2.091605 

.57210-03 

-.13120-06 

.25210-09 

-.63910-10 

11 

2.784732 

-.20560-02 

-.12420-07 

-.21920-03 

-.32050-07 

12 

3.630322 

-.175«D-02 

-.12940-05 

.11900-07 

-.33530-0S 

3.508 

13 

3.653243 

-.21730-03 

-.50630-09 

-.18970-06 

.92930-10 

14 

7.6S9C52 

.95470-02 

-.96050-08 

-.22550-05 

-.46120-09 

3.883 

15 

8.169607 

.11330-03 

.42600-09 

-.17210-07 

-.67970-09 

16 

6.466S74 

-.22020-01 

.35000-05 

.31120-07 

- . 2062D-06 

5.771 

17 

10.813030 

-.4S153-06 

-.13200-10 

.13380-08 

.23620-10 

18 

11.427218 

-.14400-02 

-.37310-03 

.71660-09 

-.37940-08 

19 

11.427431 

•476&0-03 

.54880-08 

-.11270-03 

-.26180-09 

20 

11. 86 9«40 

-.80460-11 

.28050-24 

.28310-25 

-.24360-25 

21 

14.850258 

-.56653-03 

-.12630-07 

-.7752D-05 

-.51140-08 

6.348 

22 

18.733630 

.54340-01 

.73110-05 

.14310-06 

-.S2C2D-C6 

3.891 

23 

19.955063 

.11280-03 

-.49110-09 

-.29340-07 

.281^0-09 

24 

21.2S2791 

- . 35430-01 

-.20030-07 

.16300-04 

. 153C0-07 

6.497 

25 

32.431671 

-.75850-10 

.13760-23 

. 3336D-25 

-.19730-24 

26 

33.051529 

-.84530-01 

.28700-05 

.59290-07 

- .4012D-06 

.1366 

27 

49.492233 

.10640-12 

.11040-27 

.77260-29 

-.10230-28 

23 

50.959C69 

-.45130-05 

.19290-08 

.76750-08 

.6SCCD-09 

29 

52.527542 

.83250-02 

.25740-05 

.42140-03 

-.32100-07 

1.273 

30 

53.850S15 

-.14520-03 

-.27120-05 

.16050-08 

- . 1556D-C6 

1.425 

31 

55.376  6  0  2 

.20360-03 

-.79730-03 

-.608=0-10 

-.78340-09 

32 

55.376617 

-.59850-03 

.26910-07 

-.4659D-10 

- .57170-09 

33 

71.290451 

-.31230-07 

-.11750-12 

-.28090-12 

.16360-13 

34 

72 . 242966 

. 794^0-05 

-.97330-08 

.14860-07 

.31630-03 

35 

79.953420 

-.33540-04 

-.36810-26 

.16590-07 

.53440-06 

36 

85.345627 

.43780-02 

.91820-07 

.32720-05 

-.30050-06 

3.659 

37 

86.168518 

.31040-03 

-.22170-06 

-.31390-05 

.35460-06 

3.538 

38 

68.972637 

.28010-04 

.11080-06 

. 1663D-C6 

-.30380-07 

39 

93. 346039 

.12270-04 

-.17360-06 

.39430-07 

-.60950-07 

40 

100.985977 

-.10640-04 

-.20260-07 

-.11110-06 

.24400-03 

41 

103.825430 

-.53260-05 

.36110-07 

.11820-07 

-.75230-07 

,42 

105.214050 

.43300-05 

.16660-07 

-.4130D-03 

-.59520-07 

43 

107.791077 

. 7116D-05 

.60590-09 

-.59830-07 

-.12360-07 

44 

112 . 0183S7 

-.4S43D-05 

-.15720-07 

.56620-08 

- . 3477D-07 

45 

119.826273 

-.20770-05 

.13370-03 

.30030-07 

.32100-08 

46 

149.361542 

.29780-06 

-.67460-03 

.49700-08 

.15540-07 

47 

153.397614 

-.11850-05 

-.10740-07 

-.43=90-07 

.27160-07 

48 

162.790466 

-.26850-06 

-.72160-10 

-. 3236D-08 

.8907D-09 

49 

165.640636 

-.50270-06 

.72580-08 

■524CD-03 

- .  1197D-03 

50 

166.059479 

- . 12o50-16 

-.11910-30 

. 7916D-31 

. 1879D-31 
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